F/G  9/5 


AD-A036  321  NAVY  UNDERWATER  SOUND  LAB  NEW  LONDON  CONN 

MATCHED  FILTERING  AND  OPTIMAL  USE  OF  AN  ANTENNA# (U) 

JUN  66  H MERMOZ 

UNCLASSIFIED  USL-PUB-745  NL 


I cf2 

«DA0363^l 


ftp  HOST  EroJect 


Matched  Filtering 


and  Optimal  Use  of  an  Antenna 


HENRyMERMOZ  j 

P oly technique,  0 renoble 

'ran slated  from 


Institut 


“^Signal  Processing  — with  Emphasis  on  Underwater  Acoustics 


NATO  Advanced  Study  Institute, 
f 26  Sep«e«rt»ri4*64,  fpi  162-299- 

1 SV, 

7 by  / 

| RoBERf/O’ Neill,  Jr,  l 


i£! j^Fi(l±(l>  3 It  ^ 


U.  S.  Navy  Underwater  Sound  Laboratory 
Fort  Trumbull,  New  London,  Connecticut 


DISTRIBUTION  STATEMENT  A 


Approved  for  public  release; 
l‘  Distribution  Unlimited 


ABSTRACT 

Originally  published  in  the  French  language  in  "Signal 
Processing  — with  Emphasis  on  Underwater  Acoustics," 
by  the  NATO  Advanced  Study  Institute,  14-26  September 
1964,  pp.  162-299,  this  study  has  been  translated  to 
serve  as  a ready  reference  work.  The  study  itself  lies  in 
the  area  of  research  that  is  concerned  with  the  improve- 
ment of  long-range -detection  techniques.  Subjects  covered 
include  (1)  a formal  solution  to  matched  filtering  with  N 
inputs,  (2)  the  case  of  identical  signals — proper  filter- 
ing— orthogonal  images  of  a system  with  N noise  inputs, 
(3)  narrow-band  approximation  — variation  of  proper 
filtering,  (4)  case  of  two  inputs  — coherent  noise,  (5) 
matched  filtering  and  directivity,  and  (6)  generation  of 
proper  filtering  — autoadaptive  systems  — practical 
stationarity  conditions . 


ADMINISTRATIVE  INFORMATION 

The  translation  presented  in  this  document  was  per- 
formed for  USL  Problem  No.  1-410-00-00,  Navy  Sub- 
project  and  Task  No.  SF  101  03  1^-11224. 


REVIEWED  AND  APPROVED:  29  June  1966 


F.  H.  Hunt 

Assistant  Technical  Director  for  Administration 


TABLE  OF  CONTENTS 


Page 


GLOSSARY  OF  SYMBOLS  iv 

CHAPTER  I.  SOME  PRACTICAL  PROBLEMS  AT  THE  BEGINNING 

OF  THE  STUDY.  1 


Summary.  I-I.  General  Outline.  1-2.  Matched  Filtering  and 
Multiple  Antennas.  1-3.  Some  Comments  on  Notation 

CHAPTER  II.  RESTATEMENT  OF  THE  CLASSICAL  MATCHED 

FILTER  THEOREM 

Summary.  II-l.  History  of  the  Concept  of  Matched  Filtering. 

II -2 . Matched  Filtering  with  One  Input. 

CHAPTER  III.  MATCHED  FILTERING  WITH  N INPUTS  - 

FORMAL  SOLUTION 

Summary,  lll-l.  Status  of  the  Problem.  I1I-2.  Review  of  the 
Classical  Properties  of  the  Crosscorrelation  of  Two  Noise  Inputs. 

LII-3.  Review  and  Physical  Interpretation  of  a Fundamental 
Property  of  the  Crosscorrelation  of  Two  Noises.  UI-4.  Investi- 
gation of  the  Solution  to  the  Problem  Posed  in  Paragraph  HI-1. 
ni-5.  Demonstration  - Temporal  Form  of  the  Solution. 

IH-fi.  Spectral  Form  of  the  Solution  - Value  of  the  Optimized 
Parameter.  I1I-7.  Matrix  Notation  - Characteristic  Property  of 
Matched  Filtering.  III-8.  Properties  of  the  Output  Signal. 

Appendix  I.  Appendix  II. 

CHAPTER  IV.  REMARKS  CONCERNING  THE  FORMAL  SOLUTIONS 

OF  SOME  SPECIAL  CASES 

Summary.  IV-1.  Limitations  of  the  Cases  Considered  in  this 
Chapter.  IV-2.  The  Case  of  a Single  Input.  IV -3.  Matched 
Filter  Rule  for  Uncorrelated  Noise.  Signal  Energy  and  Signal- 
to-Noise  Ratio.  IV -4.  General  Case  - Usefulness  of  the  Arbitrary 
Parameter  t IV-5.  Limited  Possibilities  for  Elimination  of 
Noise  Having  a Particular  Structure. 

CHAPTER  V.  THE  CASE  OF  IDENTICAL  SIGNALS  - PROPER  FILTERING  - 
ORTHOGONAL  IMAGES  OF  A SYSTEM  WITH  N NOISE  INPUTS 

Summary.  V-l.  Nonrestrictive  Character  of  the  Case  of  Identical 
Signals.  V-2.  Formal  Solution  for  the  Case  of  Identical  Signals. 

V-3.  Proper  Filtering  of  an  Antenna  and  its  Gain.  First  Charac- 
teristic Property.  V-4.  Review  of  the  Known  Properties  of  the 
Correlation  Matrix.  V-5.  Second  Characteristic  Property  of 
Proper  Filtering. 

DISTRIBUTION  stTiemcnt  a 

1 Approved  lor  public  release; 

Dirtiibution  Unlimited 


TABLE  OF  CONTENTS  (Cont'd) 


d 


Page 

CHAPTER  VI.  NARROW-BAND  APPROXIMATION  - VARIATION 

OF  PROPER  FILTERING  (NPF) 51 

Summary.  VI-1.  Narrow-Band  Hypothesis.  Physical  Meaning 
of  the  C^(vD)  and  of  the  Matrix  Notation.  VI-2.  Normalized 
Proper  Filtering  - Technical  Advantages.  VI-3.  Invariance  of 
the  NPF  in  the  Presence  of  the  Same  Noise  Superimposed  on 
Each  of  the  :,.(t)  . VI-4.  Matrix  Notation  for  NPF.  VI-5. 

Limiting  Behavior  of  the  NPF.  Eliminable  Noise  and  Cutoff. 

VI-(j.  Limiting  Case  of  a Nulled  Input  Noise.  VI-7.  Influence 
of  Independent  Noises  Superimposed  on  the  N Inputs. 

CHAPTER  VII.  TWO-INPUT  CASE  - COHERENT  NOISE f>3 


Summary.  VII-1.  General  Equations.  VII-2.  Eliminable 
Interference  and  Cutoff.  VII-3.  Advantages  of  Matched  Filtering 
over  the  Direct  Sum.  VII-4.  A System  of  Graphs  Representing 
the  Gain  of  the  NPF.  VII-5.  Coherent  Interference  (Jamming)  - 
Mechanism  for  its  Elimination.  VII-(>.  Return  to  a Remark  Made 
in  Chapter  I - The  Special  Case  of  Coherent  Interference.  VII -7. 

Optimum  Use  of  an  Antenna. 

CHAPTER  VIII.  MATCHED  FILTERING  AND  DIRECTIVITY 80 

Summary.  VII-1.  Limits  of  the  Validity  of  the  Directivity  Concept. 

VIII-2.  Definition  and  Physical  Meaning  of  Directivity.  VIII-3. 

Expression  for  the  "Directivity  Factor"  as  a Function  of  the  Gain 
of  the  Multiple  Filters.  VIII-4.  Optimization  of  the  Directivity 
Factor.  V1II-5.  Correlation  Matrix  for  Omnidirectional  Noise. 

VIll-(i.  Matched  Filtering  in  Omnidirectional  Noise  - Similarity 
to  Optimization  of  Directivity  Factor.  VIII-7.  Application  to 
Two  and  Three  Inputs. 

CHAPTER  IX.  GENERATION  OF  PROPER  FILTERING  - AUTOADAPTIVE 

SYSTEMS  - PRACTICAL  ST  AT  ION  A R IT  Y CONDITIONS 94 

Summary.  IX-1.  Representation  of  a Crosscor relation  (Narrow  Band) 
of  Two  Slowly  Varying  Voltages.  EX-2.  Practical  Stationaritv  ;uid  Time 
Constants  in  an  Autoadaptive  System.  EX-3.  Representation  of  a Cross- 
correlation  by  ;in  AC  Voltage.  IX-4.  Principles  of  Producing  a Product 
of  Complex  Numbers  in  the  Form  of  an  AC  Voltage.  IX-5.  Narrow-Band 
Linear  Filtering  Using  Multiplication  by  a Voltage  Representing  the 
Filter.  IX-<>.  Flow  Charts  for  Two-and  Three-Input  Systems.  IX-7. 
Verification  of  One  Property  of  the  NPF  in  the  Preceding  Diagrams. 


ii 


<_ 


TABLE  )i  CONTENTS  (Cont'd) 


Page 


CONCLUSIONS 

REFERENCES  

INITIAL  DISTRIBUTION  LIST 


Ill 

113 

Inside  Back  Cover 


D D C 

:-.:^rr.nrLQG 

MAR  3 I9T7 


GLOSSARY  OF  SYMBOLS 


a 

a(  t ) 

B(t),b(t) 

b 

i w 

bj 

e(t), eD(t) 

cjk<u) 

Cn(T> 

Cn(v) 

c 

c 

r 

p 

,jk(,) 

Vv) 

Yr  (v  ) 

1 0 

U(  w,  v ) 

dI 

A 

6 ( t ) 

E x 
E 


J ’ k 


E. 


P(a) 


separation  of  receiver  elements 
column  matrix  of  ones 
unitary  signal 
noise  inputs 

ratio  of  S/i\(NPF)  to:3/N  (simple  summation) 

. th 

elementary  noise  component  from  direction  w at  .1  receiving 
element 

noise  at  filter  output 

, th  th 

crosscorrelation  oi  J and  k noise  components 
cross  power  spectrum  of  i L 1 and  kL1  noise  components 
autocorrelation  of  noise  n ( t ) 
power  density  spectrum  of  noise  n(t) 

matrix  of  C . , 

Jk 

velocity  of  propagation 
normalized  proper  filter  gain 
crosscorrelation  of  filtered  noise  components 
cross  power  spectrum  of  filtered  noise  components 

spectral  density  of  noise  at  output  ^ 
complex  directivity  coefficient 
spectral  density  of  summed  noises 
determinant  of  matrix  of  C . ^ 
unit  impulse  function 

mathematical  expectation  of  x 

. t h i t h \ 

receiver  elements  (.  j > K > 

signal  energy 

frequency  (cps)  of  reference  voltage 

spatial  correlation  of  noise  inputs  to  receiver  elements 
with  separation  a 


iv 


hk(v) 

hk(t) 


M^(.) 

Mj^v) 


NJk(t) 

N,JkCv) 


pk(v) 


frequency  (eps) 

system  function  of  proper  filter 

th 

system  function  of  k filter 
impulse  response  of  k^’1  filter 
matrix  of  hk(v) 

imaginary  part  of  complex  number  x 

& 

complex  correlation  coefficient 
time  constant 

average  duration  of  stationarity 
(subscripts)  indices 
constants 

eigenvalue  of  matrix  C 
matrix  of  eigenvalues  X . 
cofactor  of  determinant  A 
matrix  of  cofactors  M . , . 

crosscorrelation  of  noise  components  from  direction  to  at 
,•  th  and  kth  receiving  elements 

CO 

cross  power  spectrum  derived  from  Mj  k ( t ) 
constant 

number  of  elements  or  channels 
crosscorrelation  of  and  kltl  noise  components 
cross  power  spectrum  of  ,1  ^ ' and  noise  components 

arbitrary  noise  input 
frequency  (cps) 
mean  noise  power 

noise  power  minimized  by  optimal  filtering 
system  function  of  filter 


v 


p 


"k 

Q 


q 

R ( t ) , r . ( t ) 
J J 

V'0-  rj<u) 

R,  r 
R(  x ) 
p 


S(t) 


a ( v ) 
T 


i 


x 


(v) 


X 

V 

*(  v) 

4> 


»j(.) 

X 


matrix  of  pk 

th 

k proper  t liter 
homogeneity  factor  (constant) 

V°  11 

- — for  2 -input  system 

L 12 

impulse  response  of  filters 
system  functions  of  filters 
matrices  of  R ,,  fj 
real  part  of  complex  number  x 
signal-to-noise  ratio  (S/N) 

S/M  for  simple  summation 
optimized  (maximized)  S/N 
signal 

Fourier  transform  of  output  signal 
time  duration  of  signal 
time 
time  lag 

. . th  th 

normalized  cross  power  spectrum  of  j and  k noise 
inputs  after  signal  compensation  - case  of  omnidirectional 

noise 

matrix  of  V* 

J k 

system  function  of  filter  with  unit  gain 
matrix  of  $ ( v ) 

bandwidth  of  narrow-band  filter 

th 

system  function  of  J filter 
column  matrix  of  * 


vi 


VT) 


Vv) 


w 


crosscorrelation  of  noises  at  and  1 receiving 
elements  in  omnidirectional  noise 

cross  power  spectrum  derived  from  y ( t ) 

spatial  parameter  defining  direction  of  a plane  wave 

unit  vector  in  direction  w 

numbers  in  parentheses  correspond  to  references 


CHAPTER  I. 


SOME  PRACTICAL  PROBLEMS  AT  THE  BEGINNING  OF  THE  STUDY 


This  study  lies  in  that  area  of  research  that  is  concerned  with  the  improvement 
of  long-range  detection  techniques.  By  extending  the  well-known  concept  of  a 
matched  filter  to  include  a receiver  with  several  elements  and  by  using  a few 
simple  examples,  it  is  shown  that  current  ideas  about  the  optimum  use  of  an  antenna 
provide  answers  to  only  a few  special  aspects  of  the  general  problem. 

1-1.  General  Outline 

M’he  problem  of  extracting  weak  signals  from  interfering  noise  may  be  found  in 
the  general  area  of  Information  Theory.  ,We  may  consider  two  different  approaches 
to  the  problem,  each  arising  from  a different  practical  situation.  The  first  approach 
is  concerned  with  the  faithful  reproduction m, an  input  signal  at  the  output  of  a com- 
munication channel.  Here,  the  applicable  thfeory  is  that  of  information  transmission; 
the  signal  is  unspecified  except  for  the  variety  of  forms  or  values  which  it  may 
assume  and  the  probability  of  occurrence  of  each.  Optimization  criteria  are  based 
upon  the  reconstruction  of  the  most  probable  form  of  the  signal,  that  is,  the  best 
"guess"  at  the  original  form  and  content  of  the  transmitted  signal.  This  first  ap- 
proach is  that  of  optimizing  the  reconstruction  of  a message  that  has  been  deterio- 
rated by  noise. 

The  second  aspect  of  the  processing  of  weak  signals  - the  aspect  with  which  this 
paper  is  concerned  - has  its  origins  in  practical  problems  associated  with  the  develop- 
ment of  long-range  detection  methods,  using  both  electromagnetic  (radar)  and 
acoustic  (sonar)  techniques.  Here,  we  wish  to  be  informed  of  the  presence  of  an 
obstacle  or  target  as  soon  as  possible  and  at  as  great  a distance  as  possible.  The 
"echo"  returned  by  such  a target  or  obstacle  is  weak,  immersed  in  the  inevitable 
interfering  noise,  but  it  is  assumed  to  be  of  a specified  or  known  foi'm  - the  same 
form  as  a specific  model  (i.e.,  identical  to  the  model  except  for  some  constant  ampli- 
tude factor  and  time  lag).  The  model  might  be  the  original  emitted  signal,  which  may 
be  controlled  or  even  derived  from  it  by  a known  transformation  (e.g.  doppler  shift). 
Conceivably  each  possible  doppler  shift  value  may  yield  a different  "model."  We  may 
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group  the  variables  in  the  received  signal  in  order  to  form  a reasonable  number  of 
specific  models  and  treat  all  models  separately  but  simultaneously.  The  basic 
problem  is  reduced  to  the  question:  what  are  the  criteria  which  optimize  the  detec- 
tion of  the  presence  of  a signal  of  a given  form?  It  is  not  required  that  the  signal 
be  accurately  reproduced;  it  is  even  permissible  to  transform  and  alter  it.  What  is 
required  is  the  best  possible  indication  of  its  presence;  this  indication  might  be,  for 
example,  a maximum  contrast  in  the  output  level  of  the  receiver. 

1-2.  Matched  Filtering  and  Multiple  Antennas 

A particular,  but  essential,  solution  to  the  problem  stated  in  the  preceding 
paragraph  may  be  found  using  the  Matched  Filter  Theorem.  This  study  deals  with 
the  scope  of  that  theorem  and  with  the  criteria  which  it  utilizes. 

In  its  classical  form,  the  Matched  Filter  Theorem  treats  a single  signal  and  a 
single  noise  and  does  not  deal  with  the  multi-parameter  aspects  of  the  detection  of 
weak  signals  using  an  antenna.  An  antenna  is  an  ensemble  of  N sensors  whose  out- 
puts are  collected  in  some  more  or  less  complex  fashion.  For  each  sensor  there  is 
a corresponding  noise.  When  a signal  reaches  an  antenna,  it  is  distributed  in  some 
specific  maimer  into  N particular  signal  outputs. 

What,  then,  is  the  best  possible  utilization  of  the  N signals  and  the  N noise 
outputs  for  optimal  detection  of  the  presence  of  an  incident  signal?  What  is  the 
best  arrangement  of  the  sensors?  And  what  gain  might  be  expected  from  using 
N sensors  as  compared  to  using  a single  sensor?  It  may  be  noted  that  the  usual 
ways  of  evaluating  antenna  gain  answer  the  last  question  in  special  cases  only.  In 
some  cases,  this  gain  is  expressed  simply  by  the  number  N itself. 
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Classical  reasoning  assumes  that  at  the  output  the  signal  power  increases  as 
while  that  of  the  noise  increases  as  N This  reasoning,  however,  postulates 
statistical  independence  among  the  noise  outputs  from  the  U elements  or  sensors 
and  is  therefore  dependent  upon  a particularly  simple  hypothesis  about  the  statis- 
tical "relationships"  of  the  N noise  outputs. 

In  another  case,  antenna  gain  is  computed  in  terms  of  "directivity;"  thereby 
introducing  the  more  or  less  implicit  hypothesis  of  a noise  field  with  particular 
spatial  properties,  known  as  "omnidirectional  noise,"  which  completely  defines  the 
statistical  relationships  of  the  N noise  components. 
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Two  methods  of  improving  an  antenna  are  linked  to  these  two  means  of  evaluating 
antenna  gain.  In  the  case  of  iue  -pendent  noise  components,  N is  made  very  large  and 
the  outputs  of  the  N elements  are  summed.  This  "direct  sum"  is  the  antenna 
processing,  that  is,  the  way  in  which  the  outputs  of  the  elements  are  assembled. 

In  the  case  of  omnidirectional  noise,  there  are  more  subtle  solutions,  consisting, 
for  example,  of  assigning  to  each  element  a suitable  weighting  factor  before  summing 
the  outputs,  such  that  the  "directivity"  is  optimized. 

Methods  for  the  calculation  of  directivity,  derivation  of  directivity  patterns  for 
specific  arrays  of  elements,  and  methods  for  reducing  or  equalizing  "secondary 
lobes"  have  provided  abundant  literature  (references  (6),  (15)  to  (22))  because  of  the 
usefulness  of  the  directivity  concept  in  both  transmission  and  reception.  The  under- 
lying hypothesis  of  omnidirectional  noise  is,  however,  always  present  when  this 
concept  is  applied  to  a receiving  antenna. 

The  variety  of  techniques  for  estimating  antenna  gain  suggests  that  optimal 
antenna  processing  depends  upon  the  statistical  relationships  between  the  noise 
components  associated  with  the  elements  of  the  antenna. 

For  a given  antenna,  the  noise  outputs  are  known.  Their  properties  are,  in 
general,  stable  (stationary),  at  least  for  a time  duration  on  the  order  of  the  signal 
duration.  Thus,  they  may  be  measured  and  the  best  methods  of  handling  them  may 
be  determined.  These  methods  may  include  direct  summation  or  optimization  of 
directivity,  if  the  nature  of  the  noise  corresponds  to  these  special  cases,  or  they 
may  be  something  quite  different  and  may  not  even  require  that  N be  large.  Such 
is  the  case,  for  instance,  with  "coherent"  noise  interference,  which,  like  the 
signal,  propagates  in  a single  plane  wave  (from  a different  direction  from  that  of 
the  signal).  This  coherent  interference  is  a case  of  still  another  statistical  rela- 
tionship: the  noise  components  in  the  antenna  elements  differ  from  one  another  by 
a time  delay  only.  Using  only  two  elements,  it  is  possible  to  null  the  noise  without 
losing  the  signal,  which  is  obviously  the  optimum  procedure.  All  that  is  required 
is  to  take  the  difference  between  the  outputs  of  the  two  elements  after  having 
introduced  into  one  of  them  a delay  corresponding  to  the  difference  in  path  length 
for  the  noise  itself.  A limited  case  such  as  this  illustrates  the  way  in  which 
optimum  antenna  processing  may  occasionally  be  achieved  when  the  statistical 
relationships  between  the  noise  components  are  taken  into  account.  Hence,  if  the 
N noise  components  are  known,  the  problem  is  to  define  and  construct  the  optimal 
processing  system. 
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A logical  second  step  is  to  allow  the  statistical  relationships  between  noise 
components  to  vary  slowly  with  time  in  an  irregular,  unpredictable  manner.  We 
must  then  look  for  auto  adaptive  systems  that  evolve  the  optimum  antenna  processing 
as  a function  of  the  statistical  relationships  - relationships  which  are  continually 
tested  in  the  course  of  that  evolution. 

A great  deal  of  work  to  date  has  been  devoted  to  autoadaptive  systems,  but  in  a 
rather  different  theoretical  context.  In  general,  it  has  been  concerned  with  the 
construction  of  receivers  that  progressively  "adjust"  themselves  to  the  carrier 
frequency  of  a repetitive  signal  (as  in  the  case  of  radar).  The  evolution  of  the  re- 
ceiver is  thus  dictated  by  the  signal  itself  (references  (25),  (26),  (27)). 

Here,  on  the  other  hand,  we  want  to  modify  the  receiver  as  a function  of  the 
noise,  which  is  always  present,  in  such  a way  that  reception  will  be  optimized  for 
the  signal  whenever  it  arrives.  Note  that  the  idea  of  adaptation  does  not  have  the 
same  meaning  in  a "matched  filter"  as  it  does  in  an  "adaptive  system".  In  the  first 
case,  it  has  a "spectral"  sense  and  corresponds  to  a well-defined  criterion;  in  the 
second,  it  takes  on  a "temporal"  meaning  and  is  found  to  be  tied  in  witii  a practical 
conception  of  stationarity. 

These  practical  problems  are  mentioned  here  because  they  will  be  considered  in 
the  course  of  the  first  few  chapters,  employing  more  theoretical  techniques.  They 
will  not  be  resolved  in  all  their  generality,  but  only  in  those  special  cases  which  are 
of  particular  interest.  Thus,  the  ensuing  study  is  developed  in  the  direction  of  de- 
cre-  sing  generality,  from  the  theoretical  basis  formed  by  generalizing  the  Matched 
Filter  Theorem  to  technological  methods  and  principles  useful  in  the  construction 
oi  equipment.  Meanwhile,  some  general  properties  of  the  statistical  relationships 
of  the  noise  components  and  their  illustration  through  the  important  concept  which  wc 
have  called  proper  antenna  filtering  will  be  examined,  especially  with  respect  to 
matched  filtering.  The  cases  in  which  it  is  possible  to  obtain  very  large  antenna  gain 
will  be  stated  precisely,  and  it  will  be  shown  that  optimization  of  "directivity"  is 
actually  a special  case  of  matched  filtering. 

1-3.  Some  Comments  on  Notation 

a.  As  a general  rule,  two  functions  l'elated  by  a Fourier  transformation  are 
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designated  by  the  same  letter,  for  example1 

h(t)  + h(  v) 

This  convention  does  not  lead  to  any  ambiguity,  since  the  variable  is  always  indicated. 
It  has  the  advantage  of  economy  of  symbols. 

b.  The  symbol  J indicates  an  integral  between  the  limits  -<*>  and  +°°  . 

c.  We  have  made  use  of  some  results  of  harmonic  analysis  of  unspecified 
functions  (references  (1)  and  (7)). 

Generally  the  variable  t is  indicated  for  auto-  and  crosscorrelation  functions. 

The  variable  v may  be  omitted  from  the  notation  for  the  Fourier  transforms  of  the 
correlation  functions,  especially  when  dealing  with  matrices.  The  complex  number 
designates  either  C^(v)  or  C.^(vo)  for  a particular  value  v Q of  v 
The  convention,  in  this  case,  is  always  defined  in  this  text. 
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RESTATEMENT  OF  THE  CLASSICAL  MATCHED  FILTER  THEOREM 


Summary 

This  theorem,  taken  from  the  area  of  prediction  theory,  is  concerned  with 
specific  signals  (of  known  form)  mixed  with  stationary  noise.  It  defines  the  linear 
filtering  process  that  optimizes  a certain  parameter  called  "signal -to -noise  ratio." 
This  parameter  summarizes  the  best  possible  information  about  the  presence  or 
absence  of  the  signal  without  being  concerned  with  preserving  its  "shape"  or  form 
at  the  output.  Although  this  "best  possible  information"  is  limited  to  the  case  of 
Gaussian  noise,  it  is  frequently  found  in  practice.  The  optimizations  derived  during 
the  course  of  the  discussion  are  in  agreement  with  matched  filter  theory. 

H-l.  History  of  the  Concept  of  Matched  Filtering 

The  concept  of  matched  filters  was  introduced  in  the  technical  literature  around 
1943  by  D.  O.  North  (10)  in  a study  of  the  detection  of  weak  signals  of  some  known 
form  S(  t ) in  noise  that  is  stationary  and  has  a uniform  spectral  density.  North's 
essential  result  is  that  the  filter  which  maximizes  at  its  output  a certain  parameter 
known  as  the  "signal-to-noise  ratio"  is  the  filter  whose  impulse  response  is  the 
image  of  the  signal:  S(-t ) . This  result  justifies  the  expression  "matched  filter" 

later  used  by  Van  Vleck  and  Middleton  (11),  who  obtained  the  same  result  inde- 
pendently. The  extension  of  this  theorem  to  the  case  of  noise  having  a nonuniform 
spectral  density  is  found  in  the  work  of  B.  M.  Dwork  (5),  as  well  as  in  that  of 
L.  A.  Zadeh  and  J.  R.  Ragazzini  (12). 

These  latter  authors,  in  particular,  demonstrated  that  the'signal -to -noise  ratio" 
criterion  is,  in  fact,  a criterion  of  the  "separation"  between  the  signal  and  the  noise. 
This  observation  connects  matched  filtering  with  the  Wicner-IIopf  equation  (8). 

Criteria  of  "separation'/  however,  are  non-statistical  criteria,  which  depend 
only  upon  the  spectral  density  or  autocorrelation  of  the  noise,  that  is,  on  a single 
moment  of  the  probability  distribution.  The  noise  is  entirely  specified  by  this  moment 
only  in  the  case  of  a Gaussian  distribution. 

Another  approach  to  the  problem  of  optimizing  weak  signal  detection  is  through 
statistical  decision  theory.  The  literature  on  the  subject  is  abundant,  and  we  intend 
here  only  to  skim  its  surface  in  order  to  specify  the  criteria  for  a matched  filter 
(references  (3),  (4),  (9),  (28),  (29),  (30),  (31)).  This  theory  takes  into  account  the 
ensemble  of  statistical  properties  of  the  combination  of  signal  and  noise  in  order  to 
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deal  with,  especially,  the  a posteriori  probabilities  (for  a given  combination)  of  the 
presence  or  absence  of  signal.  The  ideal  receiver  must,  then,  construct  the 
"probability  ratio"  connected  with  these  probabilities.  It  may  be  shown  that  if  the 
noise  is  Gaussian  and  has  a uniform  density,  this  probability  ratio  is  completely 
described  by  the  convolution  of  the  "mixture"  (of  signal  and  noise)  with  the  "pure" 
signal.  This  convolution  is  nothing  more  than  the  output  of  the  matched  filter 
(having  impulse  response  S(-t)  ).  Thus,  matched  filtering  appears  as  a "statistical 

optimum"  in  the  Gaussian  case.  This  case  is  especially  important  in  practice.  The 
signals  to  be  processed  are  often  "narrow  band",  that  is,  lying  in  a small  band 
around  a center  frequency.  It  may  then  be  assumed  that  the  spectral  density  is 
uniform  in  the  band. 

On  the  other  hand  (reference  (13)),  for  a very  large  class  of  non -Gaussian  noise, 
narrow-band  filtering  tends  to  restore  some  of  the  Gaussian  characteristics.  It 
may  be  said  that  optimization  in  the  sense  of  matched  filter  criteria  is  often  very 
close  to  optimization  in  the  statistical  sense.  One  may  then  think  of  matched  filtering 
as  being  very  close  to  an  ideal  receiver. 

These  considerations  underline  the  importance  of  the  matched  filter  criterion  that 
may  at  first  seem  a little  arbitrary.  This  criterion  is  the  maximizing  of  the  signal- 
to-noise  ratio  defined  in  the  following  way: 

p = instantaneous  power  of  the  signal  at  arbitrary  time  t 

average  power  of  the  noise 

As  stated  more  precisely  in  the  succeeding  pages,  this  criterion  corresponds  to  the 
requirement  of  producing  at  the  filter  output  maximum  contrast  between  the  presence 
and  the  absence  of  a signal.  This  is  to  be  accomplished  by  collecting  all  the  energy 
of  the  signal  in  order  to  produce  a peak  that  is  as  narrow  and  as  high  above  the 
average  noise  power  as  possible. 

II-2.  Matched  Filtering  With  One  Input 

Let  us  recall  that  the  theoretical  solution  of  the  matched  filter  may  be  presented 
in  two  equivalent  forms  related  by  a Fourier  transformation:  the  temporal  form, 
which  furnishes  the  impulse  response  h(t)  of  the  filter;  and  the  spectral  form,  which 
gives  its  transfer  function  h(v) 

h(t)  i h(v)  . (TI-1 ) 


7 


Let  S ( t ) be  the  signal  having  a spectrum  S ( v ) , let  B ( t ) be  the  stationary  noise  whose 
autocorrelation  is  C ( t ) derived  from  its  spectral  density  C ( v ),  and  let  1 0 be  an 
arbitary  time. 


a.  The  impulse  response  h(t)  is  specified  by  the  integral  equation 


/h(e)  • c(t-e)  de  = ks( t o -t ) , 

where  K is  an  arbitrary  real  factor. 

b.  The  transfer  function  h(  v ) is  given  by 


(H-2) 


-2irlvt  o 

h(v)  = Ke  . S* ( v ) , 


( H-3 ) 


where  * designates  the  complex  conjugate. 

The  factor  e °ofEq.(II-3)  corresponds  to  a time  lag  of  t 0 • 

Since  t 0 and  K are  arbitrary,  it  may  be  said  that  the  matched  filter  is  defined 
except  for  a real  factor  and  a time  lag.  This  is  physically  obvious.  The  real  factor 
acts  in  the  same  way  upon  both  signal  and  noise  and  can  not  alter  their  behavior. 

As  for  the  time  lag  t 0 , if  it  varies,  the  time  of  the  appearance  of  the  signal  "peak" 
is  more  or  less  displaced,  but  its  height  is  not  changed.  It  is  preferable  generally 
to  have  the  time  lag  as  small  as  possible  in  order  that  the  observer  may  be  informed 
of  the  presence  of  the  signal  as  soon  as  possible  after  it  reaches  the  receiver  input. 
However,  it  may  happen  that  considerations  relating  to  the  realization  of  the  filter 
will  lead  to  a compromise  in  this  area. 

Nevertheless,  we  will,  in  the  course  of  this  dissertation,  again  encounter  the 
factor  :.e  V which,  as  we  shall  see,  has  no  effect  upon  the  filtering  processes. 
On  the  basis  of  the  preceding  considerations,  it  is  not  difficult  to  find  that,  in  the 
presence  of  white  noise,  the  matched  filter  is  defined  by  (except  for  the  above  factor) 

h ( v ) = S*(v) 
h ( t ) = S(-t) 


and  that,  as  a result,  it  performs  the  convolution  of  its  input  with  the  signal  (more 
precisely  with  the  signal  "reflected"  in  time  — the  image  of  the  signal). 
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Having  reviewed  these  results,  it  will  be  possible  to  deduce  them  as  well  as  all 
the  other  properties  of  matched  filtering,  from  the  generalization  of  the  theorem 
to  a system  of  N inputs  - a generalization  which  we  shall  now  consider. 


CHAPTER  HI 


MATCHED  FILTERING  WITH  N INPUTS  - FORMAL  SOLUTION 

Summary 

This  chapter  treats  the  problem  of  generalizing  matched  filtering  to  N inputs 
(multiple  filtering).  An  antenna  with  N sensitive  elements  exhibits,  at  all  of  these 
elements,  N parasitic  noise  components  whose  statistical  relationships  are  defined 
by  the  crosscorrelation  functions  of  the  noises  taken  two  at  a time.  It  also  delivers 
a "multiple  signal"  consisting  of  N particular  time  functions.  What  is  the  multiple 
filtering  of  these  data  - the  ensemble  of  N filters  - such  that,  when  the  N outputs  are 
summed,  the  characteristic  parameter  oi  matched  filtering(  S/N )is  optimized? 

The  formal  answer  to  this  question  is  presented  in  two  equivalent  forms  - the 
list  of  N impulse  responses  of  the  desired  filters,  and  their  N complex  gains 
(system  functions). 

In  the  second  form,  the  solution  is  seen  to  be  unique  (except  for  a real  factor 
and  a time  delay)  and  suitable  for  expression  in  convenient  matrix  notation. 

Two  characteristic  properties  are  established: 

a.  The  optimized  parameter  is  independent  of  the  reference  time  t 0 

b.  Except  for  a real  factor  and  a time  delay  (representing  an  arbitrary  choice 
of  the  time  origin) , the  spectrum  of  the  signal  and  the  spectral  density  of  the  noise 
are  identical  at  the  output  of  the  matched  filter. 

III-l.  Status  of  the  Problem 

Let  us  consider  a system  of  N inputs  Ej  ,E2  . . . . E„  where 

a.  the  signal  is  represented  by  N specific  real  functions  of  time, 

Sj(t),  S2(t),  ....  sN(t)  , 

b.  the  interference  (parasitic  noise)  is  represented  by  N unspecified  real 

functions,  Bj(t),  B2(t),  Br,(t)  , 

which  are  stationary  and  whose  correlations  are  stationary. 
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By  definition,  let" 


Cjr(t)  = E j B.(t)  Bk(t  + T)  } , (III-l) 

where  E designates  the  mathematical  expectation  (of  the  function  in  the  brackets). 


The  real  functions  C , , ( i ) are  autocorrelations  (even  functions)  when  j = k 

2 

and  crosscorrelations  when  j f k .The  N equations  of  the  type(  III-l ) express 
the  statistical  relationships  between  the  N noise  inputs.  These  noise  inputs  are 
independent  of  the  signals  S j ( t ) , S 2 ( t ) , . . 


We  then  ask: 


What  is  the  linear  filtering  process  represented  by  N filters  having  impulse 
responses  Ri(t),R2(t),....R^(t)  which,  when  applied  to  the  inputs  E i , E 2 . . . . I--,. , 
respectively,  have  the  following  property:  by  taking  the  sum  of  their  N outputs,  a 
common  output  ^ is  formed  in  which  the  following  ratio  is  optimized  (maximized): 

P=  instantaneous  power  of  the  signal  at  arbitrary  t0 


-.(111-2) 


average  noise  power 

Thus,  the  characteristic  criterion  of  matched  filtering  is  utilized. 

III-2.  Review  of  the  Classical  Properties  of  the  Crosscorrelation  of  Two  Noise  biputs 

Before  proceeding  to  the  proof,  let  us  recall  a few  classical  properties  of  the 

functions  C ( t ) and  C ( v ) . 

J k j k 


a.  First  of  all,  according  to  the  definition, 


cJkO>  - ckJ(-,>  . 

b.  Expressing  the  fact  that  the  power  of  a real  noise 


B (t)  + XBR ( t ) (x  real) 


(HI-3) 


2 


CXher  authors  adopt  a convention  which  corresponds  to  exchanging  i and  - i 
in  the  relation  TTT-1  . This  results,  for  the  remainder  of  this  paper,  in  the 
exchanging  of  the  correlation  matrix  with  its  transpose. 
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is  a positive  quantity  regardless  of  the  value  of  X , it  may  be  shown  (Schwartz 
inequality)  that 

[Vo]  <Cjj(0)  . Ckk(0),  for  all  t,  (III-4) 

the  right  side  of  the  inequality  being,  moreover,  the  product  of  the  powers  of  B ( t ) 
and  of  Bk  ( t ) . 

c.  The  existence  of  cjj(u)  , the  Fourier  transform  of  Cj  ( r ) , will  be  assumed. 
It  is  the  spectral  density  of  ( t ) , a real,  even,  non-negative  function  of  v 

for  any  value  of  v . 

d.  The  existence  of  C j k ( v ) > the  Fourier  transform  of  C j k ( t ) , will  be  assumed. 
As  the  transform  of  a real  function,  it  satisfies  the  equation 


C III— 5 } 


which  expresses  its  Hermitian  symmetry  with  respect  to  the  variable  v . Moreover, 
these  functions  have  Hermitian  symmetry  with  respect  to  their  indices  since  Eq.(  1 1 1-3  ) 
implies  that 


Cjk(v)  = Ckj(v)  • 


(iir-6) 


Since  the  Cj  ^ ( v ) are  obviously  spectral  densities  we  will  call  the  C . k ( v ) 
crosscorrelation  spectra. 

e.  If  B j ( t ) and  Bk  ( t ) are  filtered  by  filters  having  impulse  responses  Rj  ( t ) 
and  Rk  ( t ) } respectively,  two  new  noises  b ^ ( t ) and  t»k  ( t ) will  be  obtained  whose 
spectral  densities  are 

Y, .(v)  = C..(v)  | R , ( v ) |? 


>kk(v)  = CUk(v>  • 


( T 1 1 — 7 ) 


The  crosscorrelation  function  of  b,(t)  ;Uul  b,  (t)  is,  by  definition, 


12 


>Jk(T)  = E | bj(t)  . bk(t+x)  | . 


( II 1-8 ) 


A classical  calculation  (Appendix  I)  shows  that 


/ 


Rk(v) 


Vv) 


2tt1vt  , 
e dv 


( 1 1 1 — 9 ) 


from  which  we  derive  the  Fourier  transform  of  y ^ ( t)  ; 


* 

'j 


Yjk(v)  = R.,(v)  . R,,(v) 


(III-10) 


f.  Let  two  filters  R . ( v ) and  R^  ( v ) be  identical  with  real  impulse  responses 
(with  Hermitian  transfer  functions)  having  gains  equal  to  unity  in  a band  $ around 
v o (and  -vo  ) and  zero  elsewhere.  Let  us  assume  that  this  band  is  narrow  enough 
so  that  CJR(v)  does  not  differ  from  C . k ( v 0 ) within  the  band. 

The  crosscorrelation  function  obtained  for  the  two  noise  inputs  is,  according  to 

( III-9 ) and  (III-6)  , 

2*  . (CJk(v0)|  C0S{  2-it  v o t + Arg  [cjkCv0)]  \ , 

which,  as  4>^0  , approaches  the  limit 

y.  k ( t ) = 2ij)|Cjk(vo)|  COS  | 2 n v o t + Arg  (HI-11) 

(see  application  in  Chapter  IX); 

that  is,  it  approaches  a real,  sinusoidal  function  of  t whose  phase  is  equal  to  the 

argument  of  C ( v 0 ) and  whose  amplitude  is  pi'oportional  to  | C ( v 0 ) | . 

J k J k 

The  Fourier  transform  is  a spectral  "line"  at  the  frequency  v o characterized  by 
the  complex  number  C , . ( v 0 ) .Thus,  C , . ( v 0 ) is  a complex  number  which  describes 
the  crosscorrelation  of  noises  ( t ) and  Bj.  ( t ) by  two  identical  narrow -band 

filtering  processes  at  the  frequency  v o • 


The  function  y . k ( t ) given  by  ( 1 1 1 - 1 1 ) is  justified  by  the  inequality  ( 1 1 1 - 4 ) , 


which  becomes  here 

2 

[ YJk^r)  ] - Yjj(0)  • Ykk(0)  ’ 

from  which  may  be  deduced 


for  all  t 


lCjk(vo)l  1 Cjj(v0)  . Ckk (v0)  . 


( 1 1 1-12  ) 
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The  right-hand  side  is  the  product  of  spectral  densities,  for  v = v0  , of  the 
noisesB  ( t )and  B,  (t ).  Since  the  preceding  equation  is  true  for  all  o 0 ,we  may  write 


cjk(v)l  1 Cjj(v)  Ckk(v)’  for  a11  v’ 


(III-13) 


or  in  determinant  form 


Cjj 


1 0 » (III-13a) 

^k,J  ^kk  for  all  v . 

ni-3.  Review  and  Physical  Interpretation  of  a Fundamental  Property  of  the 
Crosscorrelation  of  Two  Noises 

Let  us  consider,  for  the  present,  the  problem  posed  in  Paragraph  III-l  and 
illustrated  by  Fig.  III-l  (multiple  filtering).  The  signal  at  the  output  of  filter  R,  (t ) 


J R^ta)  Sj(t-e)  de  , 

and  as  a result  the  signal  at  the  output  Y is 


o(t)  = Y f 


R (e)  Sj(t-e)  de 


The  noise  at  the  output  of  filter  R j ( t ) is 
y*Rj(u)  B . (t-u)  du  , 
and  the  noise  at  Y is 

B(t)  = Y [ R«  di 


(111-14) 


Its  average  power  is 


Pp  = E 6 (t)  = E 


R (u)  R ^ ( t-u ) di 
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The  [£]'  in  Eq  (III-16)  may  be  put  in  the  form  of  a double  £ : 

pe  = E E f /"Rj(u)  Rk(v)  * E|  * Bk(t-v)  j du  dv  , 

J k J J * ( III-17 ) 

obtained  by  permuting  the  symbols  } and  E,  which  is  legitimate  under 
many  realizable  conditions  in  physics.  We  have,  then, 


pe  -EE 

J k 


Rj(u)  Rk(v)  Cjk(u-v)  du  dv 


(III-18) 


from  the  definition  of  C j ^ ( t ) (Eq.  ( 1 1 1 -1 ) ) . 

The  right-hand  side  of  ( 1 11-18  ) being  a power,  is  non-negative,  regardless  of 

what  filtering  processes  R j ( t ) are  considered.  This  fact  constitutes  a characteristic 

property  of  the  ensemble  of  cross-  and  autocorrelation  functions  of  any  N real  noises 

B ( t ) . We  have,  then, 

J /•  /• 


E E 

j k 


//■*“ 


) • R (v)  C (u-v)  du  dv  > 0 , (III-19) 


for  any  R . . 

Equation  ( 1 11-19  ) will  be  later  used  to  advantage.  1 

HI -4.  Investigation  of  the  Solution  to  the  Problem  Posed  in  Paragraph  III-l. 

The  ratio  P which  is  to  be  maximized  by  the  proper  choice  of  the  filters  R j ( t ) 
is,  according  to  equation  (III-2)  , 


o(t0  ) 


(111-20) 


It  is  evident  that  this  ratio  P » except  for  a constant  factor,  depends  only  on  R,  (t)  , 
We  may  recall  the  remark  made  in  Paragraph  II-2  that  matched  filtering  would  be 
defined  except  for  a constant  factor.  In  order  to  maximize  p , this  arbitrary  factor 
must  be  chosen  in  such  a way  that  the  numerator  of  p , that  is,  ( 1 0 ) j , 
is  given  an  arbitrary  non-zero  value  k and  then,  taking  this  constraint  into  con- 
sideration, the  denominator  is  minimized. 

The  problem  is  now  reduced  to  finding  filters  R . ( t ) such  that 


°(to)  = E / Me>  s ( 1 0- e ) de  = k , 


( II 1-21 ) 
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and  such  that  the  quantity  Pg  in  Eq.  ( III-8)  is  minimized. 

Let  us  consider  the  ensembles  of  N impulse  responses  r*  1 C t ) , r2(t),  ..  r^  ( t ) 

such  that 

(III -22) 


r.(e)  Sj ( t o — 0 ) d6  = 0 ; 


that  is,  filtering  processes  which  null  (or  reduce  to  zero)  the  signal  at  the  output 
2 at  time  1 0 . 

It  is  clear  that  all  filtering  of  the  form 

R,  (t)  + a r.(t) 

with  some  given  real  a satisfies  Eq.  ( III-21  )}as  does  R ^ ( t ) itself.  By  varying  a 

and  r^(t)  (constrained  only  by  Eq.  ( II 1-22  ))  we  may  associate  with  any  filtering 

process  R . ( t ) an  infinite  number  of  filtering  processes 
J 

Rj(t)  + ai'j(t)  , 

2 

each  one  of  which  assigns  the  same  value  k to  the  numerator  of  P and  different 
values  to  the  denominator  P g . 

If  the  chosen  filtering  R ( t ) is  indeed  a solution  to  the  problem,  that  is,  if  it 

minimizes  P it  means  that  the  value  of  PQ  for  that  R (t)  is  less  than  or  equal  to 
p p J 

all  other  values  obtained  for  filtering  processes  R ^ ( t ) + ar  . ( t ) 

with  arbitrary  real  a andr^  ( t ) constrained  by(  II 1-22  ),  Now,  calling 

h j ( t ) > one  of  the  desired  filtering  processes, 

, the  value  of  P g for  R ^ ( t ) = h ( t ) , and 

D the  value  of  P Q forR.(t)  = h.(t)  + ar.  (t). 
n » r p J j J 

the  desired  condition  may  be  written 


D,  < ' Du 
h — h ,r 


(HI-23) 

for  arbitrary  a and  for  an  arbitrary  value  of  r constrained  only  by  ( 1 1 1-22  ) . 
Substituting  from  Eq.  ( III-18  ^the  inequality  ( III-23  )becomes 


E // 

Jk 


h ^ ( u ) hk(v)  Cjk(u-v)  du  dv 


(III-24) 


1 E//[h- 

jk*7*7  L 


j ( u ) + ar^( 


u)  hk(v)  + ark(v)  C-k(u-v)  du  dv 
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for  any  a and  any  r constrained  by  ( II 1-22 ),  or  rearranging  the  above  in 
terms  of  a , 


E // 
Jk  J J 

£//[ 


r.(u)  . r,  (v)  C.,  (u-v)  du  dv 
J k 


+ a V II  I r j ( u ) hk(v)  + hj(u)  rk(v)  I Cik^u-V^  du  dv  >,Q 
j k ' ' 


(HI-25) 


for  any  a and  any  r constrained  by  (HI-22). 

2 

The  coefficient  of  a in  ( III-25 ) is  of  the  form  found  in  Eq.  (III-19)  ; 
hence  it  is  non-negative.  It  represents,  moreover,  the  noise  power  at  the  output  £ 
for  the  filtering  process  r . ( t ) . Since  the  first  term  of  Eq.  ( III-25 ) remains 

u 

positive  for  all  a , in  order  to  prevent  the  sign  of  the  second  term  from  prevailing 
for  small  values  of  | a | , it  is  necessary  to  set  the  coefficient  of  a equal  to 
zero.  This  condition  is  then 


r,(u) 

J 


hR(v)  + h . ( u ) rR (v) 


C . k ( u-v ) du  dv  = 0 , 

(III-26) 


for  any  r constrained  by  ( II I -2 2 ) . 

HI-5.  Demonstration  - Temporal  Form  of  the  Solution 

The  first  term  (left  side)  of  ( 1 1 1-2  6 ) is  composed  of  two  terms.  Because  of 
Eq.  ( 1 1 1 - 3 ) , which  implies 


cJk(u-v)  = Ckj.(v-u)  , 

these  two  terms  may  be  written  in  symmetric  form  with  their  two  indices,  as  well  as 
u and  v permuted. 

It  is  evident,  in  this  form,  that  they  are  equal.  Condition  ( 1 1 1-26  ) is  then 
reduced  to 


(u) 


hk(v)  C^k(u-v)  du  dv  = 0 , 


(III-27) 


for  any  r constrained  by  (II I -2 2) 
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r 


Equation  (III— 27)  may  be  written 


£/r  j(u)  [e  / \ 


(v)  C (u-v)  dv  du  = 0 


(III-28) 


and  from  ( II 1-22 ) 


.1 


( u ) S . ( t 0-u ) du  = 0 


( III-22 ) 


In  order  to  simplify  the  notation  for  the  moment,  let  us  set  the  bracketed  term 

of  (III-28)  equal  to  f,(u).  Among  all  the  possible  systems  of  functions  r ,( u ) 

J J 

let  us  consider  those  which  possess  only  a single  function  not  identically  equal  to 

zero,  for  example  rj(u).  This  function  is  constrained  by  Eq.  (III-22),  which  is 

expressed  as 

/ r i (u)  Si  (t0-u)  du  = 0 . (III-22a) 


Condition  ( II 1-2 8 ) defines  only  f j ( u ) ; for  any  function  r x ( u ) constrained  by 
(III-22a),  fi(u)  must  satisfy 


J ri(u)  fi(u)  du  = 0 . (III-28a) 


Equations  (III-22a)  and  ( II 1-2  8a)  are  .seen  to  result  in 


f 1 (u)  - Ki  Sj(tQ-u), 


(111-29) 


where  K1  is  a factor  independent  of  u (Appendix  EE).  The  preceding  relation  is  valid 
for  all  values  of  the  index  j , and  as  a result 

fj  (u)  = K.  (t0-u) 

Applying  these  necessary  conditions  to  ( 111-28  ) , we  obtain 


£ 


( u ) S, ( t o— u)  du  = 0 . 


Comparing  with  ( 1 11-22  ) we  observe  that  all  K . must  be  equal  to  one  another  and  to 
the  same  factor  K , which  is  independent  of  the  index  j . 

Finally,  the  system  of  equations  which  gives  the  impulse  responses  hR(t ) , 

the  solution  to  the  problem,  is  the  ensemble  of  equations  of  the  type 


1!) 


( 1 11  — 30 ) 


32  hR(v)  C j k ( u-v)  dv  = K S , ( 1 0-u)  , 

k 

for  J = 1,  2,  N . 


These  N integral  equations  define  the  N impulse  responses  h j ( t ) ,h2  ( t ) h[;  ( t ) 

and  give  the  desired  solution  to  the  problem  in  its  temporal  form  (impulse  responses). 

Optimal  filtering  comprised  of  the  filters  ( t ) is  defined  by  Eq.  (HI-30 ) , 
except  for  a factor  K . This  factor  is  arbitrary  but  it  is  still  related  to  the  other 
arbitrary  factor  k = a ( t 0 ) by  Eq.  ( 1 1 1 -2 1 ) by  which  the  numerator  of  p was 
normalized;  it  is  related  as  well  to  the  noise  power  at  the  output  ^ (denominator 
of  o ) , which  is  minimized  for  a given  K but  which  is  nevertheless  proportional  to  K * . 

In  fact,  Eq.  (II I -21)  is  written,  for  the  filtering  process  h^(t), 
k = o(t0)  = 

which  is,  multiplying  both  sides  by  K, 


E/ 


h,(u) 


Sj ( t 0-u ) du 


(III-31) 


Kk 


32  f h ^ ( u ) £ K Sj(t0-u)J 
j 


du 


(III-32) 


and,  replacing  the  bracketed  term  of  (III-32  ) by  its  value  derived  trom  (III-;), 


kk  - E £ 

j k 


//h’<u 


) h.  (v)  C..  (u-v)  du  dv  . 

k (111-33) 


Thus,  the  second  term  (right  side)  of  (111-33)  is  the  value  of  power  P for 

P 

the  optimum  filtering  process  ( t ) , that  is,  the  noise  power  at  output  32 
-power  which  is  minimized  by  that  optimum  filtering.  The  minimum  power  is 
designated  by  . 

We  have  then 


or 


K° ( t o ) = Pm  . 


As  for  the  value  of  the  ratio  p , which  is  made  a maximum  p^  by  the  optimal 
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filtering  h ( t ) , it  is  equal  by  definition  to 


o ( t n ) 


PM 


o(tn) 


These  equations  may  be  summarized  by 

a{  t0) 


o(tn  ) 


(111-3*0 


and  the  proportionality  of  P to  K may  be  made  obvious  by 


P = pM  K . 

m M 


( TIT-35 ) 


The  coefficient  K has  the  dimensions  of  signal  amplitude. 

m-6.  Spectral  Form  of  the  Solution  - Value  of  the  Optimized  Parameter 

By  replacing  both  sides  of  Eq.  ( 1 1 1- 30  ) by  their  Fourier  transforms,  we 
obtain  the  solution  to  the  optimum  filtering  problem  in  "spectral”  form.  The  left 
side  of  (III-30)  contains  the  convolution  of  h^(t)  with  ^ ( t ).  Its  transform 
is  the  product  of  the  Fourier  transforms  ( v ) and  C j ^ ( v ) . The  transform 
of  the  first  term  is  then 


£ Mv)  c (v) 


If  S ( V ) is  the  Fourier  transform  of  Sj  ( t ) , , the  Fourier  transform  of  the  right 
side  of  (HI-30)  is 


S* ( v ) e-2'lvt°  , 

where  the  asterisk  denotes  the  complex  conjugate.  Accordingly,  the  system  of 
Eqs.  (II 1-30)  yields,  through  transformation,  the  equations 


V C,,  (v)  h„(v)  = K e“2irivt0  s*(v),  for  j = l,2...N 

~ Jk  k J C 111-36 ) 

k 

In  Eq.  ( 111-36  ) the  quantities  ( v ) (system  functions  of  the  filters)  are 


dimensionless  since  they  are  ratios  of  two  amplitude  spectra.  The  C , , ( v ) 


r 


are  homogeneous  in  spectral  density,  that  is,  in  power  per  cycle  (power  x time); 
hence,  in  energy. 

The  S ( y ) are  in  terms  of  amplitude  per  cycle  (Fourier  transf onus  of 
amplitudes),  and  K is  the  unit  of  amplitude  of  the  signal. 

One  may  write,  introducing  the  unit  of  time  T , 

£ CJk(v)  hk(v)  = K2T  e'2nlvt° 

2 k 

where  K i is  the  unit  of  energy  or  spectral  density. 

In  terms  of  dimensionless  quantities,  Eq.  ( - -3h  ) may  be  written 


SJ(v) 

KT 


£ 

k 


2 

K T 


hk(v) 


-2  n i vt  o 
e 


KT 


(III-36a) 


However,  this  form  is  less  concrete,  and  we  will  retain  the  prerogative  to  speak 
of  system  functions,  densities,  and  signals. 


The  system  (I 11-36)  is  a system  of  N linear  equations  in  N unknowns  consisting 
of  the  h ^ ( v ) . . Solving  explicitly,  we  have  the  N equations 

M»)  = Z S*N)  ->JkCv)  , (HI-37) 

J for  k = 1,2,  N, 

where  a ( v ) denotes  the  determinant  of  the  C . ( v ) ( j is  the  line  index  and  k 
the  column  index) , and  where  M .J’J)  denotes  the  cofactor  of  the  element!  j , k ) 

K 1-  j-j  th 

of  A(  v ) (that  is,  the  determinant  obtained  by  suppressing  the  j line  and  the  k 

i +k  \ 

column  and  assigning  the  sign  (-1)  J ) . 

The  system  of  equations  given  by  (II  I- 37)  constitutes  the  "spectral"  form  of 
the  solution.  It  defines  the  optimum  filtering  process  by  the  system  functions  of  the  N 
filters,  , hjfv),  h2(v),  hN(v). 

Let  us  note  that  the  solution  of  the  linear  system  ( 111-36)  is  generally  unique 

provided  that  A ( v ) ^ 0 ) . 

Thus,  the  ( ^ ) are,  in  general,  completely  determined,  as  are,  consequently, 
the  h k ( t ) . One  may  then  speak  of  the  solution  to  the  problem  of  matched  filtering 
with  N inputs  - a unique  solution  defined  by  ( 111-36  ) . 

The  case  of  A ( v ) = 0 will  be  investigated  later.  The  C ^ k ( v ) are  Ilermitian 

with  respect  to  the  variable  v (Eq.  (II 1-5)  ) . The  Sj  ( v)  are  also,  as  transforms 
of  the  real  signals  Sj  ( t ) . Thus,  all  the  system  functions  h^,  ( v ) are  Ilermitian  with 
respect  to  v,  and  the  impulse  responses  h (t)  are  real.  Driven  by  real  signals  or 
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noise,  an  (t)  delivers  a real  signal  or  a real  noise. 

An  interesting  expression  for  the  maximum  value  of  PM  of  the  ratio  p , 

, may  be  drawn  from  the  Eqs.  ( 1 1 1 - 3 7 ) . 

In  fact,  Eq.  (1 11-31)  may  be  written 


optimized  by  the  filtering  process  j h. 


(v) 


o(t0)  = E f Sk(u)  hR(t0-u)  du  , (III-31) 

k J 

and  when  transformed  by  Parseval's  identity  gives  the  result 


o(t0)  = E f Vv)  e2"ivt0  dv 

k 


(111-38) 


(the  Fourier  transform  of  h 

>.  -2l.lvt0 
h,  ( v ) e u 

k 


Substituting  the  value  of 


R(t0-t)  being  hR(-v)  e?  * 1 vt  0 , , that  is, 

, since  h^(t)  is  real,  as  we  nave  seen  previously). 

h^Cv)  given  by  (111-37)  into  Eq.  (111-38),  we  have 


j 


°Uo)  » £ f Sk(v)  ^ [ £ 

k 

or  from  ( 1 1 1 -3^  ) 

0 ( t 0 ) 


S.(v)  MJk(v)  dv 


M 


k 


. Vv>  M,k(v) 


a v 


J k 


A ( v ) 


‘ ( I T I — 3 9 ) 


Note  that  since  pm  is  a ratio  of  powers  the  right  side  of  ( 1 1 1-39  ) is  always 
non-negative  for  any  noise  B ^ ( t ) and  any  signal  •'  5 ( t ) • 

The  value  0 r,i  is  the  integral  of  a function  of  v : 


D(v) 


* M1k(v) 

EE  sj(v)  Vv)  — — 

j k A ( v ) 


(III-HO) 


whose  real,  non-negative  character  for  all  v will  become  evident  shortly 
(Eq.  111-52)). 

The  output  signal  0 ( t ) has  an  amplitude  spectrum  given  by 

o(v)  = Yj  Sk(v)  \(v)  (111-41) 

k 
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or,  taking  (I IT-37)  and  II 1-4 0)  into  account, 

o(v)  = ^Ke_27rivt°^  . D(v)  . (111-42) 

III— 7-  Matrix  Notation  - Characteristic  Property  of  Matched  Filtering 

We  will  now  express  the  preceding  equations  in  a particularly  convenient  matrix 
notation  that  is  suggested  quite  naturally  by  Eq.  (HI - 36). 

In  this  notation,  the  table  of  Cj^(v)  constitutes  a square  matrix  of  order  N,  t 
which,  according  to  Eq.  (III-6)  „ is  Hermitian.  The  ensemble  of  system  functions 
Rj(v)  or  h^(v)  may  be  represented  by  column  matrix  R or  h.  The  ensemble 

of  signals  may  be  designated  by  a column  matrix  S . 

It  is  obvious  that  the  matrices  S and  R represent  different  physical  quantities 
(signals  and  system  functions).  Thus,  all  operations  between  matrices  S and  R 
are  not  justified  from  a physical  point  of  view.  For  example,  a linear  combination 
aS+BR  has  no  meaning  and  will  not  be  used.  On  the  other  hand,  a matrix 
multiplication  such  as  Sn  or  hS  represents  the  right-hand  side  of  (III- 41)  , 
that  is,  the  amplitude  spec  trum  of  the  output  signal  of  the  matched  filter.  More 
generally, 

SR  = RS 

designates  the  signal  output  of  the  multiple  filtering  process  of  Fig.  (III-l).  Moreover, 
using  only  dimensionless  quantities  (Eq.  ( III-36a)  ) , only  physical  quantities  of  the 
same  nature  may  be  manipulated. 

Anticipating  its  justification  given  in  Paragraph  V-4,  let  us  use,  for  the  moment, 
Eq.  ( V-21 ) , which  gives  the  spectral  density  of  the  multiple  filtering  process  R 
at  the  output  ^ , 

For  any  arbitrary  R^(t)^R  (v)  , we  have 

yE(v)  = ££  Rj(v)  Rk(v)  C)k(v)  * (V-21) 

J k 

which  is  expressed  in  matrix  notation  as 
3 

The  symbols  of  matrix  calculations  are  those  recommended  by  Standard 


NFX02-110  of  April  19G2. 
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(111-^3) 


Y £ ( v ) = R+CR  . 

The  right-hand  side  is  the  Hermitian  form  of  the  matrix  C . In  the  case  of  matched 
filtering,  the  spectral  density  of  the  output  noise  is 


Y_  (v)  = h Ch 
1 0 


(III-44) 


E quation  (111-36)  is  written 


= ^Ke-2nivt°]  S*  . 


Ch  = | Ke-"  " -1- ' “ u | s . (111-45) 

(In  matrix  expressions,  coefficients  will  be  placed  in  brackets  for  greater  clarity.). 
Equation  (III-37)  gives,  when  A(v)  / 0 

h = jKe-  7Tivto  J C"1  S*  . ( III-46 ) 

Matrix  C being  Hermitian,  we  have 

C+  = C and  C"1+  = C"1.  (111-47) 

Setting  the  matrices  associated  with  both  sides  of  ( ITT  — 46)  equal  to  one 
another,  we  have 


i+  = ^Ke+2nlvt0] 


#+  - 1 + 

S C 


(ITT- 4 8) 


.*  + 


Since  S = S , the  transpose  of  S , we  have 


i+  = [^Ke+2lllvt0]  S C"1  . 


Equation  (III-44)  is  then  written 

YEq(v)  = h+  Ch  = [k?]  S C"1 

M,,  (v) 


(111-49) 


(TIT-50) 


Substituting  in  ( 1 1 1 - 4 0 ) and  noting  that 


A(v) 


is 
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the  term  corresponding  to  "line  k,  column  j"  of 


l .we  may  write 


D( v ) = S C'1  S*  , 

and,  consequently,  from  (II 1-50) 


(HI-51) 


Y_  (v)  = K2  D(v) 
1 0 


(III-52) 


As  a result  of  (HI-52),  D ( v ) is  indeed  a real  and  non -negative  function  of  v 
and  is  homogeneous  at  a given  time. 


Comparing  (IT  1—42)  and  ( 1 11-52  ) , one  is  led  to  the  characteristic  rule  of 
matched  filtering,  generalized  here  to  the  case  of  N inputs:  except  for  the 


factor 


[ 


-2  iti  vt  n 

e u 


J * , the  spe 


spectrum  of  the  output  signal  is  identical  to  the 


spectral  density  of  the  noise  at  the  output. 


Y 


Z 


0 


(v) 


0 ( V ) 

-2 it i vt  o 


D(v) 


( TTI-53) 


The  identity  becomes  complete  when  one  adopts  the  conventions 

K = 1 and  t o = 0 . 


The  first  convention  consists  of  defining  the  unit  of  amplitude;  it  destroys  the 
homogeneity  of  the  preceding  equations.  The  second  expresses  the  fact  that  the  time 
origin  is  taken  to  be  the  instant  when  the  "peak"  of  the  output  signal  appears.  This 
convention  of  notation,  which  does  not  alter  the  generality  of  the  problem,  will 
ultimately  be  adopted. 

A convenient  integration  of  ( 1 1 1 — d 2 ) yields  the  result 


o(t9) 


2u i vt  o 
e 


dv 


/ 


KD(v)  dv 


K PM  ’ 


and  integration  of  ( n I -52  ) gives 


P = K pM  . 
m M 


Thus,  we  again  encounter  Kqs.  (Ill  — 3;0  and  ( 1 I 1-35)  . 


20 


(iri-54) 


Note  that  the  matrix  form  of  Eq.  (111-41)  is 

o(v)  = Sh 

and  that  the  comparison  of  (III-44),  (113-53),  and  ( ! ! I - 54  ) gives 

yEq(v)  = h+  Ch  = ^Ke‘  11 1 vt  0 j sh  . ( II 1-55 ) 

This  equation  may  be  demonstrated  directly  by  beginning  with  Eq.  ( 111-45 ) , 
setting  the  two  matrices  associated  with  the  two  sides  equal  to  each  other,  and 
multiplying  on  the  right  by  h . 

m-8.  Properties  of  the  Output  Signal 

Equation  (111-34)  shows  that  o(t0)  has  necessarily  the  same  sign  as  K,  , 
which  we  will  assume  to  be  positive  in  order  to  conform  to  the  current  practice  which 
displays  the  "peak"  of  the  signal  "upward";  however,  the  opposite  convention  would  be 
just  as  useful.  Let  us  justify  the  qualitative  term  "peak"  showing  that  o (t  0 ) 
is  indeed  the  maximum  value  of  the  output  signal  o ( t ) . This  is  equivalent  to 
showing  that  o ' ( t ) = ( t 0+t ) is  less  than  or  equal  in  magnitude  to 

o ' ( 0 ) = o(t0) , 

or  |0'(t) | i o '(0)  . (III-56). 


Then, 

„’(t)  • - /.(,)  + dv  , 

and  using  equation  ( 111—42)  , 

o * ( t ) = K /*D(v)  e2lTlvt  dv  , 

which  is,  as  we  know,  real  and  non-negative,  and  the  obvious  inequality 

j J D(v)  e‘  7,^vt  dv|  < J"  D(v)  dv 

precisely  demonstrates  Eq.  (111-56)  • Note  that  KI.i(v),  being  even  in  v, 
has  a Fourier  transform  o’  ( t ) that  is  even  in  t . Thus,  0 ( t ) symmetric 
about  time  t = t o . 
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APPENDIX  I 


We  have,  in  fact 


3 


b«(t)  = J B,(e)  R,(t-e)  de  , 


r ' j 


bR(t)  = J Bk(e')  Rk(t-e')  de'. 


from  which  we  get 


VT>  ■ ff 


E j B.(e)  B,  ( e ’ ) | R.(t-e)  r (t+x-e')  de  de’  , 


YJk(t)  = JJ  c . R ( e ' - e ) Rj(t-e)  Rk(t+x-e’)  de  de  ’ . 
Letting 


we  have 


e ' - e = o , 

t - e = to  , 


■ ff 


C.,  (o)  R . ( (D ) R,  (t  + u-o)  do  dw 
0 k J K 


Then, 


Rk(t)  = J Rk(v)  e2TTivT  dv  , 

from  which  we  get 


Y j k ( t ) = J j ' ^*Cjk(o)  Rj(u)  Rk(v)  e'Triv^T  lb  du  do  dv 


a u v 


Y * ,,  ( T ) 


■/ 


if- 


Rk(v)  e27,ivT|  fRj(o)  e2"1™  du 


][/Ctk(a)  e'2lTiV0  doj 


d v 


28 


The  first  bracketed  term  is  R j ( v ) • 
The  second  bracketed  term  is  C j ^ ( v ) • 


Then, 


Tjk(T) 


-fR' 


(v)  Rk(v)  CJk(v) 


2iTivT 

e 


This  example  is  analogous  to  that  of  reference  (7),  page  39. 


APPENDIX  II 


S ( u ) being  a given  function,  let  us  determine  a function  f ( u ) such  that, 
for  any  function  r(u)  constrained  by  the  relational* ( u ) S(u)  du  = 0 , 


we  have  / r ( u ) f(u)  du  = 0 . 


There  are  two  points  u =u,and  u = U2  for  which  the  values  of  S(u)  are  S(uj ) 
and  S(u2  ) (which  will  be  assumed  to  be  different  from  zero).  Consider  the  function 

6 ( u-u  j ) 6 (u-u2 ) 


r(u)  = 


S ( u i ) S(u2) 


This  satisfies  the  first  equation  since 


J’ 5 ( u-u 


j)  S(u)  du  = S(uj) 


We  must  have  therefore 


— — f 6(u-u 

Ml)  J 


j)  f(u)  du  - 


S(u2) 


6(u-u2)  f(u)  du  = 0 


f(Uj) 

S ( u i ) 


f (u2 ) 
S ( U 2 ) 
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Considering  u2  as  an  arbitrary  point,  it  is  clear  that  s|yy  is  a constant, 

independent  of  the  point  considered  and  thus  independent  of  u . 

Therefore,  f(u)  = K S(u)  . 


CHAPTER  IV 


REMARKS  CONCERNING  THE  FORMAL  SOLUTIONS  OF  SOME  SPECIAL  CASES 
Summary 

The  preceding  formal  solution  may,  of  course,  be  applied  to  the  classical  problem 
of  a single  input. 

a.  In  the  special  case  where  the  N noise  inputs  are  uncorrelated,  it  is  shown 
that  multiple  matched  filtering  is  achieved  by  summing  the  outputs  of  separate 
matched  filters  in  each  channel. 

b.  H the  N noise  inputs  are  both  uncorrelated  and  of  the  same  uniform  spectral 
density,  the  generalization  of  a weU-known,  essential  property  of  the  classical 
theorem  establishes  the  fact  that  the  value  of  the  optimized  parameter  depends  only 
upon  the  total  energy  of  the  multiple  signal  (of  the  N signal  inputs). 

c.  The  parameter  pm  expresses  the  performance  (S/N)  of  the  matched  filter  for 
given  signal  and  noise  inputs.  If  the  noise  is  defined  and  the  choice  of  a signal  is 
arbitrary  (a  case  frequently  found  in  practice),  the  value  of  Pm  depends  upon 

that  choice. 

Examination  of  the  expression  for  Pm  shows  that  it  is  advantageous  to  choose 
signals  in  selected  frequency  bands  whose  specification  is  based  upon  the  correlation 
matrix  of  the  N noise  inputs.  This  qualitative  observation  introduces  the  idea  of 
"eliminable  interference"  or  "infinite  signal -to-noise  ratio,"  which  wiU  be  developed 
later. 

IV-1.  Limitations  of  the  Cases  Considered  in  This  Chapter 

The  systems  of  Eqs.  (HI-30)  and  ( 1 T 1 — 37)  provide,  in  temporal  form  and 
spectral  form,  respectively,  the  solution  to  the  problem  posed  in  Paragraph  HI-1. 

First,  let  us  note  that  according  to  ( III  — 37 ) the  system  functions  of  the 
filters  hR  ( v ) are  defined  only  for  i(v)/0  . Therefore,  it  is  a good  idea  to 
examine  the  case  in  which  this  condition  may  not  be  fulfilled.  This,  however,  is 
part  of  the  study  of  the  pi-operties  of  a ( v ) that  will  be  undertaken  later. 
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IV -2.  The  Case  of  a Single  Input 


With  the  aid  of  Eqs.  (HI-30)  and  (III-37  ) , the  solution  of  a matched  filter 
with  a single  input  may  be  found  - a solution  summarized  in  Paragraph  II-2. 

Equation  (III-3)  yields  Eq.  ( 1 1-2 ) directly  for  the  case  of  N = 1 , by  letting 

C 1 1 ( t ) = C(t)  , 

the  autocorrelation  function  of  the  single  noise  input. 

Equation  (III- 37 ) no  longer  has  any  meaning  in  the  case  where  N = I since 
the  cofactor  j ( v ) is  not  defined,  but  Eq.  (111-36)  is  reduced  to 

-2-rrivto  * 

h.  (v)  Cn(v)  = Ke  S j ( v)  , 


-2  -rrivt 

h j ( v ) C 1 x ( v ) 

= Ke 

-2tt1  vt 

hj  (v) 

= Ke 

sj(v) 


(IV-1) 


which  is  identical  to  ( II-3 ) since  C j j ( v ) = C ( v ) is  the  spectral  density  of  the 
single  noise  input  considered  here. 

The  rule  given  in  Paragraph  III— 7 is  also  valid  since  the  spectrum  of  the  output 
signal  is 

-2TTivt0  I s ( V ) I 2 

o(v)=h(v)S(v)=Ke  • y ' — , 

while  the  spectral  density  of  the  output  noise  is 

2 

>£„  ‘ c(v)  • lh(v)l  ’ K ‘ 1 cCv]  ' 

Wlth  '>('■)  ■ 1 C(  v j1’  • 

IV -3.  Matched  Filter  Rule  for  Uncorrelated  Noise.  Signal  Energy  and  Signal -to-Noise 
Ratio 

Let  us  consider  the  simple  case  where  the  N noise  components  are  not  correlated. 
In  this  case,  since  the  crosscorrelation  functions  are  identically  zero,  all  the 
elements  C , (v)  of  the  determinant  A(v)  are  zero  for  j / k . Thus,  A(v) 

J K 

becomes  diagonal  ;uid  takes  on  the  value  of  the  product  of  the  J spectral  densities 


The  cofactors  M j ^ ( v ) are  themselves  zero  for  j / k , and  one  cofactor  NL  j ( v ) 
is  the  product  of  N - 1 spectral  densities  obtained  by  omitting  Cjj(v)  . As  a result, 
Eqs.  ( III-37 ) are  reduced  to 
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hk(v)  = Ke 


Sk<»)  Mkk(„) 


that  is,  to  the  N equations 
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hR(v)  = Ke 


S*(v) 

* 


( IV-2 ) 


for  k=l,  2,  . . . . N . 

The  system  of  Eqs.(  IV-2  )points  out  that  for  each  input  the  optimum  filter  is  the 
same  as  the  case  in  which  that  input  is  considered  to  be  isolated.  The  following  rule 
may  then  be  stated  (referring  to  Fig.  III-l): 

A matched  filtering  process  with  I inputs,  whose  noise  inputs 


uncorrelated,  may  be  obtained  by  taking  the  sum  of  the  outputs  of 


the  matched  filters  for  each  input. 


In  the  special  case  considered  here,  the  value  of  Pm  given  by  expression ( ' ! ‘ 3 ' ^ 


is  reduced  to 


f / 


| S j ( v ) | 


( IV-3 ) 


Furthermore,  if.  the  N spectral  densities  are  uniform  and  equal  s 


CjjU)  d ’ 


and  we  have 


Pm  = d 
m 


f/ 


| 3 . ( v ) J d v . 


( IV— h ) 


The  integral  J \ S , ( v)  | dv  represents  the  energy  of  the  signal  S(t) 
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(Parseval's  identity).  Then,  the  ratio  p ^ depends  only  upon  the  sum  of  the  energies 
of  the  N signals,  that  is,  upon  the  total  signal  energy  at  the  inputs. 

This  statement  generalizes  the  equivalent  rule  valid  for  the  case  of  one  input  as 
long  as  the  noise  has  a uniform  density.  We  have  then 

E 


s 


where  Eg  is  the  energy  of  the  signal  (see  reference  (8),  Chapter  IV) . 

IV -4.  General  Case  - Usefulness  of  the  Arbitrary  Parameter 

A matched  filter  with  N inputs  is  defined  except  for  a constant  factor  - designated 
here  by  K - for  all  filters  hj.  t and  an  arbitrary  time  lag  represented  by  the 

-2  TTivt  n 

term  e u , a factor  in  the  N system  functions  h^(v)  . This  time  lag 

depends  upon  the  instant  of  time  t 0 chosen  for  the  appearance  of  the  "peak"  output 
signal.  As  has  been  already  suggested  in  Et-2,  the  performance  of  the  system  rep- 
resented by  Pm  - which  may  be  defined  as  the  signal -to-noise  ratio  of  the  process  - 
is  independent  of  K and  of  t 0 , which  is  evident  from  Eq.  ( III-39  1 . It  is  equally 
clear  physically  that  no  change  in  the  detection  ability  of  the  system  is  brought  about 
by  inserting  the  same  ideal  delay  line  in  series  with  each  input.  The  appearance  of 
the  signal  peak  is  only  changed  in  time  without  modifying  the  height  of  that  peak  above 
the  mean  noise  level.  By  this  observation,  we  touch  upon  the  question  of  the  possi- 
bility of  achieving  filters  defined  by  Eqs.  ( III-30)  or  ( III-37  ) . We  know  that, 
in  order  to  be  realizable,  a filter  must  have  an  impulse  response  which  is  zero  for 
t £ 0,  bounded,  and  unconditionally  integrable.  If  a theoretical  filter  satisfies  only 
the  two  latter  conditions,  its  unit  impulse  response  goes  to  zero  at  the  limits 
t = +°°  and  U = . It  is  possible  then  to  satisfy  the  first  condition  in  an  approxi- 

mate way  by  a translation  of  the  unit  impulse  response  to  the  right  on  the  time  axis, 
that  is,  by  a supplementary  time  delay.  This  amounts  to  saying  that  the  physical  ap- 
proximation of  a filter  is  often  facilitated  by  adding  adequate  time  delay.  The  approxi- 
mate realization  of  the  matched  filter  will  thus  be  facilitated  by  permitting  an  appro- 
priate value  of  t0  ; this  time  delay,  although  important,  certainly  constitutes  an  incon- 
venience since  it  postpones  the  instant  when  the  observer  is  informed  of  the  presence 
of  the  signal,  but  it  allows  us  to  obtain  a good  approximation  of  the  ideal  performance 
expressed  by  o . 

IV-5.  Limited  Possibilities  for  Elimination  of  Noise  Having  a Particular  Structure 
Let  us  consider  Eq.  ( III-39). 
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(Ill— 39) 


pm  - fz  Z SJ*U)  V»>  Vv) 

J j k 

A ( v ) 

The  possibility  appears  here  of  having  pm  become  "infinite"  if  there  exist 
frequencies  for  which  A ( v ) becomes  zero  (provided  that  the  signals  indeed  possess 
spectral  components  at  these  frequencies  and  that  the  sum  of  the  cofactors  does  not 
go  to  zero  as  well). 

Such  a possibility,  which  corresponds  to  the  total  elimination  of  noise,  naturally 
constitutes  a limiting  case.  It  is  apparent,  however,  that  it  may  be  interesting  in 
practice  to  concentrate  the  signal  energy  in  those  frequency  bands  where  A ( v ) 
is  very  small,  if  such  bands  exist. 

Formally,  structures  having  "eliminable  noise"  may  bo  defined  in  terms  of  matched 
filtering.  The  case  of  "coherent  interference"  which  belongs  to  this  category  will  be 
treated  later  (Paragraph  VII-5).  It  can  be  immediately  seen  that  no  such  situation 
ever  arises  in  the  case  of  N uncorrelated  noises;  in  fact,  as  we  have  seen,  A ( v ) 
is  reduced,  in  such  a case,  to  the  product  of  spectral  densities;  that  is,  to  a non- 
negative quantity.  Further,  the  spectral  density  of  noise  at  an  input  is  never  zero  in 
practice  (because  of  thermal  noise).  Thus,  A ( v ) may  not  be  zero  at  any  frequency 
(Paragraph  VI-8). 

Formally,  the  case  a(v  ) = 0 may  be  encountered  if  the  noise  inputs  are 
statistically  related.  However,  in  practice,  there  will  always  be  one  part  of  the  total 
interference  - thermal  noise  - which,  because  of  its  independence  from  one  input  to 
another,  will  prevent  A(v)  from  being  reduced  to  zero.  Thus,  A(v  ) = 0 appears 
as  a limiting  case,  interesting  in  its  formalism  but  physically  inaccessible,  as 
is  the  "infinite"  signal-to-noise  ratio. 

This  restriction  does  not  take  away  any  of  the  practical  interest  from  the  following 
statement:  the  signal-to-noise  ratio  tends  to  become  very  large  when  A(  v ) 

becomes  very  small.  In  order  to  use  this  fact,  we  must,  given  the  choice,  use  signals 
in  the  band  or  bands  where  A ( v ) is  the  smallest;  as  a result,  we  must  know  A ( v ) . 

This  preceding  statement  is  obviously  a generalization  of  the  very  simple  ease  of 
one  input  for  which  the  expression  in  ( 111  — 3^)  becomes 


35 


Applying  the  theorem  of  the  mean,  it  may  be  shown  that  for  a given  signal  energy, 


Es=  J | S( v ) | dv  . 

The  maximum  value  for  pm  will  be  obtained  by  placing  this  energy  in  that  band  - or 
in  the  limiting  case,  at  that  frequency  - where  C ( v ) is  minimized.  Since  it  is  a 
question  of  a limiting  case,  it  is  particularly  interesting  to  know  the  nature  of  the 
statistical  relationships  - let  us  say,  the  spatial  structure  of  the  noise  at  the 
inputs  - which  will  result  in  A ( v ) being  zero.  In  fact,  the  performance  of  matched 
filtering  is  then  considerable  since  it  permits  a theoretically  infinite  signal -to -noise 
ratio.  In  other  words,  it  allows  us  - in  this  case  only  - to  suppress  the  noise.  A 
more  complete  discussion  of  this  case  will  appear  later. 

At  present,  let  us  be  satisfied  to  take  note  of  the  formal  difficulty  involved  in 
practical  use  of  the  case  A ( v ) = 0 .In  fact,  K being  considered  an  arbitrary, 
fixed  constant  - easily  normalized  to  the  value  1 - in  Eq.  ( III  — 37 ) , the  h ^ (v  ) 
are  found  to  be  unbounded  and  thus  unrealizable  in  just  that  case  where  they  yield 
the  best  results.  The  output  signal  0 ( t 0 ) from  filters  with  "infinite"  gain  is  itself 
infinite.  Its  observation  would  thus  be  inconvenient.  We  will  see  in  Chapter  VI  how 
the  concept  of  normalized  matched  filtering  permits  us  to  avoid  this  difficulty 
in  the  case  of  a narrow  band. 
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CHAPTER  V 


THE  CASE  OF  IDENTICAL  SIGNALS  - PROPER  FILTERING-ORTHOGONAL  IMAGES 

OF  A SYSTEM  WITH  N NOISE  INPUTS 

Summary 

To  assume  N identical  signals  does  not  reduce  the  generality  of  the  problem. 

This  is  a case  to  which  we  may  always  return  and  one  to  which,  in  general,  even  the 
most  advanced  techniques  are  effectively  reduced.  In  this  case,  the  following  may  be 
established: 

a.  Matched  filtering  is  divided  into  (1  ) a proper  antenna 
filtering  (PF),  which  depends  only  on  the  N noises  and 
involves  one  filter  per  channel,  and  (2  ) a unique  filter 
defined  by  the  signal  alone,  which  is  nothing  more  than 
the  matched  filter  for  the  signal  in  noise  with  uniform 
density. 

b.  The  direct  sum  of  the  N elements  of  an  antenna  forms  a part  of 
a PF  - now  an  optimal  process  - only  if  the  N noises  are 
uncorrelated  and  have  the  same  spectral  density. 

The  concept  of  proper  filtering  is  very  important.  The  properties  of  a PF  are 
associated  with  the  noise  alone.  The  "system  function"  (complex  gain)  of  a PF  may  be 
defined  for  the  signal  between  the  input  and  output.  It  is  always  real  and  non-negative, 
and  the  first  characteristic  property  of  the  PF  is  the  following: 

The  gain  of  the  PF  is  equal  (except  for  a homogeneity  factor)  to  the 

spectral  density  of  the  output  noise. 

The  study  of  a PF  reduces  naturally  to  a study  of  the  correlation  matrix,  whose 
known  classical  properties  will  first  be  reviewed.  The  eigenvalues  of  this  matrix 
are  functions  of  frequency  and  have  the  properties  of  a spectral  density. 
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The  initial  argument  assuming  N correlated  noise  inputs  B j ( t ) may  be 
replaced  by  a simpler  argument  assuming  II  uncorrelated  noise  inputs  derived 
from  the  first,  which  represent  the  ensemble  B ^ ( t ) equally  well.  This  trans- 
position offers  the  same  convenience  as  the  reduction  of  a surface  to  principal  axes. 
Thus  it  is,  in  itself,  interesting  and  will  be  used  effectively  in  the  remainder  of  the 
study. 

The  N "reduced"  noise  inputs,  which  we  will  call  the  "orthogonal  images"  of 
the  B u ) , are  obtained  from  the  latter  with  the  aid  of  a collection  of  linear  filters 

cl 

comprising  a "matrix  a This  means  of  generation  illustrates  the  abstract 
operation  of  diagonalizing  the  correlation  matrix  C 


The  "orthogonal  images"  are  uncorrelated,  have  for  their  spectral  densities  the 
eigenvalues  of  the  correlation  matrix,  and  are  such  that  to  any  multiple  filtering  of 
the  B ( t ) there  corresponds  a multiple  filtering  of  the  images  that  yields  the 
same  signal  and  the  same  noise  spectral  density.  The  latter,  however,  is  easier 
to  study. 


In  particular,  to  the  matched  filtering  of  a signal  in  the  B , ( t ) there  corresponds 
the  matched  filtering  of  the  signal  transformed  (through  the  matrix  o ) to 
"orthogonal  image"  form.  Proper  filtering  of  the  B ( t ) corresponds  to  matched 
filtering  of  unlike  signals  in  uncorrelated  noise  . 


Thus,  from  now  on,  we  will  be  able  to  transpose  the  entire  study  into  the  realm 
of  "orthogonal  images",  which  will  permit  us  to  establish,  in  a simple,  "physical" 
way,  certain  properties  of  matched  filters  in  the  limiting  case  where  the  correlation 
matrix  becomes  singular. 


Finally,  the  second  characteristic  property  of  proper  filtering  is  established; 
it  is  a multiple  filtering  process  where  the  crosscorrelation  between  a single  noise 
input  and  the  output  noise  is  the  same  for  all  inputs.  Moreover,  this  mutual 
crosscorrelation  is  zero  for  any  value  of  t except  for  one  ( t = u ) but  dependent 
upon  the  choice  of  the  time  origin.  This  is  a "microscopic"  crosscorrelation. 

V-l.  Non-Restrictive  Character  of  the  Case  of  Identical  Signals 

We  will  consider,  for  the  present,  the  case  in  which  the  N signals  are  identical, 
first  emphasizing  the  practical  interest  of  this  special  case;  it  is  immediately  clear 
that  we  may  always  come  back  to  it.  Let  us  assume,  as  we  have  done  before,  that 
the  N inputs  are  the  N elements  of  a receiving  antenna  with  fixed  elements  arranged 
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in  some  arbitrary  manner.  The  signal  to  be  received  is  carried  by  a plane  wave  from 
a known  direction.  Let  us  put  on  each  element  an  ideal  delay  line  of  such  a value  that 
the  relative  time  lags  of  the  elements  for  that  particular  direction  are  equalized.  The 
output  signals  of  the  N delay  lines  become  identical.  Naturally  the  statistical 
relationships  between  the  noise  components  are  not  the  same  at  the  outputs  as  at  the 
inputs,  but  this  does  not  restrict  us  much  since  we  have  assumed  them  to  be  arbitrary. 
We  have  thus  constructed  a new  system  with  N inputs  (the  delay  line  outputs)  at  which 
the  N signals  are  identical.  The  operation  may  be  repeated  for  each  direction  of  the 
signal  plane  wave.  This  is  exactly  what  is  done,  in  general,  in  a technique  aptly 
named  "preformation  of  beams",  designed  to  provide  for  each  direction  an  output 
whose  signal-to-noise  ratio  is  improved  by  putting  all  the  signals  in  phase.  Let  us 
say  "improved"  and  not  "optimized"  since  it  will  be  seen  later  that  it  is  necessary  to 
consider  this  simple  "putting  in  phase"  and  the  limits  of  its  efficacy.  (See  Paragraph 
V-2.) 

Be  that  as  it  may,  making  signals  at  antenna  elements  identical  is  already  in 
current  practice  and  is  related  to  the  angular  separation  of  the  different  plane  waves 
which  are  received. 

Strictly  speaking,  the  process  of  delay  compensation  is  sufficient  to  render  the 
signals  identical  only  if  one  is  dealing  with  point  antenna  elements,  which  do  not  alter 
the  incident  sound  field  by  their  presence. 

If  there  exists  a rigid  structure  around  which  sound  waves  are  diffracted,  delays 
eomputed  from  the  geometry  are  not  adequate.  For  a given  element,  however,  this 
diffraction  appears  as  a linear  filtering  which  modifies,  at  the  input,  the  type  of 
signal  carried  by  the  incident  wave.  This  filtering  is  defined  by  the  relative  geometry 
of  the  element  and  the  structure  and  by  the  boundary  conditions  of  the  acoustic  field 
at  the  structure.  At  any  rate,  one  may  at  least  formally  assume  each  input  to  be 
affected  by  an  inverse  filtering  which  nullifies  the  diffraction  and  restores  identical 
signals  to  all  the  inputs.  The  compensation  process  is  complex  but  still  possible. 

Let  us  summarize: 

a.  At  the  formal  level,  we  may  always  reduce  the  problem  to  the 
case  of  identical  signals.  Properly  speaking,  there  is  no  restriction  here  upon  the 
general  study  of  matched  filtering  with  N inputs. 
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b.  From  a practical  standpoint,  it  is  often  very  close  to  the  situation  of 
N identical  signals. 

V-2.  Formal  Solution  for  the  Case  of  Identical  Signals 
Assume  that  the  N signals  are  identical: 


Si(t)  = s2(t)  = sM(t)  = S(t)  , 


(V-l) 


Si(v)  = S2(v)  = 


SN(v)  = S( v ) . 


Equation  ( III  — 37  ) is  then  reduced  to  the  following  N equations: 

-2 ni vt  o * -I  V' 

hk(v)  = Ke  • S <v>  . L-  MJk<»>  , 

J 


(V-2) 


for  k = 1,  2 N 


As  a result  of  (V-2  ) each  filter  h,  ( v ) is  obtained  by  placing  in  series  a 

4 K 

filter  whose  gain  is 

-2irivt  n * 

Ke  S (v)  , 

which  depends  only  on  the  signal  and  arbitrary  time  t q , with  a filter 


Pk(u)  * 


M^v) 


( V-3 ) 


which  depends  only  upon  the  noise  at  the  N inputs. 

-2irivto  * 

The  filtering  defined  by  Ke  S ( v ) is  nothing  more  than  the 

matched  filter  for  the  signal  S ( t ) in  noise  having  a uniform  density.  It  is  common 
to  all  N channels  of  the  system.  If  we  refer  to  Fig.  III-l  we  will  see  that,  in  the 


^Hereafter,  except  for  a homogeneity  factor,  the  expression  "complex  gain" 
is  understood. 
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case  of  identical  signals,  it  is  reduced  to  Fig.  V-l,  in  which  the  filtering  of  the 
signal  is  accomplished  after  summing  the  filtered  output  of  each  channel  p ( v ) . 

I 1 


I J 

(FILTERING  DEPENDENT  ON  NOISE  ONLY) 


Fig.  V-l.  Identical  Signals  - Arbitrary  Noise  Inputs 

Let  us  consider,  for  the  moment,  the  special  case  where  the  N noise  inputs  are 
uncorrelated  (see  Eq.  IV-2 ) . Equation  (V-2)  becomes 


hk(v)  = Ke 


-2TTivt 


S*  ( v ) 


(V-K) 


and  Eq.  (V-3)  becomes 


P*'V)  ' C"^T 


( V-5 ) 


If,  moreover,  the  N spectral  densities  are  identical  to  C ( v ) , all  the  p.  ( v ) 

1 K 

are  identical  to  gv  vy  . Figure  V-l  becomes  Fig.  V-l(a)  in  which  all  filtering 

processes  are  reduced  to  a single  filter  ~^ttivu0  g (v)  placed 

Ke  • ffTT 

after  the  direct  summation  of  the  inputs.  This  single  filtering  process  is,  moreover, 
the  matched  filter  for  an  arbitrary  input. 
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Fig.  V-l(a).  Identical  Signals  — Uncorrelated  Noise  Inputs 
Having  the  Same  Spectral  Density 

The  preceding  line  of  reasoning  presents  all  the  necessary  conditions  in  order 
that  the  direct  summation  of  antenna  elements  may  be  interposed  in  the  matched 
filtering  process  (see  Chapter  I).  As  attractive  as  the  replacement  of  N filters 
by  a single  one  may  be,  it  is  well  to  realize  that  an  optimal  process  may  be  realized 
in  this  wav  only  under  the  following  conditions: 

For  identical  signals,  with  uncor related  noise  inputs  having  the  same 
spectral  density,  the  matched  filtering  process  is  reduced  to  direct  summation  of 
inputs  followed  by  the  matched  filtering  for  an  arbitrary  input  signal, 

V-3.  Proper  Filtering  of  an  Antenna  and  Its  Gain.  First  Characteristic  Property 

We  will  now  define  the  important  concept  of  proper  antenna  filtering.  Let  us 
recall  the  general  Fqs.  (V-2)  and  (V-3)  , referring  to  identical  signals  and 
arbitrary  noise  (Fig.  V-l).  The  ensemble  of  filters  P^(v)  and  the  subsequent 
summation  will  be  called  the  proper  filtering  for  the  system  having  N inputs  in 
the  presence  of  noises  B j ( t ) , B ? ( t ) , ....  ( t ) 

The  letter  P designates  the  output  of  the  proper  filtering  process  where  the 
signal  will  be  o^Ct)  t °p(v)  and  where  the  noise  will  be  ?p(t) 

The  numerator  of  p,  ( v ) in  Eq.  ( V-3  ) is  the  sum  of  the  cofactors  of  the  terms  of 

t"  h ^ 

the  K 1 column  of  the  determinant  A(v)  . In  other  words,  it  is  a determinant 
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obtained  by  starting  with  A ( v ) and  replacing  the  terms  of  the  k column  by  ones. 

Proper  filtering  has  been  defined  from  the  concept  of  identical  signals,  that  is, 
by  a single  unique  signal  applied  to  all  elements  of  the  antenna.  For  such  a unique 
signal,  it  is  equivalent  to  a linear  filter  whose  system  function  (complex  gain)  is 

p (v)  since  the  spectrum  of  the  signal  at  P is  S(v)  *E  P^(v^  • 

k k 

The  system  function  G ( v ) of  the  proper  filtering  process  will  be  defined  as 

2 2 MIk(v)  / . (V-6) 

This  is  the  spectrum  of  the  response  obtained  at  point  P when  the  same  unit  impulse 
6(t)  (Dirac  delta  function),  whose  spectrum  S(  v)  = 1,  is  applied  at  the  inputs. 

The  ratio  Pm  (optimum  signal -to-noise  ratio)  at  the  output  becomes,  in  the 
case  of  identical  signals  (see  Eq.  ( III-39  ) ) , 


Gp(v) 


; ^ p,  (v  ) = 


P 


m 


( ls(v)|2  E E 

J j k 


MJk(v)  dv  , 


( V-7 ) 


or 


m 


G ( v ) dv  . 
P 


(V-7  a) 


The  density  D(v)  defined  by  ( II 1-40  ) becomes  in  this  case 

D ( v ) = | S ( v ) | 2 G (v)  , (V-8) 

and  as  a result  G ( v ) is  a real,  non-negative  function  that  is  even  in  v . 

Thus,  the  proper  gain  of  an  antenna  is  real  and  non-negative;  that  is, an  arbitrary 
spectral  component  of  the  input  signal  appears  at  point  P (Fig.  V-l)  multiplied  by  a 
real,  positive  number  and  therefore  is  not  changed  in  phase. 

The  matrix  notation  of  Paragraph  III-7  is  modified  because  column  matrix  S 
may  be  replaced  by  |^S  ( v )j  ot,  where  a is  the  column  matrix  whose  N elements 
are  equal  to  1.  Equation  ( II 1-46  ) becomes 

[-2 n i vt  o * _ 1 . , . „ x 

Ke  . S ( v ) J G a • 
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P95JI 


Equation  (III-50)  (spectral  density  at  53  ) becomes 


y y (v) 

2 o 


t2  2 I 1 

K | S ( v ) | a C"  a , 


from  which  we  get 

D ( v ) = £ | S ( v ) | J a C a 


(V-10) 


(V-ll) 


and 


G ( v ) = a C a 
P 


( V-12 ) 


The  spectrum  of  the  signal  at  output  P is  obviously  G ( v ) S ( v ) 

P 

but  at  the  output  53  , this  spectrum  is  (see  Eq.  (111-5^)) 
o(v)  = £s ( v ) j a h , 


or 


-2nivt0  2 

o(v)=  Ke  | S( v ) | 


1- 

Ja  C a , 


and,  of  course,  Eq.  ( II 1-5  3)  remains  valid  in  the  form 

o(v)  -•  ■ - ■ 2 ' -1 


yr  ( v ) 
1 o 


K 


-2  it i vt  i 


Ke 


a ( v ) 
p * 


(V-13) 


= D(v)  = | S ( v ) | "a  C"  a . (V-14) 


The  concept  of  proper  filtering  and  especially  of  antenna  gain  G ( v ) effectively 
expresses  the  properties  of  the  ensemble  of  N noises  at  the  n inputs.  This  is  a re- 
sult of  Eq.  (V-12  ) . The  spectral  density  of  the  noise  at  the  output  is  given  by  (V-10  ) , 
which  from  (V-12  ) may  be  written 


y_  (v)  = K2 |S( v) | 2 G ( v ) . 


(V-15) 
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b 


I 


Now,  this  density  is  the  product  of  the  spectral  density  at  P or  yd  ( v ) , 
with  the  square  of  the  magnitude  of  filtering  gain,  with  the  filter  placed  between  P 

and  S (Fig.  V-l).  It  may  be  seen  that 

YPo(v)  = Gp(v)  ' (V-16) 


Thus,  the  system  function  (complex  gain)  of  proper  filtering  is  equal  to  the 
spectral  density  of  the  noise  at  the  output  of  that  filtering.  (See  the  note  of  Paragraph 
V-2.)  This  property  is  nothing  more  than  the  characteristic  rule  of  matched  filtering 
(Paragraph  ni-7)  applied  to  proper  filtering. 


This  is  the  same  as  saying  that  everything  takes  place  as  if  the  antenna  were 
equivalent  to  a single  input  at  which  the  noise  has  a spectral  density  of 


A ( v ) 


1 

jvy 


A(v) 

^ E M . . ( v ) 

j k 


(v-17) 


the  gain  of  the  antenna  being  G ( v ) . and  the  input  signal,  S ( t ) . The  matched 

P 

filtering  of  such  a system  (see  Eq.  ( II-3)  ) is  then  comprised  of 


Ke 


-2  TT 1 vt  0 


* 

S ( v ) 


Gp(v) 


that  is,  precisely  as  in  Fig.  V-l,where  the  antenna  was  replaced  by  a single  filter  with 
a gain  G ( v ) . The  preceding  property  may  also  be  expressed  as  follows: 

Proper  filtering  is  such  that  the  spectral  density  of  the  noise  at  the  output  is 
equal  to  the  amplitude  spectrum  of  the  impulse  response  (signal  G ( t ) at  all  *; 
inputs) , which  is  nothing  more  than  the  system  function.  We  are  dealing  here 
with  an  intrinsic  property  of  the  antenna  in  the  presence  of  K given  noise  inputs. 

V-4.  Review  of  the  Known  Properties  of  the  Correlation  Matrix 

X 

Since  proper  filtering  of  a system  is  independent  of  the  signal  considered,  its 
study  is  reduced  to  the  study  of  a(  v)  and  its  cofactors  and,  more  generally,  of  the 
N statistically  related  noise  inputs.  Hence,  for  the  remainder  of  the  chapter,  we 
will  depart  somewhat  from  the  viewpoint  of  optimal  signal  detection  - to  which  we  will 
return  in  the  following  chapters  - in  order  to  review  the  mathematical  properties  of 
A ( v ) , and  we  will  relate  them,  by  their  physical  interpretation,  to  the  properties  of 
N statistically  related  noises. 
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Let  us  first  recall  that  each  term  of  the  determinant  A ( v ) is  a C . • ( v ) 

J 

having  Hermitian  symmetry  with  respect  to  v and  with  respect  to  its  indices,  ac- 
cording to  Eqs.  (II 1-5)  and  II 1-6): 


Cjk(v)  ' Cjk(  v)  ’ 


( I II -5 ) 


Cjk(v)  Ck,j^v)  • 


(III-6) 


As  a result,  the  determinant  A ( v ) itself,  being  a polynomial  in  C . , of  degree  II , 

j K 


has  Hermitian  symmetry  with  respect  to  v : 


, ( v ) = A ( -v ) . 


Furthermore,  thanks  to  Eq.  ( III-6  ) , reflection  about  its  principal  diagonal  leaves 
A ( v ) unchanged  and  equal  to  its  complex  conjugate.  Thus, 


A ( v ) = A ( v ) 


(V-19) 


As  a result,  A ( v ) is  real  and  an  even  function  of  v . We  know  that  it  is  also  non- 
negative, as  are  the  eigenvalues  of  the  matrix  C . 

In  the  following  discussion,  it  is  intended  to  define  exactly  how  this  property  ex- 
presses an  obvious  physical  fact;  namely,  the  sum  of  arbitrary  filtering  of  the  II 
noise  inputs  under  consideration  is  itself  a noise  and  by  virtue  of  this  fact  possesses  a 
necessarily  non-negative  spectral  density. 

Equation  (III-18)  , 


R.(u)  Rk(v)  C.R  (u-v)  du  dv  , (III-18) 


ps  ‘£  £ ff 

j k J J 

gives  the  noise  power  of  the  sum  of  the  arbitrary  filtering  R j ( t ) , R z ( t ) , ....  Rj,  ( t 

of  the  N noise  inputs  Bj(t),  B?(t),  ...  B,.(t)  (multiple  filtering) . 

This  relation  may  be  written  in  spectral  form.  Let  us  first  put  it  in  the  form 


P8  = E E 


Rk(v)  C , k ( u-v ) dv 


dU 


The  bracketed  term  is  the  convolution  of  R^r  ( t ) with  C . , . ( t ) . Its  F 
transform  is  then  Rk  ( v ) C . C v ) , and  it  may  be  written 


ourier 


f 


„ , s ^ / \ 2nivu 

Ri/V)  Cjk^v)  e 


d v 


or 


PB  = £ £ 


/ / 


Rk(v)  CJk(v) 


(„  / 2nivu  1 

Rjtu)  e j 


du  d v . 


* 

In  the  last  expression,  integration  in  u results  in  H . ( u ) , and  we  have  finally 


P 8 = 


/[ 


E E RJ(v)  Vv)  cik(v)  ] dv  • (v_20) 

j k -* 

We  would  also  have  been  able  to  obtain  Eq.  (V-2  0)  from  the  autocorrelation  of 
he  sum  of  the  filtering  processes,  that  is,  by  g moralizing  to  'J  noise  inputs  Eqs. 

( 1 1 1-9  ) and  ( III-10  ) , which  were  established  for  two  noise  inputs. 

Thus,  it  would  appear,  according  to  (V-2  0)  , that  the  spectral  density  of  the 
aultiple  filter  output  is 


yz(v)  = E E Rj^v)  Rk(v^  ci 

j k 


j'"  “kvv/  Cjkfv)  " R CR  » (V— 21) 


that  is,  the  Hermitian  form  of  the  correlation  matrix  C , It  is  non-negative  for  any 
multiple  filtering  R (positive  definite  matrix).  Let  us  recall  that,  in  a linear  trans- 
formation which  diagonalizes  the  matrix  (',  the  expression  of  the  Hermitian  form 
becomes 


. ( v ) = E * ^ ( v ) |x,(v) 


(7-.  2) 


where  the  X ( v ) are  the  eigenvalues  of  matrix  C - real  eigenvalues  since  C is 
Hermitian  - and  where  the  x j ( v ) correspond  to  the  F . ( v ) in  the  transformation 
(reference  (14),  Chapter  III). 

The  matrix  notation  for  (V-22)  is 


Yz( v)  = X *X  , 


( V - 2 3 ) 
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where  X is  the  column  matrix  of  X j ( v ) and  where  A (without  indices)  is  the 
diagonal  matrix  of  eigenvalues. 

Equation  ( V-2  3 ) is  true  for  any  arbitrary  X j ( v ) (resulting  from  any  arbi- 
trary R ( v ) . Thus,  it  is  evident  that 


A ( v ) >_  0 for  any  j and  v . 


(V-2U) 


The  eigenvalues  of  the  matrix  are  non-negative.  Their  product  is  non-negative, 


A ( v ) >0  , 


(V-25) 


all  these  properties  being  consequences  of  the  fact  that  Y j.  ( v ) for  an  arbi- 

trary multiple  filtering  process. 

Any  cofactor  M..(v)  plays  the  same  role  as  A(v)  for  N-l  noises  among 

J J 

the  N noises  considered.  Hence 


M j j ( v ) >0  . 


(V-26) 


More  generally,  every  minor  of  A ( v ) obtained  by  suppressing  a certain  number 
of  lines  and  columns  of  the  same  rank,  is  non-negative.  If  the  minor  is  of  order  two 

we  get  Eq.  (III-13)  . 


To  finish  this  review,  note  that 


Mjk(v)  ' Vi(v) 


(V-27) 


because  of  the  Hermitian  symmetry  of  A ( v ) . Thus,  the  sum  of  all  M , , ( v ) is 


We  have  seen  that  the  gain  for  proper  filtering  defined  by  (V-6  ) , 

M (v) 

Vv)  = £ E "irvr  * 


( V-6 ) 


j k 

is,  according  to  (V-8),  real,  even  in  v , and  non -negative.  Since  this  is  also 
the  case  for  A ( v ) , we  have 

YL  which  is  real,  non-negative  (V-28), 

j k and  even  in  v . 
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V-5.  Second  Characteristic  Property  of  Proper  Filtering 

Let  us  make  a final  observation  concerning  the  physical  interpretation  of  the 
equation  defining  proper  filtering,  which  may  be  written  (from  ( V-9  ) ) 


CP  = a , 


(V-29) 


which  summarizes  N equations  of  the  form 


C j 1 P1  + Cj  2P2  + 


C j N PN  ~ 1 » 

for  j = 1 , 2,  ...  N . 


(V-30) 


An  arbitrary  term  of  the  sum  ( V-30  ) is  in  the  form 

CJkP.  > 0r  ♦'  CJkPK  > 

where  <t  ( v ) designates  the  system  function  of  a filter  of  gain  1 over  the  whole  spec- 
trum. It  may  be  seen  by  analogy  with  Eq.  ( III-10  ) that  ^ (a  function  of  v) 
is  the  crosscorrelation  spectrum  0 of  two  noises,  the  first  resulting  from  the  filtering 
of  Bj(t)by  <t> ( v ) (that  is,  unfiltered  Bj(t))  and  the  second  resulting  from 
the  filtering  of  ( t ) by  ( v ) . 

The  sum  of  (V-30)  represents,  therefore,  the  crosscorrelation  of  Bj(t) 
with  the  sum  of 

B1(t)  filtered  by  Pj(v)  > 

B2(t)  filtered  by  P2(v),  etc,  ...  , 

in  other  words,  the  crosscorrelation  of  B(  t ) with  the  output  of  the  proper  filter. 


Thus,  Eq.  (V-29)  expresses  the  fact  that  the  crosscorrelation  spectrum  of  any 
one  of  the  noise  inputs  Bj  ( t ) with  the  output  of  the  filter  is  equal  to  1,  or  rather  to 
Q where  Q represents  a homogeneity  factor.  The  corresponding  crosscorrelation 
function  is  thus  in  the  form  Q 6 ( r ) . 

Hence,  proper  filtering  is  characterized  by  the  following  property:  it  is  a multiple 
filtering  process  such  that  the  crosscorrelation  between  any  one  input  and  the  output  is 


the  same  and  is  zero  for  all  t except  t = 0.  This  is  a "microscopic"  cross- 
correlation. 


That  is,  as  we  will  recall,  the  Fourier  transform  of  the  crosscorrelation  function. 
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In  the  special  case  of  a single  input,  the  proper  filter  has  a gain  proportional  to 

vhrr ' p(”)  ■ (v-3i: 

where  y ( v ) is  the  spectral  density  of  the  input  noise  (see  Eq.  ( IV- 1 ) ) . 

The  crosscorrelation  spectrum  of  noise  between  input  and  output  is  that  of  B(t ) un- 
filtered with  B ( t ) filtered  by  p ( v ) , 

or  according  to  (III-10),  [l]  [y(v)]  [p(v  )]  = 1.  (V-32 


The  corresponding  crosscorrelation  function  is  thus  equal  to  6 ( t ) except  for 
a homogeneity  factor. 
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CHAPTER  VI 


NARROW-BAND  APPROXIMATION  - VARIATION  OF  PROPER  FILTERING  (NPF) 


Summary 

This  chapter  describes  applications  and  develops  the  realization  ol  matched 
filters.  We  will  be  limited  here  to  the  important  and  practical  case  of  narrow-band 
signals.  It  will  be  assumed  that  in  the  frequency  band  of  interest,  properties  of  noise 
are  independent  of  frequency.  The  useful  gain  of  a linear  filter  is  reduced  to  a 
complex  number  that  represents  a "phase  shift"  accompanied  by  a "weighting"  of 
amplitude;  these  operations  are  always  realizable  in  the  sense  of  Network  Theory, 
provided  that  the  weighting  remains  bounded,  that  is,  provided  that  the  gains  remain 
finite. 

Thus,  the  narrow-band  approximation  offers  new  possibilities  of  normalization 
since  the  parameters  that  depend  upon  the  spectral  properties  of  the  noise  are 
reduced  to  constants.  These  possibilities  permit  us  to  avoid  a previously  described 
difficulty  (Paragraph  IV -5):  for  broad  bands,  the  gains  of  the  filters  are  not  finite 
(constant  K being  chosen)  if  the  correlation  matrix  becomes  singular,  that  is,  in  the 
limiting  case  where  the  signal -to-noise  ratio  might  be  "infinite."  For  a narrow 
band,  the  same  difficulty  may  occur,  but  it  is  avoidable  to  some  degree  by  using  a 
convenient  normalization.  Thus,  we  are  led  to  define  the  variations  ol  proper 
filtering  which  perform  in  the  same  way  for  all  ordinary  cases  but  offer  additional, 
useful  behavior  in  the  limiting  cases.  These  variations  are  characterized  by  the 
disappearance,  in  the  expression  for  filter  gains,  of  the  determinant  a(  v 0 ) 
of  the  correlation  mat r Lx. 

It  is  found,  moreover,  that  this  disappearance  provides  an  economic  advantage. 
The  technology  of  proper  filters,  as  will  be  outlined  in  Chapter  L\,  consists,  in  fact, 
of  constructing  the  C j and  then  the  M . , the  latter  being  polynomials  of  degree 

(N-l)  in  C.,  . But  A is,  itself,  of  degree  N.  By  avoiding  the  construction  of  A, 
we  reduce  the  technological  difficulty  by  one  degree. 

Here,  consequently,  are  the  points  which  are  to  be  considered.  Calling  tin' 
previously  defined  proper  filtering  simple  proper  filtering  (SPF),  let  us  consider 
a variation  of  it,  normalized  proper  filtering  (NPF).  which  has  the  following 
properties: 
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a.  It  is  always  defined,  even  when  the  correlation  matrix  becomes  singular, 
unlike  SPF. 

b.  In  this  latter  case  it  may  act  in  two  ways: 

- in  the  case  of  eliminable  noise,  the  noise  is  cancelled 
at  the  output  without  losing  the  signal. 

♦ 

- in  the  case  of  cutoff,  all  the  filter  gains  go  to  zero 
simultaneously. 

c.  A special  case  of  eliminable  noise  is  that  in  which  the  noise  goes  to  zero  at  a 
single  input;  all  channels  cut  off  except  the  single  input. 

d.  These  cases  are  physical  limits.  The  inevitable  presence  of  thermal  noise 
prevents  their  strict  realization. 

e.  Similarly,  it  will  be  established,  in  the  case  of  NPF,  that  no  change  occurs 
if  the  same  uncorrelated  noise  is  superimposed  upon  the  initial  noise  at  each  input. 

VI-1.  Narrow -Band  Hypothesis.  Physical  Meaning  of  the  Cj  ^ ( v 0 ) and  of  the 
Matrix  Notation 

We  have  seen  that  the  hypothesis  of  identical  signals  does  not,  strictly  speaking, 

constitute  a restriction  on  the  general  study  of  matched  filtering,  with  which  we  will 

be  concerned  from  now  on.  On  the  other  hand,  we  will  for  the  moment  consider  a 

G 

special  case  where  the  spectral  band  of  interest,  obviously  containing  the  signal  , 
is  such  that  the  statistical  properties  of  the  noise  (spectral  density  and  cross - 
cor relation  spectrum)  can  be  considered  as  being  independent  of  frequency.  Although, 
strictly  speaking,  such  a hypothesis  may  be  made  in  a band  that  is  arbitrarily  large, 
it  is  clear  that,  in  practice,  the  best  chance  of  realizing  the  hypothesis  lies  in  the  use 
of  narrow  bands.  Thus,  it  (the  hypothesis)  is  very  desirable  in  the  practical  cases  of 
radar  or  sonar  signals,  which  are  in  principle  narrow  band  - that  is,  of  a bandwidth 
that  is  small  compared  with  the  center  frequency. 

Let  us  assume  that  we  are  dealing  with  signals  lying  in  a narrow  band  centered  at 
vq  . That  is,  to  a first  approximation,  only  the  statistical  properties  of  noise  at 


Or,  more  precisely,  the  quasi-total  of  the  signal  energy,  since  a signal  of  finite 
duration  has  no  spectral  limits. 


the  frequency  v 0 will  enter  into  the  problem;  in  other  words,  we  are  interested  only 
in  the  values  C k ( v 0 ) . The  matched  filtering  is  thus  composed  of  the  filters 
nk  ^ vo  ^ distributed  over  the  • inputs  (Fig.  V-l),  then  the  summation  of  these  N 
filter  outputs,  followed  by  the  filter  S* ( v ) . This  last  filtering,  narrow  band  by 

hypothesis,  assures  to  some  extent  the  exclusiveness  of  the  frequency  v 0 

between  the  input  and  output  of  the  system.  The  filters  ( v 0 ) must  be  interpreted 
in  the  following  way:  sinusoidal  signal,  with  frequency  v0  , unit  amplitude,  and 
zero  phase,  denoted  by  a(t ) , applied  at  input  E,  , comes  out  with  amplitude 

I ( v 0 ) I and  phase  Arg  ( v 0 ) . Likewise,  the  proper  gain  G^  ( v 0 ) of  the 

system  represents  the  amplitude  and  phase  of  the  response  at  P for  the  same 
unitary  signal  a(t)  applied  at  the  N inputs.  Describing  all  the  transformations 
by  the  modulus  and  phase  obtained  under  the  above  condition,  all  equations  previously 
written  and,  in  particular,  the  matrix  equations  remain  valid  when  the  functions  of 
v are  replaced  by  their  values  at  v = vq  . For  example,  the  unitary  signal 
a(t ) applied  at  the  N inputs  may  be  described  by  the  same  column  matrix  a 
(all  elements  = 1)  that  previously  represented  the  spectrum  of  6 ( t ) . We  may 
write  then 

Gp(v0)  = aP  , 

where  p is  the  column  matrix  of  the  Pk  ( v 0 ) . We  will  thus  be  able  to  examine 
what  the  preceding  formalism  becomes  in  terms  of  the  narrow-band  approximation 
and  to  deduce  from  it,  should  the  occasion  arise,  its  special  properties. 

At  this  point,  we  have  no  need  to  isolate  the  narrow  band  by  means  of  filters 
since  S*(  v)  itself  takes  that  responsibility.  On  the  other  hand,  the  "filters" 

( vo  ) may  always  be  realized  since  they  consist  only  of  a "phase  shifter"  and  an 
"amplifier",  that  is,  very  simple  electronic  circuits  (provided,  however,  that 

I Pr  ^ v 0 ) I is  bounded).  *| 

Nevertheless,  if  we  had  to  construct  the  C k ( v 0 ) themselves,  we  would  need 
to  filter,  or  rather  prefilter,  the  N inputs  with  narrow-band  filters  <t>  centered  at  v p , 
which  are  identical  and  have  a gain  of  unity  in  the  narrow  band  and  zero  gain  elsewhere. 

If  the  B,  ( t ) were  put  into  these  filters  , the  outputs  would  be  narrow -band 

noise  b j ( t ) . We  have  seen  in  Paragraph  III-2  that  the  complex  number  C . f ( v 0 ) 
is  the  "spectral  line"  in  the  crosscorrelation  of  b , ( t ) with  t . ( t ) . In  fact,  the 
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narrower  the  filters  <J> , the  closer  the  crosscorrelation  function  is  to  a sinusoidal 
function  of  x with  frequency  v 0 . This  function  has  a phase  and  an  amplitude 
specified  exactly  by  C^(  v0  )(see  Eq.  III-ll)  . Thus,  technically,  it  is  the  cross- 
correlation of  b.(t)  and  b^(t)  which  C, k(v0)  will  represent,  and  this  same 
complex  number  defines  the  crosscorrelation  function  and  the  "crosscorrelation 
spectrum." 


Since  this  crosscorrelation  may  be  constructed  only  with  nonlinear  operations 
(products)  produced  by  interactions  between  spectral  bands,  preliminary  isolation 
(of  bands)  by  prefiltering  is  indispensable.  We  will  return  to  this  point  when  it  is 
required  to  consider  the  C^,{(\j0)  in  order  to  construct  the  pk ( v 0 ) (see  Chapter  IX) . 

VI-2.  Normalized  Proper  Filtering  - Technical  Advantages 


Multiple  filtering  (in  the  sense  which  we  have  outlined)  which  defines  proper  filtering 
of  the  system  (see  Fig.  V-l),  is  given  on  the  basis  of  Eq.  (V-3) : 


Pk(v0) 


MJk(v0) 


A ( v o ) 


(VI-1) 


NARR0W-3AND  SIGNAL 
WITH  CENTER  FREQUENCY  vq 


Fig.  Vl-1.  Normalized  Matched  Filtering 
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The  proper  gain  of  the  system  (see  Eq.  ( V— 6 ) ) becomes 


GpN0) 


Mj  R ( v o ) 


(VI-2) 


A ( V n ) 


We  have  seen  (Eqs.  (V-8)  and(  V-28  ) )that  G ( v 0 ) is  real  and  non -negative, 

r 

its  numerator  and  its  denominator  each  being  real  and  non -negative.  The  optimum 
signal-to-noise  ratio  is  (Eq.  (V-7  . a)  ) 


(VI-3) 


PM  Gp(vo)  • Ks  > 

where  E denotes  the  signal  energy  at  an  arbitrary  input.  We  know  that  the 

filters  in  a matched  filtering  process  are  definable  except  for  a constant  I'eal 
factor.  This  property  applies  here  also.  Thus,  the  filters  Py  (the  notation  vq 

will  be  understood)  of  Eq.(VI-l)  all  contain  a common  factor  i.  This  factor 

A 

depends  upon  the  properties  of  the  noise.  However,  in  all  cases  where  A is  not 
zero,  the  pk  may  be  replaced  by  the  filters 


"k  - S M.k  • 


(VI-4 ) 


Henceforth,  we  will  distinguish  between  simple  proper  filtering  (SPF),  defined 
by  Eqs.  (VI-1 ) and  (VI-2  ) , and  normalized  proper  filtering  (NPF),  defined 
by  Eq.  ( VI- 4 ) , to  which  corresponds  a normalized  proper  gain. 


r = 

p 

J k 

which  is  retd  and  non -negative. 


jk 


(VI-5) 


It  is  understood  that  these  limitations  apply  only  in  the  case  of  the  "narrow-band" 
approximation.  The  performance  of  the  NPF  remains  the  same  as  that  of  the  SPF 
and,  following  Eq.  (VI-3)  , depends  in  particular  upon  A: 


M 


E = G E 
s p s 


(VI-6) 


Now  we  must  determine  the  way  in  which  the  NPF,  as  defined  by  Eq,  ( V ! — • • ) , 
= 0 . This  will  be  the  objective  of  Paragraph  VI-5. 


behaves  when  A 


It  may  be  noted  that  the  NPF  presents  an  interesting  technological  advantage; 
although  Pn,  always  depends  upon  A,  the  filters  ( ) themselves  do  not 
depend  upon  it.  Thus  (see  Chapter  IX),  A is  a polynomial  of  order  N in  terms  of 
the  C . ^ whereas  the  are  only  of  order  N-l.  Computations  which  involve 
taking  the  product  of  complex  numbers  and  which  lead  to  the  "fabrication"  of  the 
tt,^.  have  a lower  degree  of  difficulty  than  those  required  for  the  computation  of  A . 

VI-3.  Invariance  of  the  NPF  in  the  Presence  of  the  Same  Noise  Superimposed  on 

E ach  of  the  B . ( t ) 

J 

A property  common  to  both  SPF  and  NPF  is  the  fact  that  the  phase  of  the  signals 
(narrow  band)  at  the  output  of  the  summation  of  the  N channels  is  unchanged.  This 
is  tied  in  with  the  fact  that  and  Gp  are  real  and  non -negative. 

A second  special  property  of  NPF  is  connected  with  the  structure  of  the 
Equation  (VI-4  ) shows  that  ^ is  the  sum  of  the  cofactors  of  a given  column, 
i.e.,  a determinant  of  order  N formed  from  A by  replacing  the  elements  of  the 
column  by  ones. 

It  is  easily  verified  that  if  the  same  constant  a is  added  to  each  of  the 
elements  of  such  a determinant  other  than  those  of  the  column,  the  value 

of  the  determinant,  and  hence  that  of  , is  not  changed.  The  physical  inter- 

pretation of  this  property  is  as  follows. 

An  NPF  is  unchanged  in  the  presence  of  a stationary  supplementary  noise  that 
is  identical  at  all  inputs,  not  correlated  with  the  B ( t ) , and  superimposed  upon 
them.  Let  us  assume,  in  fact,  that  the  same  noise  n ( t ) is  added  to  all  the  B ( t ) 
The  crosscorrelation  functions  C . ( t ) become 


NJk(T)  = E 


S j [Bj(t)  + n(fc)]  [Vt  + T)  + n(t  + T)]  J ’ 


which  is,  because  of  the  lack  of  correlation  between  n(  t ) and  the  other  noise, 

V’>  ’ CJk(t)  + Cn<’>  • 


where  C ( t ) is  the  autocorrelation  of  n ( t ) . 

n ' 


Hence,  the  Fourier  transform  is 


NJk(v)  = CJk(v)  + Cn(v)  , (VI-8) 

which  is,  for  a narrow  band, 

Njk(v0)=  Cjk(v0)  + Cn(v0)  . (VI-9) 

Thus,  all  the  Cj  k ( v 0 ) are  increased  by  the  same  amount  a = ( v 0 ) } 

which  is,  moreover,  the  spectral  density  of  n(t)  at  v = v0.  This  modification 
leaves  the  value  of  irk  unchanged;  hence  the  NPF  is  unchanged,  as  we  have  already 
seen.  It  is  evident,  however,  that  the  performance  (ratio  p ) _is  affected. 

Equation  (VI-6)  shows,  in  fact.that  if  is  unchanged  (Eq.  (VI-5)  ) , A 
itself  is  altered;  it  may  be  shown  that  it  increases,  becoming 

A + a rp 

and  that,  consequently,  as  would  be  expected,  the  signal-to-noise  ratio  decreases. 

It  may  be  said  that  everything  capable  of  constructing  the  n , from  the  C, 

K j k 

will  be,  a priori,  unaffected  by  a noise  n ( t ) at  the  inputs. 

An  analogous  reasoning  allows  us,  with  the  same  result,  to  replace  the  noise  n(t ) 
by  a sinusoidal  signal  of  "infinite"  duration  and  frequency  lying  in  the  narrow -band 
of  interest.  We  concede  that  under  quasi-stationary  conditions  (see  Paragraph  EX -2) 
with  a long,  real  signal,  it  may  be  assumed  that  the  tt  are  not  altered  bv  the 

K 

presence  of  the  signal  at  the  inputs. 

VI-4.  Matrix  Notation  for  NPF 

The  matrix  notation  for  SPF  and  NPF  is  derived  from  the  equations  of  Paragraph  V -3 
and  from  the  definition  of  the  column  matrix  tt  for  NPF. 

n = fA]  p . (VI-10) 

For  SPF  the  notation  is  that  of  Paragraph  V-3  , and  only  the  interpretation  of  the 

functions  of  v which  appear  in  it  is  changed,  following  the  convention  of  Paragraph  VI-1. 
For  NPF,  using  Eq.  (V-3)  , the  following  may  be  derived: 

tt  = [A]  C-1a  . (VI-11) 
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A 


Let  M be  the  adjoint  matrix  of  C (the  matrix  of  cofactors)  defined  by 


M = r a 1 c' 


(VI-12) 


it  = Ma 


The  gain  of  the  NPF  may  be  written  (see  Paragraph  VI -1  ) 


that  is, 


T = an  , 
p 


Ip  - a Ma  - £Aj  a C a = £a]  a P • 


(VI-13) 


(VI-14) 


(VI-15) 


These  equations  must  be  compared  with  Eq.  (VI-5)  in  order  to  be  defined.  The 
spectral  density  at  the  output  of  the  NPF  is 

Y = TT+  c TT  = [a2]  ct  c-1a  = [a]  a M a,  (VI-16) 


Y = A r 


( VI-l6a ) 


Having  replaced  an  SPF  by  an  NPF,  let  us  complete  the  system  by  adding  a 
filter  .S*(v)  in  series  with  the  summation  of  the  N channels.  If  the  narrow-band 
signal  s ( v ) is  applied  at  the  input,  the  spectrum  of  the  signal  after  summation  of 
the  N channels  is  S(  v ) . At  the  output  ^ (see  Fig.  VI-1)  it  is  I\  | f,(  v ) | 2 . 

Its  value  (maximum)  at  time  to  = 0 is  r E . The  noise  density,  which  has 

P b 

a value  y at  point  P , becomes  equal  to  Y I S ( v ) | 2 at  ^ *]  , and  its  power  becomes 


Figure  VI-1  summarizes  the  properties  of  an  NPF  followed  by  a filter  •'*(  \) 
(normalized  matched  filtering). 

VI-5.  Limiting  Behavior  of  the  NPF-Eliminable  Noise  and  Cutoff 

The  NPF  has  been  defined,  starting  with  the  SPF,  for  A ^ 0 . Let  us  assume 
it  to  be  realized  in  terms  of  the  C in  accordance  with  Eq.  ( V 1 — 1 ) . The  values 

' r ! 
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of  TT^  and  of  are  well  defined  and  remain  defined  if  A goes  to  zero7.  Then 
what  is  the  performance  of  the  NPF  under  these  circumstances?  We  have  already 
stressed  the  characteristics  of  the  "physical  limitations"  of  the  hypothesis  A = 0 . 

We  will  treat  it  here,  however,  in  a formal  way,  in  order  to  justify  the  use  of  the  NPF. 
In  fact,  note  that,  in  this  case,  the  filters  pk  (Eq.  (VI-1)  are  not  realizable 
since  they  are  not  defined.  This  is  then  a particularly  regrettable  case  of  failure  of 
the  SPF  (see  the  end  of  Paragraph  TV-G)  since  it  corresponds  to  an  "infinite"  signal- 
to-noise  ratio  (Eq.  ( VI -6  ) ) . On  the  other  hand,  the  NPF  is  realizable  and  can 
yield  this  "infinite"  S/N  ratio  provided  that  Tn  f 0. 

Such  performance  may  be  obtained  - it  j.  remaining  bounded  and  hence  the  output 
signal  remaining  bounded  only  by  nulling  the  noise  (while  the  SPF  produces  an  infinite 
signal  since  the  p are  unbounded) . This  is  clear  from  Eq.  (VI-16)  where  Y 
is  nulled  when  A = 0,  the  matrix  M and  the  CV  . themselves  remaining  bounded. 
Thus,  the  performance  of  the  NPF  is  better  than  that  of  the  SPF  in  the  case  where 
A = 0 ; it  is  defined  and  hence  it  provides  the  limiting  performance  of  matched 
filtering  by  nulling  the  noise  at  its  output. 

This  performance  exists,  in  fact,  only  if  itself  does  not  simultaneously 

go  to  zero,  that  is,  if  the  system  does  not  "cut  off",  nulling  both  the  noise  and  the 
signal.  This  would  be  the  case  in  particular  if  the  matrix  C were  of  rank  <N-2  . 

All  the  cofactors  M . , would  then  go  to  zero,  as  would  F . 

j k ’ p 

Thus,  the  "limiting  performance"  of  an  NPF  system  is  possible  only  if  the 
matrix  C is  of  rank  N-l  . 

If  the  matrix  C is  oi  rank  N-l , it  may  be  shown  that  r,  can  be  ze ro  only  if  A = 0 . 
An  NPF,  then,  has  two  possible  modes  of  behavior; 

a.  The  case  of  eliminable  interference  characterized  by  A = 0 and  r / 

P 

The  noise  is  nulled  and  the  signal  is  not. 

b.  The  case  of  cutoff,  characterized  by  r = 0 (hence  A = 0).  The 
multiple  filter  lets  nothing  pass,  by  letting  its  gain  go  to  zero.  It  may  be  shown  that 


7r 


The  woi'd  "becomes"  contradicts  "stationary."  Let  us  assume  a very  slow  change  in 
the  noise,  and  hence  of  the  C „ , allowing  time  for  adjustments  in  the  value  of  the 
and  concede  that  in  this  case  a situation  exists  where  = 
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such  a case  occurs  only  by  simultaneous  nulling  of  all  the  ti  , , that  is,  simultaneous 
cutoff  on  all  channels.  Returning  to  the  case  where  the  matrix  C is  of  order  N— 2 , 
it  is  obviously  still  a case  of  cutoff. 

VI-G.  Limiting  Case  of  a Zero  Noise  at  One  Input 

The  case  in  which  one  of  the  noise  inputs  is  found  to  be  zero  is  a special  case  of 
systems  with  eliminable  interference.  It  is  obviously  a limiting  case  and  merely  a 
formal  one,  since  in  such  a case  it  is  evident  that  one  need  only  use  that  single  input 
at  which  the  noise  is  zero.  If,  however,  one  of  the  noises  "becomes"  zero  (see  note 
of  Paragraph  (VI-5)),  let  us  verify  that  the  NPF  immediately  makes  use  of  that  fact. 

In  fact,  if  B^(t ) is  the  noise  that  goes  to  zero,  all  C go  to  zero  for  all  j . 

The  matrix  C , provided  with  a column  and  a line  of  zeros,  is  reduced  to  rank  N-l . 

The  representation  of  the  tt  . in  the  form  of  determinants  taken  from  A (by  replacing 
the  j -n  column  by  ones)  shows  that  all  of  them  go  to  zero  except  ^ , which  becomes 
equal  to  M.  , (real,  non-negative)  and  also,  in  this  case,  equal  to  r . Thus,  the 
only  channel  that  remains  open  is  the  one  in  which  the  noise  is  zero,  which  obviously 
results  in  zero  output  noise.  Note,  in  passing,  that  the  occurrence  of  two  input  noises 
B . ( t ) f being  zero, is  a case  of  cutoff,  all  the  it  . then  being  zero.  The  matrix  C , 
moreover,  is  of  rank  H-2  . Although  this  is  of  a particularly  formal  case,  we  see 
here  an  example  where  the  NPF  fails  because  of  an  excess,  so  to  speak,  of  good  noise 
behavior.  It  is  clear  that  we  may  get  an  output,  in  the  case  of  cutoff,  only  by  trying 
systems  having  fewer  than  N inputs.  (See  Paragraph  VII-7  on  the  optimal  use 
of  an  antenna.) 

VI-7.  Influence  ol  Independent  Noises  Superimposed  on  the  N Inputs. 

We  have  just  considered  some  limiting  cases  of  "eliminable  interference"  and 
"cutoff,"  all  of  which  involve  the  reduction  of  the  rank  of  the  matrix  C . Now,  let 
us  demonstrate  that  the  presence  of  unavoidable  thermal  noise  at  the  inputs  is  sufficient, 
in  practice,  to  exclude  these  limiting  cases. 

Thermal  noise  appears  as  noise  inputs  n . ( t ) produced  in  the  input  amplifiers  or 
in  the  receiving  elements  themselves,  superimposed  upon  the  B . ( t, ) , which  may  be 
considered  "external"  interference  noise. 


The  noises  n ( t ) are  characterized  by 

a.  non-zero  densities, 

b.  no  mutual  correlation,  and 

c.  no  correlation  between  themselves  and  the  B . ( t ) . 


In  this  case,  the  matrix  C ' of  such  a system  is 

for  j / k , 


cl,  = C,,  , 
Jk  jk* 


that  is, 


Cjj  Gjj  + dj  ; 


C + d , 


(VI-17) 


calling  d the  diagonal  matrix  made  up  of  the  spectral  densities  of  the  noises  n . ( t ) 
This  matrix  d is 

a.  of  rank  N since  no  d . is  zero,  and 

J 

b.  Hermitian  since  it  is  diagonal  and  real. 

Let  us  assume  that  the  rank  of  C * is  less  than  N ; this  means  that  there  exists 
a non-zero  column  matrix  (all  u . are  non-zero)  such  that 

C'u  = 0 , 

which  results  in 

U+C'u  = 0 ; (VI-18) 

taking  into  account  Eq.  (VI-17)  , Eq.  (VI-18)  may  be  written 

0 • (VI-19) 


u+Cu  + u+du 


The  two  bilinear  terms  of  ( VI-19  ) , being  non -negative,  must  each  go  to  zero 
separately.  The  column  matrix  u would  then  have  to  satisfy  the  following: 


u+du  = 


4J 


d,  = 0 


which  is  impossible  since  all  d are  / 0 and  all  the  u cannot  go  to  zero. 


<il 


rank  less  than  N . It  is  clear  that  this 


More  generally,  it  may  be  said  that,  if  some  portion  of  the  interference  noise 
uncorrelated  throughout  the  N inputs,  the  limiting  cases  will  not  occur. 


CHAPTER  VH 


TWO-INPUT  CASE  - COHERENT  NOISE 


Summary 

The  case  of  two  inputs  is  treated  here  in  a detailed  manner,  including  some  numer- 
ical calculations.  The  effectiveness  of  an  NPF  with  two  inputs  is  compared  with  the 
"direct  sum."  The  theoretical  gain  procured  from  the  NPF  is  calculated  as  a function 
of  two  parameters: 

a.  the  ratio  of  the  two  input  powers, and 

b.  the  complex  correlation  coefficient  @ (amplitude  and  phase  of  the  nor- 
malized correlation  function). 

This  gain  is  represented  by  a diagram  which  gives  the  curves  for  level  as  a function 
of  the  position  of  6 in  the  complex  plane.  Naturally,  it  is  always  greater  than  1 and 
is  "infinite"  in  the  case  of  eliminable  noise. 

Two  limiting  cases  are  settled,  one  by  an  equality  of  performance,  the  other,  the 
case  of  cutoff,  by  an  exceptional  inferiority  of  the  NPF  compared  with  the  direct  sum. 
The  latter  case  yields  the  same  result  as  a single  input. 

"Coherent  noise"  is  eliminable  from  the  inputs,  and  the  mechanism  for  its  elimi- 
nation is  considered. 


VH-1.  General  Equations 


Let  us  anply  the  preceding 

elements  C . , are  C i i and 
jk  1 1 

and  C i 2 = 


results  to  a system  with  two  inputs.  The  four  matrix 
C22  , spectral  densities  of  the  two  noise  inputs 
at  v = v o ; 

C-f  j , which  represents  the  crosscorrelation  of  the 
two  noises  at  v = vg  . 


Let  us  recall  the  relation 

2 

| Cl  2 I 1 Cl  1 C2  2 , 

which  expresses  the  fact  that  A of  matrix  C is  non-negative.  NPF  is  defined  by 

(VI-4). 

M = '2  2 “ '12 

TT  2 = ! 1 1 “ C2  1 . 


G3 


(VIT-1) 


Note  that  if  the  spectral  densities  are  equal  the  two  filtering  processes  are  com- 
plex conjugates.  Let  us  not  consider  here  the  case  where  one  of  noises  is  zero,  a 
case  which  comes  under  general  consideration  in  Paragraph  VI-7.  We  have 


r = TT  ! + TT2  = Cn  + C2  2 - C12  - C21  • (VII-2) 

Note  that  r is  nothing  more  than  the  spectral  density  ( at  v = v 0 ),  of  the 

r' 

difference  between  the  two  noises  (see  Fig.  IX-6).  In  fact,  the  crosscorrelation  of 
that  difference  is,  by  definition, 

E j [BjCt)  - B2(t)  ] [ Bj  ( t+r ) - B2  (t  + x)  ] J , 

or 

Cji  (t)  - Ci2(t)  - C2 1 ( t ) + C22(t)  , 

whose  Fourier  transform,  according  to  (VII-2),  is  equal  to  r^Cv)  . 

VH-2.  Eliminable  Interference  and  Cutoff. 

Let 


Cl  1 = 

5 1 

^ C1 1 

c2  2 > 

1 

j C i } J 

c22  = 

1 

q 1 

^ cl  1 

C2  2 » 

real,  positive  q = ^ 

c22  1 

(VII-3) 

c 1 2 - 

0 1 

^ C i i 

c22  , 

|e | < 1 (Eq.  (III-13) 

Thus  0 is  a complex  number  that  represents  the  amplitude  and  phase  oi  the  nor 
malized  correlation  function  for  two  noise  inputs  (or  degree  of  complex  correlation). 

The  matrix  G and  the  preceding  quantities  may  be  written 


* 


0 


I 


<;4 


yCjl  c22 


q 


(VII-4) 


Setting  aside  the  case  where  one  of  the  two  densities  is  zero,  systems  with  elim- 
inable  interference  are  defined  by 


0=1 


( VII-6 ) 


\ (q+|)  - R[o]  * 0 » 

where  R designates  the  real  part. 

We  always  have  ^ (q+jy)  — -*■  (VI I -7) 

Since  | 0 | =1  in  order  that  the  inequality  (VII-6)  be  realized,  it  is  necessary 

that  0 / 1 , that  is,  that  0 have  a non-zero  argument. 

In  more  technical  terms,  let  us  say  that,  if  two  noises  are  "totally  correlated" 

( | q | =i ) without  their  "average  phases"  being  equal  (0^1 ) , we  have  a system 
with  eliminable  interference. 

Thus,  conditions  for  cutoff  of  the  system  require,  according  to  the  preceding  re- 
marks, that  0 = 1 and  as  a result 


? ■ 1 


hence,  q = 1 


(VII-8) 


Summarizing,  we  have  0 = 1 

q = 1 .) 

Thus,  the  case  of  cutoff  may  occur  only  for  noises  having  equal  spectral  densities 


at 


vo  (q=l ) , 


Moreover,  the  condition  0 = 1 itself  reduces  Eq.  (VI  I -6)  to 


C.5 


A = 0,  or  rather  C12  = C21 


Ci!  C2  2 


r 

o 


1 


tt2 


1 


(VIT-S) 


If,  besides, 
nulling  tt  1 , 


the  two  densities  are  equal,  the  system  cuts  off  by  simultaneously 
tt  2 and  r , as  shown  by  Paragraph  VI-G  . 


Thus,  in  order  for  a two-input  NPF  to  cut  off,  it  is  required  that 


a.  the  spectral  densities  (and  consequently,  the  powers)  be  equal  and 

b.  the  "average  phases’1  be  equal,  that  is,  that 


l 2 


'2  1 


■{ 


1 1 


'2  2 
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1 1 
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Let  us  return  to  the  general  case.  The  signal -to-noise  ratio  is  (see  Eq.  (VI-6)) 


p 


m 


Cl  1 +C2  2-C 1 2-C2 1 
Cl  1 C 2 2 - I C ! 2 I 


• Es 


(VII-11) 


It  is  "infinite"  for  systems  with  eliminable  interference  and  indeterminate  in  the  case 
of  cutoff. 


Finally,  if  the  two  noises  are  uncorrelated,  C 1 2 = C2 1 

TT  ! = C 2 2 
r 2 = C1  1 
rp  = C 1 1 + C2-2 
Ci  1 +C  2 2 

= 77 7! 

rn  ^11  ^22  * 


and  we  have 


( VII-12 ) 


VII -3.  Advantages  of  Matched  Filtering  Over  the  Direct  Sum 

We  may  represent  the  "benefits"  derived  from  matched  filtering  in  the  following 

way: 

The  two  inputs  E,  and  E_  may  be  considered  to  be  receiving  elements  of  an 

I C. 

antenna.  Assume  that  the  signal,  identical  at  the  two  elements,  is  carried  by  a plane 
wave  whose  wave  front  is  parallel  to  the  straight  line  i H , . The  process  of 


(if. 


A 


|k 


"direction  summation"  is  most  commonly  used  to  enhance  the  signal-to-noise  ratio, 
which  means  that  proper  antenna  filtering  is  replaced  by  a simple  sum  (Fig.  VO-1). 

The  signal,  upon  summation,  has  for  an  amplitude  spectrum  2 S(  v)  . Subse- 

2 

quently  received  by  the  filter  S*  ( v ) , it  yields  an  output  signal  2 | S(  v ) | 
whose  value  at  time  1 0 = 0 is 


2 E 

s 


2 

and  whose  instantaneous  power  at  that  time  is  ^Es  . 


The  signal -to-noise  ratio  of  the  process  is  therefore 


C i j +C2  2+(-'  1 2+^2  1 ‘ 


(VII-15) 


The  improvement  given  by  matched  filtering  (or  NPF)  over  the  "direct  sum"  is 
expressed  by 


b = 


m 


(Ci  i +C 2 2_C i 2-C  2 l ) (Cj  I+C22+C1 2 +C  2 1 ^ 


^ Ci  1 C 2 2 ~ c 


1 2 


( VTI-16  ) 


The  word  "improvement"  is  justified  only  if  the  ratio  b >_  1 in  all  cases  except 
for  that  of  cutoff  of  the  NPF.  It  is  already  clear  that  no  improvement  is  gained  ( b = l ) 
when  the  two  noises  are  independent  and  have  the  same  spectral  density  (Cj  2=C2  i = 0 
and  C 1 i=C2  2 ) • We  have  seen,  in  fact,  that  in  this  case  only,  matched  filtering 
would  be  identical  to  direct  summation  (see  Paragraph  V-2). 

Let  us  make  it  obvious  that  b >_  1 . Taking  Eq.  (VTI-3)  into  account, 
the  ratio  may  be  written 

b = - [r(°>3 

1-  [RO)]Z  - [Ke)]2 


(VII-17) 


The  denominator  is  always  >_  0 , since  | 0 | < 1 . 

Since  i (q+^)^l  for  any  q, 
we  may  let 


[ 


\ (q4} 


= m +1 

with  m real  and  positive. 


(VII-7 ) 


(vir-18) 


Hence  , 


m2+l  - [K(0)]‘ 


m + [l(0)] 


b = 


= 1 + 


1-  [r ( © ) ] 2 - [1(0)] 2 1-  [r(o) ] 2 - [r(o)] 


<>8 


(VTI-19) 
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which  is  a form  from  which  it  is  obvious  that  b >_  1 . Uncorrelated  noises  ( 0=0  ) 
with  the  same  spectral  density  (m=0)  correspond  to  the  case  b = 1 (zero 
improvement) . 

The  case  | © | = 1 , which  makes  the  denominator  of  b equal  to  zero  and  gives 
an  "infinite"  improvement,  corresponds  to  the  case  of  eliminable  interference,  as 
we  have  alreay  seen. 

The  improvement  b remains  undetermined  in  two  cases, 
a.  The  first  case: 

m = 0,  0 = 1,  (VII-20) 


which  corresponds  to  the  case  of  cutoff  for  an  NPF.  The  signal-to-noise  ratio  for 
direct  summation,  given  by  Fq.  (VII-15),  may  be  written 


P 
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^ c l l C2  2 
2E 
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1 

\ (q+i)  + R(0) 


1 

1+m2  +R ( 0 ) 


(VII-21 ) 


and  is  not  zero  in  the  case  (VII-20)  . Thus,  the  indeterminacy  is  resolved  as 
a very  striking  inferiority  of  the  NPF  compared  with  direct  summation.  It  may  be 
said  that,  with  noise  inputs  having  the  same  spectral  densities  and  "in  phase" 

( 0 = 0 ) , the  NPF  cuts  off  while  direct  summation  produces  a signal-to-noise 

U1 

ratio  equal  to  S/„  , which  is  identical  to  that  produced  with  a single  input  and 

C11 


a filter  S*  ( v ) . Thus,  the  direct  sum  does  not  contribute  anything  itself,  but  at 
least  it  does  not  "spoil"  the  result.  Such  a case  is  found,  for  example,  when  the 
noise  carried  by  a plane  wave  comes  from  the  same  direction  as  the  signal  and  hence 

* 

is  identical  at  the  two  inputs.  In  this  case  alone,  and  provided  that  such  a noise  input 
is  indeed  the  only  interference  — otherwise  the  NPF  adjusts  itself  for  the  other 
noise  inputs,  according  to  Paragraph  VI-3  - the  NPF  cuts  off.  In  practice,  the 
presence  of  independent  thermal  noise  at  the  two  inputs,  combined  with  some  of  the 
considerations  of  Paragraph  VI-3,  prevents  cutoff,  even  if  the  condition  0=1 
does  not  represent  in  itself  the  limit  of  a limiting  condition  ( | 0 | = 1 ) . 


(i«J 


b.  The  second,  equally  idealized,  case  where  b is  undetermined  corresponds 
to 

m = 0 

0 = -1  . ( VII-22 ) 

Here,  the  NPF  behaves  as  in  the  case  of  eliminable  interference  and  gives  an 
’’infinite"  signal -to -noise  ratio.  It  is  found,  however,  that  the  direct  sum  does  the 
same  thing  since  the  two  noise  inputs  are  in  phase  opposition.  Note  that  the  two 
filters  7t  j and  tt2  are  now  real  and  identical  (see  Eq.  (VII-5)  with 
q = 1 and  0 = 1)  and,  because  of  this  fact,  produce  the  sum  precisely,  except 
for  a multiplicative  factor.  Thus,  the  indeterminacy  is  resolved  as  an  equality  of 
performance. 

VII-4.  A System  of  Graphs  Representing  the  Gain  of  the  NPF 

Let  us  set  aside  the  two  extreme  cases  examined  above.  The  improvement  b , 
always  greater  than  unity,  may  be  represented  for  each  value  of  m,by  a surface 
(Eq.  (VII- 19  ) ) above  the  complex  plane  for  O . This  surface  is  entirely  contained 
in  the  vertical  cylinder  of  radius  1 and  is  tangent  to  it  at  infinity  (the  case  of 
eliminable  interference).  Furthermore,  this  surface  is  symmetric  with  respect  to 
the  two  planes  1(0)  = 0 and  Arg  0 = , and  it  is  sufficient  to 

represent  it  within  a 90  degree  dihedron  (in  a single  quadrant).  It  may  be  represented 
there  by  level  curves.  Figure  VII-2  shows,  in  four  successive  quarter  planes,  the 
level  curves  of  the  surfaces  for  four  values  of  m , that  is,  four  values  of  the  ratio 
of  the  spectral  densities  1 l l / t These  four  values  are 
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The  greater  the  ratio  C1  1/  , the  higher  the  "base  point"  of  the  surface 

C22 

, 2 

(corresponding  to  0 = 0 , b 0 = 1 + m ) and  the  more  rapidly  the 

surface  rises,  as  I 0 I increase,  toward  the  "infinite"  values  which  correspond  to 

I 0 I = 1 . The  level  curves  are  ellipses  whose  major  axes  lie  along  the  real  axis. 
These  ellipses  tend  toward  circles  when  b becomes  very  large  . We  may  say,  then, 

c 

that  regardless  of  the  value  of  * 1/  the  more  "advantageous"  of  two  values 

02  2 

of  0 having  the  same  modulus  is  the  least  real  one. 

The  critical  zone  corresponding  to  m=0  and  0 = 1 has  been  crossha tched 
in  Fig.  VII-2  . 

Vn-5.  Coherent  Interference  (Jamming)  - Mechanism  for  its  Elimination 

We  will  call  interference  noise  coherent  if  it  is  carried  by  a plane  wave  coming 
from  a direction  other  than  that  of  the  signal.  This  is  the  important  and  practical 
case  of  a jammer  or  of  a localized  distant  source.  Two  noise  inputs 
B 1 ( t ) and  B 2 ( t ) differ  from  each  other  only  by  a time  delay. 

We  will  generalize  the  definition  of  coherent  interference  a little  and  assume  that 


B2(t)  a B j ( t-u ) , (VII— < 

where  a is  a real  multiplicative  factor  and  u is  the  time  lag.  It  is  then  simple  to 
verify  that 

C l 2 ( t ) = a Cii(t-u)  , (VII-21) 

and  as  a result. 


C12(v)  = a Cji(v)  e ‘ 11  ^ v u 


(VII-25) 


Since  we  have,  furthermore, 

C 2 2 ( v ) 


2 

a C 1 1 ( v ) 


(VII-26) 
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we  may  write,  for  frequency  v = v0, 


C 1 2 ( v o ) = e 2"1V°U  ^ Cu(v0)  . C22(v0)  . (VII-27) 


The  number  0 defined  by  Eq.  (VI 1-3)  is 


g _ e-2niv0u  (VII-28) 

Thus,  we  have  the  case  (q  j = 1 which  gives  a system  with  eliminable  interference. 
Hence,  coherent  interference  is  eliminable  interference  in  a system  with  two  inputs. 

Figure  VH-3  illustrates  how  the  noise  is  eliminated.  In  fact,  the  two  filtering 


processes  tt  j and  n2  are,  except  for  the  factor  \l  C i i C22 
(see  Eq.  (VII-5)  ), 


tt  i = a - e 


-2  tt  i v o u 


n2  = - + e 


+2ni v 0U 


(VII-29) 


Each  may  be  considered  to  be  the  sum  of  two  filters  in  parallel.  The  term 

e ' represents  a delay  u,  and  its  conjugate  represents  an  "advance" 

In  the  general  sum  of  the  four  terms,  noise  inputs  cancel  each  other  in  pairs,  but  two 
opposing  terms  come  from  different  inputs. 


Likewise,  it  is  easily  seen  how  the  case  of  cutoff  may  occur,  when  simultaneously 


a = 1 
u = 0 


that  is,  when  the  two  noise  inputs  are  identical,  or  rather  when  the  coherent  interference 
comes  from  the  same  direction  as  the  signal. 

The  four  terms  of  the  noise  destructively  interfere  in  pairs  on  the  same  input. 

Thus  a two-input  system  may  eliminate  coherent  interference  provided  that  the 
interference  does  not  come  from  the  same  direction  as  the  sigmil.  It  is  clear,  in  fact, 
that  in  this  latter  case,  no  spatial  discrimination  between  signal  and  noise  being  possible, 
one  may  do  no  better  than  when  using  a single  input;  direct  summation  itself  results  in 
no  improvement  (Paragraph  VU-3a), 
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Fig.  VII-3 


VII -6.  Return  to  a Remark  Made  in  Chapter  I - The  Special  Case  of  Coherent  Interference 

The  reasoning  illustrated  by  Fig.  VII-3  does  not  necessarily  use  the  "narrow-band" 
hypothesis;  in  fact,  coherent  interference  is  interference  that  is  eliminable  over  a 
broad  band. 

Figure  VII -4  shows  a simpler  method  of  effecting  this  cancellation  in  the  case 
where  a = 1 . By  conveniently  delaying  the  noise  at  one  input  and  taking  a 
difference,  the  noise  is  effectively  cancelled  without  cancelling  the  signal,  if  the 
signal  comes  from  another  direction.  This  simple  case  was  brought  up  in  Chapter  I 
to  suggest  that  direct  summation  is  not  always  the  best  method. 
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Fig.  VII -4 

VII-7.  Optimum  Use  of  an  Antenna 

A few  of  the  preceding  remarks  may  be  very  quickly  extended  to  the  case  of 
3 inputs,  permitting  their  generalization  to  N inputs.  It  is  easy  to  verify,  for 
example,  that  "coherent  interference,"  which  is  oliminable  interference  with  two  inputs, 

g 

results  in  cutoff  with  three  or  more  inputs. 

g 

Each  of  the  N noise  inputs  are  derived  from  a reference  noise  by  (Eq.  (VII -2  3)) 
b1 (t)  = a ^ b o (t-u. ) , 

°r  Cjk(v)  = aj  ak  Co(v)  e-2,Tlv(uj-uk)  « 

where  C n ( v ) is  the  spectral  density  of  b 0 ( t ) . After  putting  them  in  convenient  form, 
it  may  be  proven  that  the  determinant  of  matrix  C is  reduced  to  a determinant  of  order  K , 
all  of  whose  elements  are  equal  to  1.  Thus  the  rank  of  the  system  is  1 in  all  cases. 
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To  eliminate  a jammer,  for  example,  assuming  that  it  is  the  only  interference, 
a system  of  only  two  inputs  would  be  required. 

Thus  with  different  N ’ s , various  possibilities  are  offered  for  eliminating 
interference,  but  these  possibilities  depend  upon  the  value  of  IJ  . 

On  the  other  hand,  the  number  of  possibilities  is  greater  for  larger  N 
(number  of  solutions  to  A = 0 ) . It  is  indicated  that,  in  order  to  derive 
the  most  advantage  from  an  antenna  with  N elements,  several  systems  should 
be  used; 

a.  a system  with  two  elements  which  are  used  directly, 

b.  a system  of  3 elements  which  may  include  the  two  preceding  elements,  and 

c.  then,  a system  with  4,  5 ...  N elements, 

the  last  making  use  of  the  entire  antenna.  We  have  at  our  disposal,  then,  all 
possibilities  offered  by  the  antenna  from  the  point  of  view  of  "eliminable  interference." 

For  example,  since  coherent  interference  will  be  eliminable  on  a system  with 
two  elements,  it  doesn't  matter  that  it  causes  cutoff  on  the  others.  Other  noise 
configurations  will  be  eliminable  with  3 elements,  etc. 


. 
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CHAPTER  VIII 


MATCHED  FILTERING  AND  DIRECTIVITY 

Summary 

After  a review  of  the  definition  and  physical  meaning  of  antenna  directivity,  as  well 
as  the  definition  of  omnidirectional  noise,  this  chapter  establishes  the  fact  that 
optimization  of  directivity  is  only  a special  case  of  matched  filtering  applied  to 
special  crosscorrelation  properties  created  at  the  antenna  by  omnidirectional  noise. 

Thus,  theories  and  processes  tending  to  optimize  directivity  optimize  signal 
detection  only  to  the  degree  that  the  more  or  less  implicit  hypothesis  of  omnidirec- 
tional noise  is  valid. 

A few  special  properties  of  plane  or  linear  antennas  are  reviewed. 

VIII-1.  Limits  of  the  Validity  of  the  Directivity  Concept 


Current  usage  in  technical  literature  assumes  that  the  way  to  favor  detection  of  a 
signal  carried  by  a plane  wave  is  to  obtain,  by  means  of  an  antenna,  "the  greatest 
directivity"  in  the  direction  of  that  wave.  We  have  just  seen,  however,  that  optimization 
of  detection  depends  upon  the  statistical  relationships  of  the  interference  noise  inputs 
to  the  antenna  elements.  If  we  claim  to  optimize  detection  by  optimizing  "directivity," 
we  become  involved  with  a more  or  less  implicit  hypothesis  concerning  the  nature 
of  the  statistical  relationships.  We  will  see,  in  fact,  that  this  hypothesis  is  that  of 
"omnidirectional  noise,"  which  corresponds  to  an  exact  definition  (as  does  coherent 
interference)  and,  consequently,  to  a well-defined  correlation  between  two  elements. 
Likewise,  we  will  see  that,  as  would  be  expected,  to  optimize  directivity  is  nothing 
more  than  performing  matched  filtering  in  the  presence  of  omnidirectional  noise,  that 
is,  treating  a special  case  of  the  spatial  structure  of  the  noise.  Thus,  from  this 
chapter  we  may  derive  ways  to  modify  some  of  the  hasty  conclusions  regarding 
the  general  character  of  this  directivity  concept,  to  which  a considerable  theoretical 
and  technical  effort  has  been  attached.  To  assume  that  a noise  is  omnidirectional 
when  we  know  nothing  about  it  is  a simple  hypothesis  but  it  is  not  an  optimum  procedure. 
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VTH-2.  Definition  and  Physical  Meaning  of  Directivity 


First,  let  us  recall  the  definition  of  directivity.  The  "complex  directivity"  is  a 
function  of  frequency  and  of  two  spatial  parameters  defining  the  direction  of  a plane 
wave.  Let  us  combine  these  two  spatial  parameters  into  the  symbol  w.  The 
complex  function  D ( u , \> ) gives  the  phase  and  amplitude  of  the  voltage  observed 
at  the  antenna  output  in  the  presence  of  a plane  wave  having  direction  u and 
frequency  v . The  output  is  defined  as  the  sum  of  the  voltages  coming  from  the 
antenna  elements  through,  if  need  be,  a linear  filter  in  each  channel  (multiple 
filtering) . 

The  quantity 

4u | D ( u o , v ) | 

f(v)  = (vrii-i) 

J" | D ( u , v ) | 2 dw 

is  called  the  directivity  factor  relative  to  the  direction  Mo  (or,  in  shortened,  but 
less  correct  form,  "directivity"),  where  dm  represents  an  elementary  solid  angle 
about  direction  <d.  To  optimize  the  directivity  is,  in  fact,  to  optimize  f(v)  . 

The  physical  meaning  of  this  quantity  is  obvious.  The  square  of  the  magnitude 
of  D(  u),  v ) has  the  meaning  of  power  per  unit  solid  angle.  The  denominator  of 
f ( v ) is  thus  the  total  power  of  independent  contributions,  each  assigned  a 
direction  w . We  may  still  say  that  it  is  the  power  at  the  output  of  an  antenna  when 
the  antenna  is  in  a field  of  plane  waves  coming  from  all  directions,  which  are 
independent  of  direction  (even  infinitesimally  separated)  and  have  the  same  spectral 
density9.  This,  then,  is  the  precise  definition  of  omnidirectional  interference . 

The  numerator  of  f ( v ) may  be  said  to  be  the  power  received  by  the  antenna 
in  omnidirectional  noise  if  the  antenna  possesses  the  same  response  for  all 
directions  that  it  has  for  the  direction  w0  • Consequently,  to  say  that  "the  directivity 
is  high"  (or  f ( v)  is  high)  is  the  same  as  saying  that  the  antenna  is  particularly  suited 


9That  is,  "perfectly  diffuse  radiation." 
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to  "avoid"  plane  waves  coming  from  directions  other  than  coq  , thereby  achieving 
a spatial  selection;  this  has  the  form  of  a linear  filtering  process  with  a gain  that  is 
a function  of  direction.  We  may  also  say  that  the  numerator  is,  except  for  some 
factor,  the  power  produced  at  the  output  by  a plane  wave  coming  from  direction  up  ■ 
Thus,  f(v)  also  represents  the  signal -to-noise  ratio  (ratio  of  mean  powers) 
at  the  output  for  a signal  plane  wave  from  direction  wo  and  an  omnidirectional  inter- 
ference noise  field. 


The  definition  of  f ( v ) , then,  is  strictly  related  to  that  of  omnidirectional 
noise,  and  the  use  of  the  concept  of  "directivity  factor"  implies  the  omnidirectional 
noise  hypothesis.  The  value  of  f ( v)  depends  upon  antenna  geometry,  which  is 
assumed  here  to  be  defined  by  the  N elements  and  upon  the  filtering,  if  required, 
in  each  channel.  It  is  the  optimization  of  f ( v ) by  means  of  this  multiple  filtering 
which  we  will  discuss  here. 


Vffl-3.  Expression  for  the  "Directivity  Factor"  as  a Function  of  the  Gains  of  the 
Multiple  Filters 

Now,  let  us  consider  a collection  of  H antenna  elements  ,E2,  ...  E^ 

(Fig.  Vm-1),  all  identical  as  far  as  electromagnetic  or  electroacoustic  receiving 
elements  are  concerned,  and  a phase  reference  point  0_. 


For  a plane  wave  propagating  in  direction  w 
amplitude,  the  voltage  received  by  element  Ej 
amplitude  and  phase  by 


with  frequency  v and  unit 
is  represented  in  terms  of 


2tt1U  (OE.  . w) 
e C J 


> 


where  w is  the  unit  vector  in  direction  w . The  scalar  quantity  OE.  . w 

represents  the  path  difference  of  the  plane  wave  between  and  E . ; when 
divided  by  c (propagation  velocity)  it  is  the  "time  lag"  or  "advance"  of  E , 
with  respect  to  0 . 

Let  us  begin  by  making  "identical  signals”  by  applying  to  each  input  a filtering 
process  which  compensates  for  the  path  differences  among  the  different  inputs  for  a 
plane  wave  coming  from  direction  wo  • Hence,  we  are  led  to  apply  to  I . a 
" filtering  process" 
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which,  over  the  frequencies  v , is  merely  an  ideal  algebraic  "time  lag" 

oSj  . w 
C 

Next  let  us  apply  to  each  input  a supplementary  filtering  q ( v ) , which  is 

unspecified  for  the  time  being  and  which  will  be  varied  in  order  to  optimize  the 
directivity  factor  f ( v ) . 


The  voltage  at  the  output  of  the  filter  q . ( v ) in  channel  E for  a plane 

J j 

wave  coming  from  direction  w is  represented  by 

[OE  . C-:o)] 
q . ( v ) e J 

Thus,  the  previously  defined  "complex  directivity  function"  D(u,v)  is 
(see  reference  (13)) 


D(w, v ) 


2 1“  [Of  . (Z-Z0)] 
2 Qj  ( v ) e c J 


(VIII-2) 


and  we  have 


D(  wg  , v ) = 2 Qj(v) 


(VIII-3) 


The  numerator  of  the  directivity  factor,  according  to  Eq.  (VIII-1)  , is 


47t|D(w  ,v)  | = 4rr|  2 Qj(v)|2  . 


(VIIT-4) 


The  denominator  is 


I = f I D ( , v)  | 2 dw  - 2 j.  Qj(v)  e 

J j k 


* / 2 l^[(OE  -OEk).  («-«„)] 


doj. 

( VIT-5 ) 


Let 


VCv)  = inr/ 


2ni^-  [(OE.-OEk).(w-w0'] 


dw 


(VITI-6) 


with  the  result 


1 = ^ 2 2 Qj(v)  qk(v)  V j R ( v ) 

j k 


( VIII-7 ) 


Note  that 


Vjk(-u)  ■ V<i,(v) 


jk 


VkJ(v)  =VJk(v) 

lvJkCv)|  < 1 


CVIII-8 ) 


VJj(“>  ' 1 

On  the  other  hand,  I is  real  and  non-negative  for  all  q . . The  expression 

J 

for  the  directivity  factor  is 


f(v)  = 


2 q1 ( v)  | 

J J 


( VIII-9 ) 


2 2 qj(v)  qk(v)  Vjk(v) 

j k 


and  the  problem  consists  of  finding  the  q^(v)  that  maximize  f(v)  , 

Vni-4.  Optimization  of  the  Directivity  Factor 

The  solution  has  an  obvious  relationship  to  the  reasoning  of  Chapter  III  and  is 
purposely  presented  here  in  an  analogous  and  brief  manner.  We  will  temporarily 
omit  the  frequency  notation  ( v ) , which  will  remain  understood. 

Equation  (VIII-9)  defines  f(v)  for  the  q , , which  are  themselves  defined 
except  for  a complex  factor.  This  factor  may  be  chosen  such  that  the  numerator  of 
f(  v ) is  normalized;  that  is, 


v 

J 


qJ  = 1 » 


(VIII-10) 
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which  is  the  condition  under  which  the  will  be  constrained.  Consider  complex 

numbers  r . for  which 

3 


2 r.  = 0 . (VIII-11) 

J 


All  numbers  xj  = q ^ + ar^ 
condition  for  the  numerator 


with  real  arbitrary  a will  satisfy  the  normalization 


2 

j 


Let  hj , 

(VIII-10) . 


in  particular,  be  a solution  sought  among  the  q j 

This  solution  must  minimize  the  value  of  I . ; 

4 7T 


I_ 
4 IT 


2 2 

j k 


* 

q>.  v 


Jk 


constrained  by 


( VIII-7 ) 


Thus,  in  particular,  the  value  of  I 

for  q*  = h.  + ar  for  any  r 

J J J J 

inequality 


for  Qj  = is  less  than  any  value  of  I 

constrainea  by  (VIII-11)  • We  may  write  the 


L < L + ar 
n — h 


(VIII-12 ) 


which  may  be  expressed  as 


2 2 hj  hk  Vjk  - 2 2 (hj+otrj)  (hk  + ark)  Vjk 

J k J k 


or 

2 


(VIII-13) 

£ J rJ  rk  Vjk  + “ 1 2 (rj  hk  + hy  rk>  VJk  i 0 • 

J k j k (VI 11-14 ) 

2 

The  coefficient  of  a is  real  and  non-negative,  since  it  is  a particular  value 


of  for  4.  = r,  . In  order  that  the  inequality  (VIII-14)  be  true 

' M 7T  J 

for  any  a , we  must  not  allow  the  sign  of  the  term  in  a to  dominate  for 


small  values  of  | a | . The  coefficient  of  a , which  is  divided  into  two  complex 
conjugate  terms,  must  therefore  go  to  zero.  Setting  one  of  the  terms  equal  to  zero 
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we  get  the  condition 


2 2 rj  nk  Vjk  “ 0 

j k 


(VIII-15) 


for  all  r . constrained  by  (VIII-11)  . Equation  (VIII-15)  may  be  written 


h,,  V ..  i = 3 


3 "J  [2  VJk] 

j k 


(VIII-16) 


Comparing  this  with  Eq.  (VIII-11)  , it  is  clear  that  all  the  coefficients 
of  the  r . on  (VIII-16 ) must  be  equal  to  the  same  number  y . The  solution, 


then,  is  given  by 


* 

2 hk  Vjk  = J = 1>2, N.  (VIII-17) 

k 


The  value  of  I , , minimized  by  the  solution  h ^is  (see  Paragraph  VIII-7) 


2 2 

.1  k 


h.  h.  V.. 

J k jk 


'3  iz 

k 


hk  Vjk]  • (VIIT-18) 


hj  = W 


taking  (VIII-17  ) into  account  and  the  fact  that  h ^ is  constrained,  as  are  all  the 
Qj  (Eq.  (VIII-10)  ) . 

The  directivity  factor  f , made  a maximum  by  the  solution  h , is 


f _ *1  IT 

m = T~ 

m 


(VIII-l'i) 


The  solution  (VIII-17)  written  in  matrix  form  is 


(VITT-20) 
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where  V is  the  matrix  whose  elements  are  the  V . , (Hermitian  with  respect  to  the 

J K 

indices) , h is  the  column  matrix  of  the  h A , and  a is  the  column  matrix  all  of 


j 


whose  elements  are  ones. 


Vm-5.  Correlation  Matrix  for  Omnidirectional  Noise 

We  have  just  calculated  the  "filtering"  that  optimizes  the  directivity  of  an  ensemble 
of  N antenna  elements.  The  proper  filtering  corresponding  to  omnidirectional  noise 
will  now  be  calculated,  in  order  to  verify  the  fact  that  it  is  identical  to  the  filtering 
process  h defined  by  (VIII-20).  The  similarity  between  Eqs.  (VI 1 1-20  ) and 
( 111-45  ) (or  rather  ( V-3  ) , since  it  is  a case  of  proper  filtering)  clearly  shows 
that  we  need  only  prove  that  the  matrix  V*  is  the  matrix  of  crosscorrelation  spectra 

( C . , ( v ) ) in  the  case  of  omnidirectional  noise  after  the  signal  is  made  identical 

J k 

at  the  N inputs  . 


Thus,  the  problem  is  reduced  to 

a.  evaluating  the  crosscorrelation  between  two  inputs  E.  and  E^, 

b.  taking  into  account  the  alteration  of  that  crosscorrelation  by 
making  the  signals  at  Ej  and  identical,  and 

c.  confirming  that  the  result  is  V.  . 

J K 

The  first  step  is  a classical  calculation  which  already  has  been  done  for  some 
special  cases  by  B.  Picinbono1  and  is  reproduced  here  using  a general  notation 
compatible  with  the  rest  of  the  paper. 


Let  us  call 


b0“(t) 


the  elementary  noise  received  at  point  0 from  direction  w 


(see  Fig.  VIII-1).  The  noise  received  at  E from  the  same  direction  is 


j 


bj  (t)  = b0“(t  + 


OEj  . w 


) . 


The  corresponding  noise  at  e,  is  b(.  ( t ) = bQU  (t  + 


OE^ . u> 


( VIII-21 ) 
) . (VIII-22) 


10 


Not  published. 
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The  elementary  crosscorrelation  (for  direction  u>  )is  then 


- E 


b “(t)-b,  "(t+T 


)|  = m0u'  ( 


t + 


OE,  -0E, ) • w 

— J ) ( VIII-2 3 ) 


designating  by  mgW(  t ) the  elementary  autocorrelation  at 
taking  Fourier  transforms, 


OJ  / \ 

mjk  v = 


m0 


ID  / \ 

( V ) 


+27Ti3-  ^(0V0Ej) 


0. 


Thus,  we  have, 


(VIII-2  *1 ) 


where  mo^Cv)  designates  the  spectral  density  in  direction  u . Then,  by 
definition,  omnidirectional  noise  is  composed  of  uncorrelated  contributions  having 
the  same  spectral  density  in  all  directions. 

Thus, 

a.  ra0  u ( v ) is,  in  fact,  independent  of  id  , and  may  be  written  m0  ( v ) , and 

b.  the  elementary  crosscorrelations  and  their  transforms  add  simply 

with  integration  with  respect  to  id  , giving  the  total  crosscorrelation  . 
Thus  the  total  crosscorrelation  spectrum  at  inputs  and  is 


V”- 


dw 


[(OE.-OE  -u)] 

K J din 


(VIII-25) 


The  scalar  product  of  ( VI I -2  5 ) only  involves  the  distance  Ej  KR  . The 
classical  result,  according  to  which  Yj  R ( v ) is  real  and  even  in  ._v,  is  easily 

justified  on  the  basis  of  this  equation.  The  result  is  that  y j k ( t ) is  even_m 

for  any  pair  of  points.  This  result  is  physically  obvious,  the  data  being  symmetrical 

with  respect  to  the  median  plane  of  EjEj<  ' APPlyin8  a delay  ^ +1  ^ at  E,1 
or  (-t  ) at  E gives  the  same  result.  Furthermore,  the  crosscorrelation  does 

K 

not  depend  upon  the  distance  E.E^  • 
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VIH-6.  Matched  Filtering  in  Omnidirectional  Noise-Similarity  to  Optimization 
of  Directivity  Factor 


The  second  step  will  take  into  account  the  fact  that  input  E 1 has  been  affected 

by  a "delay"  -2ni^-(0E  . • oj  ) and  E.  by  a "delay"  -2  n i^-(OE^.  • w ) 
e e 


Transposing  Eq.  (III-10)  , where  the  two  preceding  filtering  processes 

play  the  roles  of  R . ( v ) and  R ( v ) and  y . , ( v ) plays  the  part  of  C . .,  ( v ) , 

the  new  crosscorrelation  spectrum,  after  these  "delays"  make  the  signals  identical,  is 


y3k(v) 


2iri£(0Ej-0Ek)  .<o0 


(VIII-26) 


As  a result  (see  (VIII-25)), 


Yjk(v)  = m0 ( v ) 


f — 2 tt  i—  [ 

i c 


~ [(OE  -OEk).(u-u0)] 


dw 


(VIII-27) 


Let  us  compare  the  preceding  equation  to  Eq.  (VIII-6)  . Except  for  a real 

* 

factor  m0(v)  , the  Yjk(v)  are  identical  to  the  V. ^(v)  . These 

Y j k ( v ) are  the  crosscorrelation  spectra  of  the  noise  inputs  for  which  we  propose  to 
find  the  proper  filtering.  As  we  know,  the  solution,  except  for  a real  factor,  is 

YP  = « , (VIII-29) 

where  Y is  the  matrix  of  the  Y j j.  ( v ) and  p is  the  column  matrix  of  the  proper 
filters.  It  may  be  written  in  the  form 

Vp*  = a , (VIII-30) 

which  is  equivalent  to  the  solution  (VI 11-20)  . 

In  conclusion,  proper  matched  filtering  in  omnidirectional  noise  is  that  filtering 
which  optimizes  directivity. 


90 


in  the  same  plane  normal  to  the  reference  direction,  matched  filtering  is  real. 

VIII-7.  Application  to  Two  and  Three  Inputs 

Let  us  apply  the  preceding  results  to  a system  having  two  elements.  We  will  assume 
that  they  are  arranged  along  a line  parallel  to  the  plane  wave  front,  such  that  the  signal 
is  the  same  at  the  two  inputs.  Let  us  use  the  narrow-band  approximation.  In  omni- 
directional noise,  the  power  of  the  interference  noise  is  the  same  at  the  two  inputs. 
Furthermore,  the  crosscorrelation  spectrum  is  real.  The  two  filters  it  j and  tt2 
(see  Eq.  (VII-5)  ) are  real  and  identical.  They  may  be  replaced  by  the  direct 
summation,  which  means  that  there  is  no  efficient  matched  filter  for  two  inputs  in 
omnidirectional  noise.  This  is  physically  obvious  from  the  symmetry  of  the  parts 
played  by  the  two  inputs.  Note  that  we  have  here  a very  exceptional  case,  where  the 
direct  sum  is  a part  of  the  matched  filtering  although  the  noise  inputs  are  not 
independent  (see  Paragraph  V-2);  but  this  is  true  only  for  two  inputs  in  omnidirectional 
noise.  With  two  antenna  elements,  it  may  be  only  a question  of  their  separation 
distance  for  improvement  of  the  signal -to-noise  ratio.  It  is  well  known  that  they 
must  be  placed  with  a separation  such  that  the  crosscorrelation  of  the  two  noise  inputs, 
which  is  a function  of  this  separation,  is  minimized. 
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The  spectral  density  at  an  arbitrary  input  is  4 nm0  ( v ) 
after  the  input  signals  are  made  identical. 


. It  remains  the  same 


Thus,  the  quantity 


4m0(v) 


■ Vv) 


represents  the  normalized  crosscorrelation  of  the  noise  inputs  at  and 
after  the  signals  have  been  made  identical. 

Whenever  the  vectors  OE,  are  perpendicular  to  the  direction  w0f  the 

J 

crosscorrelation  7 j ^ ( v ) is,  except  for  a multiplicative  factor,  reduced  to  Yj  ^ ( v ) 
(Eq.  ( 7111-26  ) ),  that  is,  a real  quantity  that  is  even  in  v but  not  necessarily 
non -negative.  The  signals  are  then  identical  without  the  need  of  introducing  any 
delays,  since  all  the  inputs  are  in  the  same  plane  normal  to  u0  .In  this  case, 
the  matrices  7 and  V are  real  and  symmetric  (hence,  always  Hermitian).  The 
matched  filter  p is  real.  Hence,  in  omnidirectional  noise  and  for  inputs  situated 


-< 


If  a is  this  separation,  calculation  of  Vj^.(v)  from  Eq.  (VIII-6)  gives 


(a)  sin 


2irva 


2 a 

IT  V 


in  the  case  of  two  point  receivers; 

(b)  in  the  case  of  two  parallel  line  elements  of  undefined  length,  we  have 
a planar  situation  where  the  corresponding  expression  is 

T / 2 tt  v a % 

J o 1 > 


In  each  case  there  exists  an  optimum  separation  distance  between  the  two  elements 
(for  a given  frequency) . 

Let  us  call  P ( a ) one  of  the  two  preceding  functions  and  consider  the  case  of 


three  equidistant  elements  lined  up  in  the  order  1,  2, 
The  matrix  of  the  V^. is: 

3. 

’ 1 

F(  a) 

F(2a)" 

F(a) 

1 

F(a) 

F(2a) 

F(  a) 

1 . 

• 

Consequently,  the  filters 

11  1 » n2  » , 

and  77 3 

are 

*1  = 

TT3  = 

[l-F(a)] 

[l-F(  2a )] 

tt2  = [l-F(2a)]  [l+F( 2a  )-2F(a )]  > 
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Finally,  the  weighting  (real  filtering)  to  apply  to  the  two  end  elements  (1  and  3) 
with  reference  to  the  central  element  is 

l-F(a) 

l+F(2a)-2F(a) 

It  is  simple  to  generalize  this  method  to  on  arbitrary  number  of  elements  with 
arbitrary  spacing. 


1 
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CHAPTER  IX 

GENERATION  OF  PROPER  FILTERING  - AUTOADAPTIVE  SYSTEMS  - 
PRACTICAL  STATIONARITY  CONDITIONS 


Summary 

In  this  chapter,  we  will  deal  with  the  technology,  or  at  least  the  principles  of  the 
technology,  of  matched  filtering.  It  is  a question  of  constructing  the  complex 
quantities  representing  the  N filters,  and  then  making  use  of  them  to  do  the  filtering. 

Furthermore,  we  desire  to  ensure  that  the  filters  and  filtering  processes  will 
evolve  in  a continuous  maimer,  adapting  to  slow  variations  in  the  statistical  relation- 
ships of  the  noise  inputs.  The  expression  "slow  variations"  is  defined  by  the 
assumption  of  a "practical  duration  of  stationarity"  which  is  large  compared  to  the 
time  constants  used  and  the  duration  of  the  signal.  The  essential  steps  are: 

a.  formulation  of  the  elements  of  the  correlations  matrix  c 

in  the  form  of  electrical  voltages, 

b.  formulation  of  the  product  of  several  elements,  of 


'j  k 


cofactors  M.,  and  of  the  "filters' 

Jk 


themselves,  and 


c.  a method  of  "filtering"  each  channel,  starting  with  the 
complex  number  representing  its  filter. 

Finally,  the  block  diagrams  for  the  case  of  three  inputs  and  for  two  inputs  are 
presented.  These  diagrams  may  be  used  to  verify  the  permanence  of  the  filtering 
process  in  the  presence  of  noise,  superimposed  on  the  inputs,  which  is  not  correlated 
with  the  initial  noise  inputs. 

XI-1.  Representation  of  a Crosscorrelation  (Narrow  Band)  of  Two  Slowly  Varying  Voltages 

For  the  moment,  we  will  treat  an  aspect  of  matched  filtering  that  is  close  to  the 
technology  without  going  into  the  details  of  the  technology.  That  aspect  is  the  effective 
realization  of  matched  filtering  in  the  special  case  of  the  narrow-band  approximation. 

A filter  it  of  an  NPF,  for  example,  is  made  up  of  a combination  of  c 

J k 

(a  homogeneous  polynomial  of  degree  1 ).  Hence,  first  of  all,  these  C,.  must 

,1  k 

be  "fabricated."  We  know  that  c ( v n ) represents  the  crosscorrelation  of  two 
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narrow-band  noises  b.(t)  and  b^(t)  resulting  from  the  respective  filtering 
B . ( t ) and  E^_  ( t ) by  filters  of  bandwidth  0 centered  at  the  frequency  v 0 
(and  at  -v0  , considering  only  the  case  of  filters  with  real  unit  impulse  responses 
(see  Paragraph  HI-2)). 

We  will  assume  that  the  spectral  band  <6  contains  all  the  signal,  while  remaining 
narrow  enough  such  that  Cj  ^ ( v ) remains  practically  equal  to  j.  ( v 0 ) • Under 

these  conditions  the  crosscorrelation  function  of  b ( t ) and  b (t ) is 

J K 


TJk(,) 


20  |C^(  v°  1 | Cos 


2 IT  V q T 


arg 


Cjk(v< 


(IX-1) 


The  value  of  y . y ( t ) for  t = 0 is 

Yjk(T\  = o = ‘’^R  [Cjk(v°)]  * (IX-2) 

Hence,  y j ^ ( t ) ^ Q is  the  mathematical  expectation  of  the  product  of  the  two 
noises  and,  by  virtue  of  ergodicity,  the  time  average  of  this  product.  Thus,  the 

quantity  (_  j k ' v 0 J may  be  realized,  except  for  the  factor  2 0 , in  the  following 
way  (Fig.  IX-1): 

a.  by  means  of  a multiplier,  produce,  in  the  form  of  an  electrical 

voltage,  the  product  of  the  noises  b , ( t ) and  b ( t ) received 
at  inputs  E j and  E,^  after  prefiltering  0 ; 

b.  take  the  time  average  of  the  multiplier  output  voltage  by  means 

of  a sufficiently  narrow  low-pass  filter  (an  integrator  with  an 
adequate  time  constant),  delivering  a DC  component  plus  a 
fluctuating  component  which  is  small  for  a large  time  constant,*  and 

c.  the  output  voltage  of  the  low  pass  filter  represents,  except  for 

the  residual  fluctuation,  the  time  average  of  the  product  of  the 
two  noise  inputs,  or  R ^ j.  ( v 0 )J* 


<*> 


Fig.  IX-1 


If  we  were  also  able  to  produce  a voltage  representing  I Jc  ( v0)  1 
we  could  represent  the  complex  number  C k ( v 0 ) by  two  electrical  voltages. 

In  order  to  do  this,  let  us  take  the  two  noise  inputs  b . ( t ) and  b ( t ) and  let 

j k 

T = in  Eq.  (IX-1)  . It  is  clear  that 

4 v0 


YJk 


( — — ) = 2*1 
4v0 


( IX-3) 


Thus,  [ j k:  ( v 0 ) I ^s*  except  for  the  factor  2 0 , the  mathematical  expectation 
of  the  time  average  oT  a voltage  obtained  by  taking  the  product  of  b ^ ( t ) with  a 
noise  voltage  b^ ( t ) , itself  obtained  by  delaying  b^  ( t ) by  ; -*• 

Since  we  arc  dealing  with  a narrow  band,  this  time  delay  is  equivalent  to  a phase 
lag  of  w/2  . This  amounts  to  the  same  thing  as  realizing  an  advance  of  tt/4  in  the 
phase  of  b ( t ) and  a lag  of  it/ 4 in  the  phase  of  b^(  t ) . The  two  processes  of 


generating  I ^ ^ v illustrated  in  Fig.  IX -2  are  equivalent.  The  output  voltage, 
except  for  some  fluctuation,  represents  the  value  of  I j*C  ^ . . ( v 0 )1 


IX-2.  Practical  Stationarity  and  Time  Constants  in  an  Autoadaptive  System 


The  values  of  C „ ( v 0 ) mid  of 


C.  ( v o ) may  be  equally  well  represented  by 


^kkv  v 0 ' 

the  mathematical  expection  of  the  time  average  of  the  squares  of  b ( t ) and  b , , ( t ) . 
Although  very  different  techniques  may  be  used  to  form  these  quantities,  their 
generation,  using  methods  compatible  with  the  preceding  figures,  is  represented  in 
Fig.  IX-3. 


u = cos2nFt  u'  = sin2irFt 


Fig.  IX-3.  Representation  of  the  Complex  Numbers  C „ by 
AC  Voltages  with  Frequency  F 
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Finally,  for  two  inputs,  we  have  formed  four  voltages  which  represent 


j 

omit  the  indication  of  frequency  v 0 . 

With  all  the  possible  C . it  is  possible  to  construct  the  filters  tt,  . This 

J K K 

consists  of  obtaining  complex  numbers  (in  the  form  of  voltages)  equal  to  the  . 

Next,  properly  defined  filtering  consists  of  applying  the  filters  to  the  noise  inputs 
b ( t ) and  taking  the  sum  of  the  filter  outputs.  Later  we  will  see  how  to  realize  the 

tt  and  the  filtering  process  itself.  But,  for  the  present,  it  is  interesting  to  note  that 

K 

if  the  noise  inputs  are  strictly  stationary,  the  C ^ may  theoretically  be  obtained 
in  the  form  of  strictly  continuous  voltages  with  the  aid  of  infinitely  narrow  low -pass 
filters.  In  order  to  avoid  confusion,  from  now  on  we  will  call  these  filtering  processes 
integrations . 

In  practice,  the  C . are  obtained  to  within  only  a fluctuation  whose  amplitude 
is  dependent  upon  the  integration  time  constant.  On  the  other  hand,  since  the  noise 
inputs  vary  slowly  with  time  (limited  stationarity) , the  Cj ..  obtained  from  the  above 
methods  are  going  to  evolve  slowly.  The  filters  derived  from  them  will  also 

evolve,  achieving  an  autoadaptation  of  the  system  that  follows  the  slow  variations  of 

the  statistical  relationships  of  the  B . ( t ) . 

J 

Consequently,  if  an  "average  stationarity  time"  can  be  defined  (an  essentially 

practical  concept),  it  will  be  useless  to  integrate  the  outputs  of  multipliers  whose 
time  constants  9 are  greater  than  this  average  time  0 . On  the  other  hand,  if  the 
time  constant  0 is  not  greater  than  the  inverse  of  the  filter  bandwidth  0 , it  is 
ineffectual  in  practice. 

Thus  we  have  imposed  the  following  conditions. 

a.  The  signal  is  narrow  band,  lying  within  the  band  <p  . Thus,  its  time 
duration  T is  greater  than  1/<P  . ~ 


11  The  margin  between  T and  l/<t>  allows  the  possibility  of  several  signal  types 
in  the  case  of  doppler  effect,  for  example.  In  particular,  in  the  case  of  long  signals  with 
a "single  frequency,"  the  output  of  the  proper  filtering  may  be  handled  by  "band  division," 
which  corresponds  to  matched  filtering  for  each  possible  value  of  doppler  shift  (see 
reference  (23)). 


(v0) 


From  now  on  we  shall 


Cjj(v°)>Ckk(vo) 


[c.k(-°)],  and  I [< 
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b.  The  band  0,  although  larger  than  that  of  the  signal,  is  narrow  enough 

for  the  narrow-band  approximation  to  be  valid,  that  is,  such  that  c ( ) are 

jk  ' 

equal  to  Cj  ,r  ( v 0 ) within  the  band. 

c.  The  noise  may  be  considered  as  stationary  for  duration  e s , which  is  very 
large  compared  with  1/0  and  T.  The  integration  time  constant  is  both  large 
compared  with  1/0  and  T and  small  compared  with  e . 


1/0<T<<6<<6 

S (IX-A ) 

Thus,  the  c,,  are  defined  for  durations  on  the  order  of  6 . They  are  con - 

j k s 

sidered  as  being  perfectly  constant  over  the  duration  e , which  is  the  duration  over 
which  they  are  computed  by  an  integration  of  products. 

d.  Since  6 is  large  compared  with  1/0,  the  multiplier  - integrator  system 
may  be  considered  to  be  a correlator  "with  strong  integration"  (reference  7,  Chapter  V). 
The  power  associated  with  the  fluctuation  of  the  output  is  proportional  to  1/00, 
that  is,  to  a very  small  quantity.  This  fluctuation  is  considered  to  be  negligible. 

In  order  to  establish  the  ideas  more  solidly,  taking  inspiration  from  a practical 
case,  let  us  assume  that  0 is  on  the  order  of  100  cps  at  a center  frequency  of  10  kc, 
which  is  certainly  a narrow  band.  The  duration  of  signals  is  on  the  order  of  a second  , 
and  their  bandwidth  is  a few  cycles  per  second.  The  interference  noise  is  stable  - 
let  us  say  "practically  stationary"  - for  durations  of  from  30  to  60  seconds.  We  may 
envisage  integration  time  constants  on  the  order  of  5 to  10  seconds. 

IX-3.  Representation  of  a Crosscorrelation  by  an  AC  Voltage 

Having  thus  defined  the  practical  conditions  of  stationarity  which  allow  the  use  of 

the  preceding  theoretical  developments,  let  us  return  to  the  c represented  by 

J 

constant,  or  more  precisely,  slowly  varying  voltages.  We  will  indicate  briefly,  in 
the  rest  of  this  chapter,  how  the  principles  of  matched  filtering  with  N inputs  may  be 
implemented. 

We  will  see  later  that  in  order  to  apply  filtering  irk  to  noise  ( t ) it  is 
particularly  convenient  to  express  the  complex  number  n in  the  form  of  an  AC 

K 

voltage  of  fixed  frequency,  whose  amplitude  and  phase  are  equal  to  the  modulus  and 

argument,  respectively,  of  Moreover,  it  is  quite  natural  to  represent  any 

complex  numbers  by  AC  voltages.  Thus,  although  we  will  "fabricate"  the  C,,. 

J K 
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in  the  form  of  DC  voltages,  we  are  going  to  convert  them  to  AC  voltages  with  an 
arbitrary  frequency  P . 

Let  us  start  with  a voltage  source  of  frequency  F , which  will  serve  as  a phase 
reference  at  that  frequency.  Its  amplitude  is  normalized  to  the  value  1, 

or 

u = cos  2tt  P t 


Take  the  product  of  u with  the  DC  voltage  representing  C ^ j (a  simple  modulator 
is  adequate).  We  then  obtain  the  voltage  c . . cos  2 it  F t with  zero  phase,  which 
will  henceforth  represent  the  complex  number  Cjj  - here  it  is  real  and  non-negative. 
Similarly,  we  will  have  the  voltage  c cos  2 tt  F t. 

xCrC 

Take  the  product  of  u with  the  DC  voltage  representing  R 
(see  Fig.  IX -3)  y^  = R TCj^  J cos  2tt  F t. 

From  voltage  u , derive  a voltage  u ' of  the  same  (normalized)  amplitude  but  given  a 
phase  lag  of  tt/2  : 

u'  = sin  2tt  F t . 


M- 


We  obtain 


Take  the  product  of  u'  with  the  DC  voltage  which  represents  I 


V 


2 


sin  2n  F t 


We  have 


( IX-6 ) 


We  find  that  the  difference  between  the  voltages, 


V1  - V2 ! 


represents  in  ampli- 


tude (modulus)  and  in  phase  (argument)  the  complex  number  C ..  . The  sum  y 

* # J K 1 

and  y represents  r,  • Hereafter,  the  C.u  and  C,  are  represented  by 

2 J K.  J K*  J K 

voltages  with  frequency  F and  appropriate  amplitude  and  phase. 


EX -4.  Principles  of  Producing  a Product  of  Complex  Numbers  in  the  Form  of  an 
AC  Voltage 


The  following  operation  consists  of  realizing  products  of  the  C . k . The  products 
are  also  complex  numbers  representable  by  AC  voltages.  Let  a ^ , for  example, 
be  represented  by 


| | cos[2TrFt  + (arg  a^] 
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and  a 2 by 

| a0  | cos  [2irFt  + (arg  a0)] 

Using  a multiplier,  take  the  product  of  the  two  voltages.  Two  spectral  components 
are  obtained,  one  at  frequency  2F, 

W = | a^|  • | a.,  | cos  j^nFt  + (arg  a-^)  + (arg  * 

and  the  other  at  frequency  zero.  We  may  eliminate  the  second  component  by  a high- 
pass  filter  or,  even  better,  isolate  the  first  by  a band-pass  filtei'  wide  enough  to  let 
the  slow  variations  in  |a^|  • |a0  j and  in  (arg  a^  + arg  aR  ) pass  through. 

Now  the  voltage  W represents  the  product  of  two  complex  numbers  a ^ and  an  • 

Thus  the  product  of  two  complex  numbers  Cj  R is  represented  by  an  AC  voltage  at 
frequency  2F  . The  phase  reference  voltage  for  these  products  is  obviously  the 
voltage  with  frequency  2F  , obtained  by  squaring  voltage  u . 

We  have  at  our  disposal,  at  present,  monomials  of  degree  N-l  that  enter  into 
the  composition  of  a "R.  We  may  omit,  in  this  process,  the  intermediate  computation 
of  cofactors  Mj , which  are  sums  and  differences  of  the  above  monomials  (irR  being 
itself  the  sum  of  the  cofactors  of  the  1 column  of  A)  . The  phase  reference  at 
frequency  (N-l)  F = f is  handled  in  the  same  way. 

IX-5.  Narrow-band  Linear  Filtering  Using  Multiplication  By  a Voltage  Representing 
The  Filter 

The  final  step  consists  of  "filtering"  the  1 channel  with  tir,  that  is,  changing 
the  amplitude  and  phase  of  the  voltage  at  frequency  v o which  appears  at  each  channel 

a.  by  multiplying  its  amplitude  by  /nR/ » :Lr|d 

b.  by  increasing  its  phase  by 

Thus,  it k is  represented  by  the  "filtering  voltage"  4> : 

$ = |nk|cos(?nft  + arg  % ) , (IX_8) 

letting  f = (N-l)  F. 
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The  input  voltage  (to  be  filtered)  is 

X = acos(2TTVQt  + X)  , 

and  we  wish  to  obtain 

y = ■ a | itjJ  cos  ( 2™  t + X + arg  ir^)  . 


( IX-9 ) 


(IX-10) 


Take  the  product  of  the  voltage  0 and  the  voltage  x and  filter  the  component  at 
frequency  ( v 0 + f ) . 

z = | it  k | c o s ( 2 it  ( f'  + V Q)t  + x + arg  nR)  . (IX-11) 

The  above  voltage  has  the  amplitude  and  phase  of  y , but  it  is  not  at  frequency 
\»0  • In  certain  cases,  this  fact  may  not  be  troublesome.  Thus,  the  output  of  the 
matched  filter  is  obtained  in  the  form  of  a narrow-band  voltage  with  center  frequ^M 
f + v o . If,  however,  we  have  to  recover  the  frequency  «o  > it  may  be  done  with 
the  aid  of  the  "phase  reference"  voltage  at  frequency  f . 

r = cos2iift  . (IX-12) 

Taking  the  product  of  z and  r and  filtering  for  the  component  at  vo  (demodulation 
with  r ) , the  voltage  y itself  may  be  obtained  (except  for  an  unimportant  factor  of 
1/2). 

The  only  thing  remaining  is  to  take  the  sum  of  the  outputs  of  all  the  channels  in 
order  to  obtain  normalized  proper  filtering.  Thus  we  have  given  a very  brief  glance 
at  realization  of  an  NPF.  Numerous  variations  are  possible  in  the  technology  itself, 
but  the  basic  circuit  remains  as  a multiplier  of  two  narrow-band  voltages.  Simple 
modulators  may  be  used  to  realize  certain  of  the  multiplications  (when  one  of  the 
factors  has  a normalized  amplitude).  In  the  majority  of  cases,  we  need  either  an 
analog  multiplier  or  a digital  multiplier  of  quantized  voltages  accurate  enough  to 
minimize  inaccuracy  in  the  average  value  of  the  products  (for  example,  geometric 
quantization,  see  reference  (24)) . 

In  order  to  realize  a CSPF,  we  have  seen  that  it  is  necessary  to  divide  the  result 


of  the  preceding  NPF  by  r = In,..  To  form  presents  little  difficulty,  and 

it  will  be  obtained  on  the  form  of  a voltage  having  a frequency  f and  zero  phase  ( p 
is  real  and  non-negative) . 


i 
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Of  itself,  the  amplitude  of  this  voltage  is  of  interest  to  us  and  will  be  determined 
by  means  of  an  envelope  detector  in  the  form  of  a constant  or  slowly  varying  voltage, 
always  of  the  same  sign.  Subsequently,  the  output  from  the  NPF  must  be  divided  by 
this  voltage. 

This  operation  of  division  by  a voltage  with  a constant  sign  is  technologically 
possible,  although  a little  more  difficult  than  multiplication.  Thus  the  CSPF  is 
realized. 

IX-6.  Flow  Charts  for  Two-and  Three-Input  Systems 

Figure  IX-4  shows  the  block  diagram  for  an  NPF  in  the  case  of  3 inputs,  or  more 
precisely,  one  of  the  possible  block  diagrams.  The  complexity  of  this  equipment 
rests  primarily  with  the  large  number  of  analog  multipliers  that  it  requires.  Within 
the  present  state  of  the  art,  it  appears  difficult  to  handle  more  than  4 or  5 inputs. 

For  this  reason,  it  is  natural  to  turn  toward  less  perfect  but  more  easily  realized 
solutions,  such  as  that  of  connecting  the  CSPF's  together,  the  principles  of  which  are 
discussed  in  Paragraph  VI-9  and  Paragraph  VI-10. 

The  theory  of  a system  with  two  inputs  is  treated  in  Chapter  VII.  Figure  IX-5 
represents  a possible  realization  of  the  corresponding  circuit.  This  figure  is  similar 
to  Fig.  IX-4.  Six  analog  multipliers  are  required  to  realize  a 2-input  NPF,  and  a 
supplementary  divider  for  an  CSPF  (the  squares  may  be  realized  in  a number  of  ways) . 
The  other  elements  (modulators,  pass-band  filters,  integrators)  are  classical. 

Figure  EX-fi  shows  a simplified  circuit  useful  only  in  the  case  of  two  inputs  having 

the  same  spectral  density.  As  we  know,  the  filters  it  , and  n are  complex  con- 
12 

jugates.  Equation  VII-5  applies  to  this  case,  which  is  characterized  by 


C 1 1 = 0^2  = C , from  which  q = 1 and  C ^ = C 9 
and,  in  particular, 


(IX-13) 


12 

Recall  that  tt  , and  it  are  the  system  functions  (complex  gains)  of  the 
filters  for  the  frequency  v 0 . 
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PROPER  FILTERING-RETURN 
TO  FREQUENCY  vn 


Fig.  IX-4B.  Block  Diagram  of  Filters 


RETURN  TO  FREQUENCY 


Fig.  IX -5.  Two-Input  System  — General  Case 


OUTPUT  NPF 


Fig.  EX-G.  Two-Input  - Equ.nl  Power  Case 


r = 2 c (i-R( e) ) 


ttx  = c (l-e)  = £ rp 


7T.,  = c (l-e*)  = i r +i  ci(  e)  , 

c.  c.  p 


and  it  may  be  seen  that  the  real  part  of  both  n ^ and  u , is  nothing  more  than  r ^ > 
tliat  is,  half  the  spectral  density  of  the  difference.  This  justifies  the  structure  of 
the  circuit  of  Fig.  K-(>,  which  requires  only  4 analog  multipliers. 

IX -7 . Verification  of  one  Property  of  the  NPF  in  the  Preceding  Diagrams 

It  has  been  demonstrated  in  Paragraph  VI-3  that  an  additional  noise  input  n ( t ) , 
uncorrelated  with  B^(t)  and  Bn(t)  and  superimposed  upon  them  changes  the 
behavior  of  neither  the  NPF  nor  the  CSPF.  This  property  may  be  verified  in  the 
sample  circuits  of  Figs.  IX -4,  IX-5,  and  IX-G.  Its  verification  is  particularly  simple 
in  Fig.  L\-(i.  First  of  all,  it  is  obvious  that  everything  that  is  common  to  the  two 
inputs  - the  noise  n ( t ) or  the  contingent  signal  - disappears  in  the  difference  of 
the  two  and,  because  of  this  difference,  does  not  influence  the  process  of  synthesizing 
it  i and  tt  , . 

C. 

It  may  be  proven  that  the  noise  n(t ) disappears  in  the  multiplication-integration 
of  two  inputs  which  differ  in  phase  by  tt/2  . In  fact,  the  two  narrow-band  noises 
b ^ ( t ) mid  b ( t ) may  be  written,  making  more  evident  their  low  frequency  com- 
ponents (whose  bandwidth  is  on  the  order  of  <t>  ) , 


b^(t)  = ( t ) cos  2 it  v p t + N^(t)  sin  2tt  vq  t , 

b -,  ( t ) = M 0(t)  cos  2tr  vq  t + N 0 ( t ) sin  2 n v0  t , 


(iX-lb) 


;uid  the  noise  n ( t ) may  be  written 

n(t)  = p(t)  cos  2 it  v n 


t + q(t)  sin  2i  V5  t , 


where  neither  p ( t ) nor  q ( t ) is  correlated  with  M , ( t ) , M , ( t ) , N , ( t ) , or 

1 d 1 

Np(t ) . Furthermore,  the  mathematical  expectations  of  all  the  low  frequency  com- 
ponents are  zero  (reference  (7),  page  4(>  , Eq,  3-0-31). 


1(1!) 


t 


Following  the  procedure  of  the  diagram  in  Fig.  IX -G,  we  see  that  the  multiplier 
forms  the  product  of 


[Mx(t)  + p(t) 

imd 

• [-sin  2 n v q tj 

+ 

+ q ( t )J 

[m2  ( t ) + p(t)J 

• [cos  2v  v0  tj 

+ 

[M  2 ( t ) 

+ q(  t )] 

[cos  2V  v0  t 


[sin  2 tt  v Q 


( IX-16) 

Integi’ation  done  by  a low -pass  filter  eliminates  components  of  frequency  ;’v0. 
As  a strong  integrator  ( 9>>1/^,  . see  Eq.(IX-4)),  it  forms  the  time  average  - 

or  the  mathematical  expectation  - of  the  low  frequency  terms  in  the  product,  that  is, 
of 


V(t)  = [N1(t)  + q(t)]  [M2(t)  + p(t)]  - [Mx(t)  + p(t)]  [ N2(t)  + q(t)] 


In  V(  t ) , the  term  p ( t ) • q ( t ) is  eliminated  (its  time  average  is  zero,  two 
components  of  the  same  stationary  noise  being  uncorrelated  at  a given  instant  of  lime 
((7),Eq.  3-6-35  ). 

The  time  averages  of  all  products  such  as  q(t)  Mn(t)  or  p(t. ) 
are  zero,  since  they  are  equal  (with  no  correlation)  to  the  product  ol  the  averages  ol 
each  factor. 

Summarizing,  in  the  time  average  of  V ( t ) , at  the  output  of  the  integrator,  only 
quantities  depending  upon  b.(t)  and  b0(t)  remain.  Thus,  the  noise  1 ( t ^ is 
eliminated  in  the  process  ol  synthesizing  l and  n j * ;md  the  predicted  property  is 
therefore  confirmed.  The  preceding  demonstration  still  applies  il  n(t  'i  is  replaced 
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by  a sinusoidal  function  of  frequency  v 0 . It  is  enough  to  assume  in  this  case,  that 

q(t)  = 0 
p(t ) = Do  . 

The  voltage  V ( t ) becomes  V q ( t ) : 

V0(t)  = Nx(t)  [M2(t)  + p0]  - |mx  ( t ) + pol  N2(t)  , 

and  pq  is  eliminated  from  the  time  average  of  V0  (t ) because  the  averages  of  N^(t ) 
and  Np  ( t ) are  zero. 

If  a sinusoidal  function  of  unlimited  duration  - having  the  characteristic  of 
stationarity  - does  not  change  the  value  of  77  ^ and  71 2 , the  same  will  be  true,  certainly, 
for  the  narrow-band  signal  S ( t ) (narrower  than  <t  ) . In  fact,  this  signal  is  a 
sinusoidal  function  with  frequency  v and  slowly  varying  amplitude . Furthermore, 
its  duration  is  small  compared  with  the  integration  time  constant.  For  these  two 
'easons,  the  signal,  identical  at  the  two  inputs,  does  not  depend  upon  the  multiplier- 
integrator  output,  nor,  consequently,  upon  the  values  of  tt^  and  ir0.  Similar 
verification  of  this  property  may  be  made  using  Figs.  IX -5  and  IX -4. 

The  I fCj  are  not  altered  by  the  noise  n ( t ) or  by  a sinusoidal  time  function 
of  frequency  v ^ . The  R J^C  ^ , J are  all  increased  by  the  same  quantity,  which  sub- 
sequently vanishes  when  the  difference  is  taken. 


CONCLUSIONS 

The  general  area  covered  by  this  study  is  an  investigation  of  the  best  way  of 
detecting  the  presence  - and  only  the  presence  - of  weak  signals  in  noise  interference 
and,  consequently,  of  improving  electromagnetic  and  acoustic  detection  over  a wide 
range  of  processes. 

In  this  area  are  included  two  methods  of  approach  to  this  problem.  The  first  is 
statistical  detection  theory,  which  makes  use  of  the  concept  of  estimation.  The  other 
is  the  examination  of  the  optimization  of  certain  judiciously  chosen  criteria  (signal- 
to-noise  ratio) . 


Ill 
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Although  historically  unrelated,  the  answers  obtained  coincide  for  conditions 
usually  found  in  practice. 

The  Matched  Filter  Theorem  is  one  of  those  optimizations  and  concerns  a particular 
signal  in  a stationary  noise  field.  It  forms  the  point  of  departure  for  this  study. 

We  have  tried  to  modify  this  theorem  - which  handles  only  a single  signal  and  a 
single  noise  - to  include  the  multiple  aspects  of  detection  by  an  antenna  with  several 
receiving  elements.  Then  we  have  examined  most  of  the  forms  taken  by  the  solution 
for  different  statistical  relationships  presented  by  the  noise  inputs.  We  wished  to 
stress  the  very  general  character  of  this  treatment  by  showing  how  it  includes  and 
extends  beyond  the  concept  of  directivity,  which  is  generally  used  for  the  same  purpose. 
Finally,  we  have  pursued  the  work  to  the  point  of  block  diagrams  for  the  realization 
of  such  systems,  in  the  particular  case  of  narrow-band  signals,  and  from  the  point  of 
ot  an  automatic  adaptation,  with  gradual  changes  in  a quasi-stationary  noisefield. 
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